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THE  INFLUENCE  OF  MAGMA  SUPPLY  AND  ERUPTIVE  PROCESSES  ON 
AXIAL  MORPHOLOGY,  CRUSTAL  CONSTRUCTION  AND 
MAGMA  CHAMBERS 

Emilie  Ernestine  Ebba  Hooft 

Submitted  to  the  Department  of  Earth,  Atmospheric,  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
and 

the  Department  of  Geology  and  Geophysics,  Woods  Hole  Oceanographic  Institution, 
December  1996  in  partial  fulfillment  of  the  requirements  for  the  degree  of 

Doctor  of  Philosophy 

ABSTRACT 

Two-thirds  of  the  surface  of  the  Earth  is  created  at  mid-ocean  ridges  where  magmas 
rise  from  the  mantle  and  cool  to  form  the  oceanic  crust.  The  objective  of  this  Thesis  is  to 
examine  the  influence  of  magma  supply  and  eruptive  processes  on  axial  morphology, 
crustal  construction,  and  the  properties  of  crustal  magma  chambers  at  intermediate  and  fast 
spreading  ridges.  Variations  in  magma  supply  on  time  scales  of  -100  Kyr  generate  along- 
axis  changes  in  crustal  thickness  and  temperature.  Magma  sill  properties  and  hydrothermal 
activity  are  closely  linked  to  spreading  events  which  occur  on  much  shorter  time  scales  (ca 
10-100  yr)  than  the  longer-term  variations  in  magma  supply  reflected  in  along-axis  changes 
in  ridge  morphology.  The  seismically  constrained  depths  of  ridge  crest  magma  sills  (>l-2 
km)  are  considerably  deeper  than  the  level  of  neutral  buoyancy  (100-400  m).  The  apparent 
inverse  relationship  between  magma  sill  depth  and  spreading  rate  suggests  that  a  thermally 
controlled  permeability  boundary,  such  as  the  solidus  horizon,  controls  the  depth  at  which 
magma  ponds  beneath  mid-ocean  ridges.  Recent  tiiermo-mechanical  models  predict  that,  at 
intermediate  spreading  rates,  rift  valley  and  magma  sill  formation  are  sensitive  to  small 
changes  in  cmstal  thickness  and  mantle  temperature.  Analysis  of  gravity  at  an  intermediate 
spreading  ridge  shows  that  small  differences  in  cmstal  thickness  (300-700  m)  and  mantle 
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temperature  (10-15°C)  are  indeed  sufficient  to  produce  major  changes  in  lithospheric 
strength  and  axial  morphology.  A  stochastic  model  for  the  emplacement  of  dikes  and  lava 
flows  with  a  bimodal  distribution  of  lava  flows  is  required  to  satisfy  geological  and 
geophysical  constraints  on  the  construction  of  the  extrusive  section.  Most  dikes  are 
intruded  within  a  narrow  zone  at  the  ridge  axis.  Short  flows  build  up  approximately  half 
the  cxtmsive  volume.  Occasional  flows  that  pond  at  a  considerable  distance  off-axis  build 
up  the  remainder  of  the  extrusive  section.  This  Thesis  underlines  the  importance  of 
emption  dynamics  in  the  emplacement  of  the  uppermost  volcanic  layer  of  the  cmst  and  of 
the  crustal  thermal  structure  in  controlling  local  variations  in  magma  sill  depth  and  ridge 
morphology. 

Thesis  supervisor:  Robert  S.  Detrick 

Title:  Senior  Scientist,  Woods  Hole  Oceanographic  Institution 
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Chapter  1 

Magma  Chambers,  Crustal  Construction,  and 

Axial  Morphology 
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Evolution  of  Magma  Chamber  Models 


Two-thirds  of  the  surface  of  the  Earth  is  created  at  mid-ocean  ridges  where  magmas 
generated  within  the  mantle  rise  and  cool  to  form  the  oceanic  cmst  The  separation  of  two 
lithospheric  plates  at  mid-ocean  ridges  results  in  upwelling  of  the  underlying  mantle  which 
undergoes  decompression  melting  as  it  rises.  This  melt  is  focused  into  a  very  narrow 
zone,  on  the  order  of  a  few  kilometers  wide,  where  lava  is  emplaced  either  as  dikes,  sills 
and  extrusive  flows,  or  cools  more  slowly  forming  gabbro.  Through  these  processes  new 
oceanic  cmst  is  formed  along  mid-ocean  ridges. 

The  existence  of  cmstal  magma  chambers  at  mid-ocean  ridges  has  been  a  critical 
element  in  most  models  of  oceanic  cmstal  formation  since  the  early  1970s.  Until  the  first 
detailed  seismic  experiments  were  carried  out  on  mid-ocean  ridges  in  the  1980s,  it  was 
widely  believed  that  most  ridge  crests,  even  those  at  relatively  slow  spreading  rates,  were 
associated  with  large  (10-20  km  wide  by  3-4  km  thick),  essentially  molten,  well-mixed 
magma  reservoirs  [e.g.,  Cann,  1974;  Dewey  and  Kidd,  1977;  Casey  and  Karson,  1981; 
Pallister  and  Hopson,  1981;  Nicolas  et  al.,  1988].  These  large  magma  chambers  were 
attractive  for  several  reasons.  They  offer  a  simple  explanation  for  the  lithostratigraphy  of 
oceanic  cmst  inferred  from  seismic  studies  [Raitt,  1963;  Shor  et  al.,  1970]  and  from 
ophiolites  {Geotimes,  1972].  An  upper  layer  of  extmsive  pillow  lavas  and  sheeted  dikes 
could  be  generated  by  periodic  emptions  through  the  top  of  the  magma  chamber.  A  middle 
unit  of  isotropic  gabbro  may  form  by  freezing  of  the  margins  of  the  chamber  due  to  cooling 
from  above.  The  lower  layered  gabbros  would  then  form  by  cooling  and  settling  of 
crystals  on  the  floor  of  the  magma  chamber.  A  large  cmstal  magma  chamber  also  explains 
the  observation  that  most  mid-ocean  ridge  basalts  are  not  in  equilibrium  with  presumed 
mantle  compositions  (e.g.,  Bryan,  1976],  but  require  some  cooling  and  crystal 
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fractionation  at  shallow  crystal  levels  [e.g.,  Stolper  and  Walker,  1980;  Grove  and  Bryan, 
1983]. 


However,  seismic  studies  carried  out  at  both  fast  and  slow  spreading  ridges  over 
the  past  fifteen  years  effectively  preclude  the  existence  of  the  large,  molten  crustal  magma 
bodies  envisioned  in  these  earlier  models  [Solomon  and  Toomey,  1992]  and  have  led  to  a 
completely  new  view  of  the  size  and  longevity  of  cmstal  magma  bodies  at  ridge  crests,  and 
their  role  in  the  formation  of  the  oceanic  crust.  No  seismic  evidence  exists  for  crustal 
magma  bodies  of  any  significant  dimensions  along  slow  spreading  ridges,  like  the  Mid- 
Atlantic  Ridge  [Purdy  and  Detrick,  1986;  Toomey  et  al.,  1988;  Detrick  et  al.,  1990]  and 
crustal  magma  chambers  are  believed  to  be  volumetrically  small,  intmsive  bodies  that  are 
rapidly  frozen  by  hydrothermal  cooling  [Smith  and  Conn,  1992].  Even  at  fast  spreading 
ridges,  like  the  East  Pacific  Rise,  seismic  studies  indicate  that  magma  is  largely  confined  to 
a  narrow  (~l-2  km  wide),  thin  (<1  km),  mostly  molten  magma  lens  or  sill  [Detrick  et  al., 
1987;  Harding  et  al,  1989;  Kent  et  al,  1990;  Hussenoeder  et  al,  1996]  overlying  a 
broader  region  of  the  lower  crust,  which  has  low  seismic  velocities  and  contains  at  most  a 
few  percent  partial  melt  [Toomey  et  al,  1990;  Vera  et  al,  1990;  Caress  et  al,  1992]. 
These  seismic  observations  are  supported  by  gravity  models  of  the  fast  spreading  EPR 
[Madsen  et  al,  1984;  Wang  and  Cochran,  1993;  Magde  et  al,  1995],  which  do  not  require 
a  large  (>10  km  wide),  low  density  crastal  magma  body. 

These  geophysical  data  have  been  interpreted  in  terms  of  a  composite  magma 
chamber  at  fast-spreading  ridges  which  consists  of  a  thin,  narrow,  sill-like  body  overlying 
a  thicker,  wider  crystal  mush  zone  [Sinton  and  Detrick,  1992].  The  manner  in  which  the 
oceanic  crust  is  generated  from  such  a  magmatic  system  is  quite  different  than  for  a  large 
magma  reservoir.  At  fast  spreading  ridges  the  entire  crust  is  envisioned  as  forming  in  the 
magma  sill  with  eruption  through  the  top  of  the  sill  building  the  extmsive  section,  and  sub- 
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solidus  flow  down  and  out  from  the  base  of  the  sill  generating  the  layered  lower  crust 
[Henstock  et  al,  1993;  Quick  and  Denlinger,  1993;  Phipps  Morgan  and  Chen,  1993a]. 
However,  geological  observations  suggest  that  magma  sills  may  also  exist  in  the  lower 
crust  that  freeze  to  generate  the  layered  gabbros  [Kelemen  et  al,  1996].  Short-lived 
magma  bodies  at  slow  spreading  ridge  are  consistent  with  geological  observations  that  the 
oceanic  cmst  is  not  as  well  layered  as  was  previously  thought  [e.g.,  Cannat,  1993].  The 
existence  of  volumetrically  small  crustal  magma  bodies  is  supported  by  the  geochemical 
variability  observed  in  mid-ocean  ridge  basalts  which  is  also  difficult  to  reconcile  with  a 
well-mixed  magma  reservoir  [e.g.,  Langmuir  et  al,  1986;  Mevel,  1987;  Sinton  et  al, 
1991], 

Seismic  observations  of  magma  sill  reflectors  from  the  northern  and  southern  EPR, 
and  sparse  observations  from  intermediate  spreading  ridges  (Valu  Fa  and  Juan  de  Fuca) 
indicate  a  spreading  rate  dependence  for  the  average  depth  of  magma  chamber  reflectors 
[Purdy  et  al,  1992],  These  depth  variations  have  been  modeled  to  depend  on  the  balance 
of  heat  advected  from  below  and  that  removed  by  hydrothermal  cooling  [Henstock  et  al, 
1993;  Phipps  Morgan  and  Chen,  1993a].  The  thermal  balance  at  mid-ocean  ridges  also 
affects  the  rheology  of  the  cmst  [Phipps  Morgan  et  al,  1987;  Chen  and  Morgan,  1990a; 
Chen  and  Morgan,  1990b].  Slow  spreading  ridges  have  an  axial  valley,  while  fast 
spreading  ridges  have  an  axial  high  [Small  and  SandweU,  1989].  The  transition  in  mode  of 
deformation  of  the  plate  boundary  at  mid-ocean  ridges  is  modeled  to  be  very  sensitive  to 
variations  in  cmstal  temperature  and  thickness  [Phipps  Morgan  and  Chen,  1993b]. 

The  new  ridge  crest  magma  chamber  model  which  has  been  developed  over  the  last 
fifteen  years  raises  important,  and  still  unresolved,  questions  about  how  oceanic  crust  is 
formed.  How  is  magma  delivered  from  the  mantle?  Since  the  cmstal  magma  reservoirs  are 
volumetrically  so  small  and  short-lived  at  all  but  the  fastest  spreading  rates,  are  there 
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magma  reservoirs  in  the  underlying  mantle?  How  is  the  lower  crust  formed  from  the 
composite  magma  chamber?  Also,  an  important  problem  is  the  temporal  variability  of  these 
processes.  What  are  the  time  scales  of  input  of  magma  from  the  mantle  and  output  through 
eruption  onto  die  seafloor?  How  are  the  properties  of  the  crustal  magma  bodies  related  to 
variations  in  spreading  rate,  magma  supply,  hydrothermal  activity  or  emptive  cycles?  How 
can  a  thin  narrow  magma  sill  exist  without  being  rapidly  frozen  by  hydrothermal  cooling? 
What  controls  the  width  of  the  neovolcanic  zone  and  how  does  it  control  the  architecture  of 
the  shallow  crust?  What  is  the  size  and  quantity  of  lava  in  individual  eraptions? 

Thesis  Overview 

This  Thesis  discusses  the  importance  of  spreading  rate,  melt  supply,  hydrothermal 
cooling,  eruption  geometry,  and  the  differences  in  timing  of  these  processes  in  constructing 
the  layered  oceanic  crust.  The  objective  of  this  Thesis  is  to  examine  the  influence  of 
magma  supply  and  emptive  processes  on  axial  morphology,  cmstal  constmction,  and  the 
properties  of  cmstal  magma  chambers  at  intermediate  and  fast  spreading  ridges.  The 
magmatic  system  at  mid-ocean  ridges  is  complicated,  and  represents  the  boundary  between 
a  quasi-steady  input  of  magma  from  below  and  the  periodic  intmsion  of  dikes  and  emption 
of  lavas  onto  the  seafloor.  The  local  cmstal  thermal  stmcture  is  important  in  controlling 
magma  sill  properties  on  relatively  short  time  scales  (10s  to  1,000s  of  years)  and  this  is 
superimposed  on  the  longer  term  dependence  of  the  properties  of  the  magmatic  system  on 
spreading  rate  and  magma  supply.  The  diermal  regime  and  composition  of  the  lithosphere 
also  control  the  rheology  thereby  determining  the  axial  morphology  of  mid-ocean  ridges. 

The  lithospheric  stmcture  of  two  ridges  with  contrasting  axial  morphologies  that  arc 
spreading  at  the  same  intermediate  rate  is  investigated  in  Chjq)ter  2  using  gravity  data  from 
the  Juan  de  Fuca  and  northern  Gorda  ridges.  The  Cleft  segment  of  the  Juan  de  Fuca  ridge 
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has  an  axial  high  moq)hology  while  the  northern  segment  of  the  Gorda  ridge  has  a  rift 
valley.  The  difference  in  residual  gravity  anomalies  is  used  to  determine  the  change  in 
crustal  thickness  and  mantle  temperature  between  these  two  ridge  segments.  These 
observational  results  are  used  to  test  thermal  predictions  for  the  dependence  of  ridge 
morphology  on  magma  supply  and  spreading  rate  [e.g.,  Phipps  Morgan  and  Chen, 
1993b]. 

The  relationship  between  seismically  constrained  variations  in  crustal  structure  and 
indicators  of  ridge  magma  supply  is  re-examined  in  Chapter  3.  Observations  over  the  past 
decade  have  lead  to  the  development  of  a  magma  supply  model  for  fast  spreading  ridges. 
Broad,  shallow  ridge  segments  are  assumed  to  reflect  a  high  magma  supply  and  are 
predicted  to  have  larger,  more  robust  magma  chambers,  erupt  more  primitive  basalts,  form 
thicker  crust,  and  have  more  vigorous  hydrothermal  activity  than  sections  of  the  ridge  with 
a  deeper,  narrower  rise  crest.  Multichannel  seismic  reflection  data  from  an  800-km-long 
section  of  the  southern  East  Pacific  Rise  (SEPR)  allow  a  quantitative  examination  of  the 
correlation  between  along-axis  variations  in  cmstal  structure  and  other  indicators  of  magma 
supply  on  a  regional  scale,  and  show  that  there  are  several  scales  of  spatial  and  temporal 
variabiUty  that  modify  the  simple  correlations  predicted  by  the  magma  supply  model.  The 
observed  variation  in  magma  sill  depth  is  compared  to  different  models  proposed  to  explain 
the  ponding  of  magma  in  the  mid-crust. 

In  Chapter  4,  the  common  assumption  that  magma  ponds  at  a  level  of  neutral 
buoyancy  in  the  shallow  crust  [c.f.  Walker,  1989;  Ryan,  1993]  is  tested.  At  the  level  of 
neutral  buoyancy  the  melt  densities  are  equal  to  the  bulk  densities  of  the  surrounding  crust. 
The  depth  of  the  level  of  neutral  buoyancy  is  determined  from  estimates  of  melt  density 
based  on  the  composition  of  erupted  lavas  and  estimates  of  the  crustal  density  structure 
from  the  inversion  of  the  seismic  velocity  structure  of  the  crust  and  in  situ  borehole 


21 


measurements  of  crustal  porosity.  This  estimate  of  the  level  of  neutral  buoyancy  is 
compared  with  seismically  constrained  depths  of  ridge  crest  magma  chambers  to  test  the 
level  of  neutral  buoyancy  model  for  mid-ocean  ridges.  The  results  are  also  discussed  in 
terms  of  the  alternative  hypothesis  that  the  thermal  stmcture  of  the  rise  crest  controls  the 
depth  at  which  melt  ponds  in  the  upper  crust  [Phipps  Morgan  and  Chen,  1993a]. 

A  stochastic  model  for  the  emplacement  of  dikes  and  lava  flows  to  generate  an 
upper  oceanic  crustal  stmcture  similar  to  that  observed  in  seismic  data  from  the  East  Pacific 
Rise  (EPR),  in  ocean  drill  holes,  and  in  ophiolites  is  explored  in  Chapter  5.  Wide-angle 
seismic  reflections  from  the  steep  velocity  gradient  at  the  base  of  seismic  layer  2A  are 
interpreted  to  be  the  extmsive/sheeted  dike  contact.  Seismic  data  from  the  northern  and 
southern  EPR  constrain  the  on-axis  extmsive  layer  thickness,  the  distance  over  which  the 
thickening  of  the  extmsive  layer  occurs,  and  its  off-axis  thickness.  Ophiolites  and  ocean 
drill  holes  (DSDP  Hole  504B)  provide  additional  estimates  of  the  thickness  of  the  extmsive 
layer  and  constrain  the  thickness  of  the  transition  region  from  extmsives  to  sheeted  dikes. 
The  predictions  of  this  model  are  compared  with  measured  seafloor  lava  ages  [Goldstein  et 
al,  1994]  and  with  the  shallow  velocity  stmcture  determined  on  the  northern  EPR 
[Christeson  etal,  1992]. 

Chapter  6  summarizes  the  main  conclusions  of  Chapters  2  through  5  and  discusses 
some  questions  about  the  volcanic  processes  constmcting  the  cmst  and  the  factors 
controlling  the  crustal  rheology  of  mid-ocean  ridges  that  remain  to  be  addressed  or  are 
raised  by  this  work. 

Chapter  2  was  published  in  the  Journal  of  Geophysical  Research,  1995.  Chapter  3 
was  submitted  to  the  Journal  of  Geophysical  Research,  Sept.  1996.  Chapter  4  has  been 
published  in  Geophysical  Research  Letters,  1993.  Chapter  5  has  twen  published  in  Earth 
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and  Planetary  Science  Letters,  1996.  My  co-author  on  Chapters  2  and  4  is  Robert  Detrick, 
on  Chapter  3  my  co-authors  are  Robert  Detrick  and  Graham  Kent,  and  on  Chapter  5,  Hans 
Schouten  and  Robert  Detrick.  I  have  been  supported  in  this  work  by  NSF  grants:  OCE- 
9296017,  OCE-9402033,  and  OCE-80203300,  and  by  the  Joint  Program  Education 
Office. 
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Abstract  We  analyze  gravity  data  over  the  Juan  de  Fuca  and  northern  Gorda  Ridges  to 
understand  the  lithospheric  structure  of  two  ridges  with  contrasting  axial  morphologies  spreading 
at  the  same  intermediate  rate  (28  mm/yr  half  rate).  The  Cleft  Segment,  at  the  south  end  of  the 
Juan  de  Fuca  Ridge,  has  an  axial  high  morphology  while  the  northern  segment  of  the  Gorda  Ridge 
has  a  rift  valley.  Residual  mantle  Bouguer  anomalies  (RMBA)  on  the  northern  Gorxia  Ridge  are 
elevated  relative  to  the  Cleft  Segment  by  10-20  mGal,  indicating  thinner  crust  and/or  a  colder 
mantle.  The  minimum  value  (-50  mGal)  of  the  RMBA  along  the  Juan  de  Fuca  Ridge  is  over 
Axial  Seamount  and  gradually  increases  south  toward  the  Blanco  Transform.  The  observed 
RMBA  are  interpreted  to  result  from  along  axis  variations  in  crustal  thickness  and  mantle  density, 
both  of  which  are  controlled  by  temperatures  in  the  upper  mantle  where  decompression  melting 
occurs.  We  estimate  that  mande  temperatures  are  elevated  by  30°-40°C  beneath  Axial  Seamount, 
resulting  in  an  excess  crustal  thickness  of  -1.7  km.  The  Cleft  Segment  is  associated  with  crust 
that  is  estimated  to  be  only  300-700  m  thicker,  and  mantle  temperatures  are  only  10°-15°C  higher 
th^  beneath  the  northern  Gorda  Ridge.  However,  even  these  small  differences  in  crustal 
thickness  and  mande  temperature  appear  to  be  sufficient  to  produce  a  major  change  in  lithospheric 
strength  and  axial  morphology.  These  results  are  consistent  with  the  predicted  sensitivity  of 
recent  thermo-mechanical  models  for  rift  valley  formation  to  small  changes  in  crustal  thickness 
and  mantle  temperature  at  intermediate  spreading  rates.  We  attribute  the  systematic  differences  in 
axial  morphology,  crustal  thickness,  mantle  temperature,  and  lithospheric  strength  along  the  Juan 
de  Fuca/Gorda  ridge  system  to  the  presence  of  the  Cobb  thermal  anomaly  at  Axial  Seamount. 

The  hotter  mantle  beneath  the  Juan  de  Fuca  Ridge  results  in  greater  amounts  of  decompression 
melting  and  the  formation  of  a  thicker  crust  and  a  thinner,  weaker  lithosphere  than  along  the 
Gorda  Ridge. 


Introduction 

One  of  the  first  order  features  of  mid-ocean  ridges  is  the 
difference  in  axial  morphology  observed  between  fast  and  slow 
spreading  ridges.  Spreading  centers  with  half  spreading  rates  of 
less  than  15-20  mm/yr  are  typically  associated  with  a  1-2  km 
deep,  15-30  km  wide  rift  valley,  while  ridges  witii  half  spreading 
rates  faster  than  about  30-35  mm/yr  are  associated  with  a  2-10 
km  wide,  200-300  m  axial  high  [Macdonald,  1982].  The 
transition  from  rift  valley  to  axial  high  appears  to  be  relatively 
abrupt  and  occurs  at  half  spreading  rates  of  20-30  nim/yr  [Small 
and  Sandwell,  1989].  However,  spreading  rate  does  not  appear  to 
be  the  only  controlling  variable  on  ridge  crest  morphology.  For 
example,  the  Reykjanes  Ridge  has  an  axial  high  despite  a  half 
spreading  rate  of  only  10  mm/yr  [Searle  and  Laughton,  1981]. 
Along  the  Southeast  Indian  Ridge  the  axial  morphology  changes 
from  a  rift  valley  in  the  Australian-Antarctic  Discordance  to  an 
axial  high  east  and  west  of  this  zone  even  though  this  entire 
section  of  ridge  crest  has  a  relatively  uniform  half  spreading  rate 
of  30-35  mm/yr  [Sempere  et  al.,  1991].  Along-axis  changes  in 
ridge  crest  morphology  can  even  occur  at  the  segment  scale. 
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Along  the  southern  Mid-Atlantic  Ridge  several  cases  have  been 
reported  where  the  ridge  crest  is  absent  near  the  midpoint  of  a 
segment  while  showing  a  clear  rift  valley  near  the  ends  of  a 
segment  [Kuo  and  Forsyth,  1988;  Blackman  and  Forsyth,  1991; 
Foxetal,  1991]. 

The  strength  of  the  lithosphere  determines  how  it  deforms 
under  extensional  forces  and  thus  controls  whether  the  ridge  crest 
forms  an  axial  high  (weak  lithosphere  is  isostatically 
compensated)  or  a  rift  valley  (stronger  lithosphere  deforms  by 
brittle-plastic  necking)  [Sleep,  1969;  Lauchenbruch,  1973; 
Tapponnier  and  Francheteau,  1978;  Sleep  and  Rosendahl,  1979; 
Phipps  Morgan  et  al.,  1987;  Chen  and  Morgan,  1990].  The 
rheological  properties  of  the  lithosphere  depend  on  both 
temperature  and  composition  (crustal  rocks  are  weaker  than 
mantle  rocks  at  the  same  temperature  [c.g..  Brace  and  Kohlstedt, 
1980;  Kirby,  1983]).  The  thermal  structme  of  the  oceanic 
lithosphere  is  determined  by  the  advection  of  heat  by  magma 
emplacement  and  conductive  and  hydrothermal  cooling.  An 
increased  magma  supply  at  a  fixed  spreading  rate  increases  the 
heat  supply  and  creates  a  greater  thickness  of  weaker  crustal 
rocks  and  a  hotter,  weaker  lithosphere.  At  intermediate  spreading 
rates,  ridge  crest  morphology  is  particularly  sensitive  to  changes 
in  the  thermal  structure  and  magma  supply  [Phipps  Morgan  and 
Chen,  1993b]. 
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The  Juan  dc  Fuca  and  Gorda  Ridges  arc  a  classic  example  of 
the  difference  in  axial  morphology  that  is  found  at  intermediate 
spreading  rates.  Although  both  ridges  are  spreading  at  about  the 
same  rate  (~28  mm/yr  half  rate),  the  northern  Gorda  Ridge  is 
characterized  by  an  axial  valley  (400-600  m  deep)  typical  of  slow 
spreading  ridges,  while  the  southern  Juan  de  Fuca  Ridge  is 
generally  associated  with  an  axial  high  (up  to  300  m  high) 
characteristic  of  fast  spreading  ridges.  In  this  paper  we  compile 
gravity  and  bathymetry  data  over  these  two  ridges  and  interpret 
the  residual  gravity  anomaly  in  terms  of  mantle  temperature 
anomalies  and  crustal  thickness  variations.  By  correlating 
variations  in  axial  morphology  with  differences  in  lithospheric 
structure,  we  test  current  ideas  on  the  origin  of  rift  valleys  and 
axial  highs  and  show  that  relatively  small  changes  in  crustal 
thickness  and  mantle  temperature  are  sufficient  at  intermediate 
spreading  rates  to  result  in  large  differences  in  the  morphology  of 
the  ridge  axis. 


Tectonic  Setting 

The  Juan  de  Fuca/Gorda  spreading  system  has  been  intensely 
studied  over  the  past  fifteen  years  (see  special  issues  on  the  Juan 
de  Fuca  Ridge  published  in  the  Journal  of  Geophysical  Research 
in  1987,  1990,  and  1993).  The  ridge  crest  has  been  subdivided 
into  a  number  of  distinct  spreading  segments  bas^  on  detailed 
bathymetric  and  multibeam  mapping  {Delaney  et  ai,  1981; 
Malahoffet  ai,  1982;  Embley  et  al.,  1983;  Karsten  et  at,  1986J 
On  the  Juan  de  Fuca  Ridge  these  segments  are  called  (from  north 
to  the  south)  the  West  Valley,  Endeavour,  Northern  Symmetric 
(or  Cobb),  Axial  (or  Coaxial),  Vance  (or  Southern  Symmetric), 
and  Cleft  Segments  (Figure  1).  The  Gorda  Ridge  consists  of 
three  main  segments,  the  northern,  central,  and  southern  Gorda 
Ridge  (also  known  as  the  Escanaba  Trough).  The  Juan  de  Fuca 
Ridge  is  bounded  to  the  north  by  a  triple  junction  formed  by  the 
ridge,  the  left  lateral  Nootka  Fault,  and  the  I40-km-long  Sovanco 
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Figure  1.  Bathymetry  map  over  the  Gorda  Ridge,  the  Blanco  Fracture  Zone  and  the  Juan  de  Fuca  Ridge.  Grid 
spacing  is  1  km;  contour  interval  is  400  m.  Both  ridges  are  spreading  at  the  same  intermediate  spreading  rate  (~28 
mm/yr).  The  data  included  in  the  bathymetry  compilation  were  shipboard  bathymetry  extracted  from  the  National 
Geophysical  Data  Center  (NGDC)  and  multibeam  bathymetry  from  the  Lamont-Doherty  Earth  Observatory  RIDGE 
Multibeam  Data  Synthesis.  Between  ship  tracks  the  bathymetry  was  supplemented  with  DBDB5  data  (Digital 
Bathymetric  Data  Base  5.)  Where  multi  beam  bathymetry  existed,  the  NGDC  bathymetry  was  not  used.  The 
accuracy  of  this  map  varies  widely,  from  -50  m  over  the  ridge  crest  which  is  entirely  mapped  with  multibcam  data 
and  de^:ading  into  regions  that  were  filled  using  the  DBDB5  data.  The  location  of  the  ridge  axis  is  shown  by  the 
black  lines  as  are  the  two  across-axis  profiles  shown  in  Figure  2. 
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Fracture  Tjonc  [Atwater^  1970;  Hyndman  et  aL,  1979].  To  the 
south  the  Blanco  Fracture  Zone  offsets  the  spreading  axis  400  km 
from  the  Juan  de  Fuca  to  the  Gorda  Ridge. 

The  Juan  de  Fuca  plate  is  being  flooded  by  turbidites  derived 
from  the  North  American  continental  margin  to  the  east  [McCoy, 
1983],  and  their  low  density  affects  the  gravity  signature. 
Sediment  thicknesses  reach  -500  m,  150  km  to  the  east  of  the 
ridge  axis  versus  -250  m  at  the  same  distance  on  the  western  side 
of  the  ridge  [Ludwig  and  Houtz,  1976].  The  turbidites  also  flow, 
via  the  Cascadia  channel,  through  a  small  puli-apart  basin  in  the 
Blanco  Fracture  Zone  (the  Cascadia  Depression)  onto  the  Tufts 
Abyssal  Plain  to  the  west  of  the  Gorda  Ridge  [Embley  and 
Wilson,  1992].  Within  this  region,  sediment  thicknesses  range 
from  -400  m  in  the  NW  comer  to  -100  m  in  the  SW.  East  of  the 
Gorda  Ridge  sediments  also  thicken  very  rapidly  reaching  1-2  km 
over  the  subduction  trench  at  the  continental  margin. 

The  Juan  de  Fuca/Gorda  area  is  part  of  a  classic  magnetic 
anomaly  study  region  [Raff  and  Mason,  1961;  Vine  and 
Matthews,  1963],  and  the  spreading  rates  and  tectonic  history  of 
this  area  have  been  well  characterized  [Hey  and  Wilson,  1982; 
Wilson  and  Hey,  1984;  Wilson,  1988,  1993].  The  tectonic 
evolution  of  the  Juan  de  Fuca  Ridge  over  the  past  18  Ma  has  been 
dominated  by  a  southward  propagating  rift  that  intersected  the 
Blanco  Fracture  Zone  at  -1.5  Ma  [Wilson,  1988].  Within  the 
Brubnes  the  Vance  Segment  has  been  accreting  at  a  nearly 
symmetric  rate  of  28.5  mm/yr,  while  the  Cleft  Segment  is 
spreading  faster  to  the  west  (28.5  mm/yr)  than  to  the  east  (23.2 
mm/yr)  [Tivey,  1994].  Early  interpretation  of  magnetic 
anomalies  across  the  Gorda  Ridge  revealed  half  spreading  rates 


Cleft  Segment  Bathymetry 


Figure  2,  Bathymetry  profiles  across  the  central  portion  of  the 
axis  of  the  Cleft  Segment  of  the  Juan  de  Fuca  Ridge  and  the 
northern  segment  of  the  Gorda  ridge.  The  location  of  these 
profiles  is  shown  by  black  lines  in  Figure  1,  The  Cleft  Segment 
has  an  axial  high  of  -300  m  with  a  small  graben  at  the  summit, 
typical  of  the  morphology  of  fast  spreading  ridges.  In  contrast, 
the  northern  segment  of  the  (jorda  Ridge  exhibits  a  rift  valley  of 
-1  km  depth,  similar  to  the  morphology  of  the  slow  spreading 
Mid-Atlantic  Ridge. 


varying  from  30  mm/yr  on  the  northern  section  of  the  ridge  to  10 
mm/yr  on  the  southern  Gorda  Ridge  [Atwater  and  Mudie,  1973; 
Riddikough,  1980].  Kinematic  plate  modeling  by  Wilson  [1986] 
indicates  deformation  (crustal  shortening)  of  the  southernmost 
section  of  the  Gorda  "plate"  near  the  Mendocino  transform.  The 
northern  Gorda  Ridge  appears  to  have  escaped  this  tectonic 
complexity  [Masson  et  al,  1988],  and  its  spreading  rate  is 
comparable  to  the  southern  Juan  de  Fuca  Ridge. 

Despite  their  similar  spreading  rates  the  Juan  de  Fuca  and 
Gorda  spreading  centers  have  distinctly  different  morphologies, 
with  the  Juan  de  Fuca  Ridge  generally  characterized  by  a  rifted, 
axial  high  while  the  Gorda  Ridge  is  associated  with  a  rift  valley 
(Figure  2).  The  Endeavour.  Northern  Symmetrical  Segment, 
Axial,  Vance,  and  Cleft  Segments  of  the  Juan  de  Fuca  Ridge 
exhibit  a  rifted  axial  high  morphology,  though  the  rift  on  the 
Qeft  Segment  is  subdued.  West  Valley  at  the  northmi  end  of  the 
Juan  de  Fuca  Ridge  has  a  rift  valley.  The  along  axis  bathymetry 
of  most  of  the  central  and  southern  Juan  de  Fuca  Ridge  is 
comparatively  flat  except  for  the  presence  of  Axial  Seamount. 
Only  the  Endeavour  and  West  Valley  Segments  show  significant 
along  axis  variations  in  axial  depth.  Recent  changes  in  the 
location  of  the  spreading  axis  of  the  northern  portion  of  the  Juan 
de  Fuca  Ridge  may  complicate  the  interpretation  of  the 
morphology  of  West  Valley  [Davis  and  Currie,  1993;  K.M.M. 
Rohr  and  K.P.  Furlong,  Birth  of  a  triple  junction,  death  of  a  plate: 
Tectonics  of  the  (Jueen  Charlotte  triple  Junction,  submitted  to 
Geology,  1995].  The  character  of  the  Endeavour  segment  may  be 
affected  by  an  increased  melt  supply  as  a  result  of  overriding  the 
Heckle  melting  anomaly  [Karsten,  1988]  which  has  created  the 
Heckle  seamount  chain  on  the  Pacific  plate  to  the  west  of  the 
ridge.  In  contrast,  all  three  segments  of  the  Gorda  Ridge  have  rift 
valleys  reminiscent  of  the  slow  spreading  Mid-Atlantic  Ridge. 
The  variation  in  axial  depth  along  the  axis  of  the  Gorda  Ridge  is 
also  similar  to  a  slow  spreading  ridge  with  a  deeper  rift  valley 
near  the  segment  ends  than  at  the  segment  mid-point 

There  is  also  a  pronounced  difference  in  seismicity  between 
the  Juan  de  Fuca  Ridge  and  Gorda  Ridge  that  parallels  the 
differences  in  axial  morphology.  While  the  Juan  de  Fuca  Ridge 
is  remarkably  aseismic,  the  Gorda  Ridge  is  very  seismically 
active  with  microeaithquake  focal  depths  extending  to  the  mantle 
(7  ±  2  km  [Johnson  and  Lewis,  1984]  or  15  ±  2  km  depth  [Solano 
and  Bibee,  1985]).  The  earthquake  density  per  55  km  of  ridge 
length  normal  to  the  transform  direction  was  15-20  times  greater 
along  the  Gorda  Ridge  than  along  the  Juan  de  Fuca  Ridge  during 
the  period  1961-1969  [Environmental  Science  Services 
Administration,  1970].  Recent  monitoring  of  the  northeast 
Pacific  using  the  more  sensitive  military  Sounds  Surveillance 
System  (SOSUS)  hydrophone  arrays  has  confinned  the  relative 
seismic  quiescence  of  the  Juan  de  Fuca  Ridge  compared  to  the 
Gorda  Ridge  [Fox  et  at,  1994]. 

The  differences  in  axial  morphology  and  seismicity  are 
consistent  with  the  presence  of  a  stronger  axial  lithosphere 
beneath  the  Gorda  Ridge  than  along  most  of  the  Juan  de  Fuca 
Ridge.  This  variation  of  lithospheric  rheology,  despite  very 
similar  spreading  rates,  may  be  due  to  differences  in  crustal 
thickness  or  thermal  structure  or  some  combination  of  these  two 
effects.  Few  seismic  constraints  exist  on  the  changes  in  crustal 
structure  or  thickness  between  the  Juan  de  Fuca  and  Gorda 
Ridges.  Most  recent  seismic  work  on  the  Juan  de  Fuca  Ridge  has 
focused  on  the  structure  of  the  middle  and  upper  crust  [McClain 
and  Lewis,  1982;  Morton  et  oL,  1987;  Cudrak,  1988;  Rohr  et  at, 
1988;  White  and  Clowes,  1990;  Christeson  et  al.,  1993;  Cudrak 
and  Clowes,  1993;  McDonald  et  al.,  1994].  The  thermal 
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structure  of  the  crust  can  be  estimated  from  the  presence  or 
absence  of  a  magma  body  in  the  crust.  However,  evidence  of 
crustal  magma  bodies  similar  to  those  found  along  the 
morphologically  similar  East  Pacific  Rise  [Detrick  et  a/.,  1987; 
1993]  is  equivocal.  No  seismic  constraints  exist  on  the  presence 
of  crustal  magma  bodies  along  the  Gorda  Ridge  nor  is  seismic 
data  available  to  determine  if  a  systematic  difference  in  crustal 
thickness  exists  between  the  Juan  de  Fuca  and  Gorda  Ridges.  We 
will  use  gravity  data  in  this  paper  to  determine  the  relative 
variations  in  crustal  thickness  and  mantle  thermal  structure 
between  the  Juan  de  Fuca  and  Gorda  Ridge. 

Free  Air  Gravity  and  Bathymetry 

Bathymetry  and  free  air  gravity  anomaly  maps  for  the  Juan  de 
Fuca  and  Gorda  Ridges  are  shown  in  Figure  1  and  Plate  L  The 
bathymetry  map  (Figure  1)  was  compiled  over  the  region  within 
125®-132®W  and  4L5°-49®N.  It  covers  the  entire  Juan  de  Fuca 
Ridge  from  the  Sovanco  Fracture  Zone  down  to  the  Blanco 
Fracture  Zone  and  the  northern  half  of  the  Gorda  Ridge.  The 
trench  on  the  western  margin  of  the  North  American  continent 
forms  the  eastern  boundary  of  the  mapped  region.  To  the  west  of 
the  Juan  de  Fuca  Ridge  the  map  covers  ocean  floor  of  ages  5  My 
and  younger.  The  compiled  data  include  shipboard  bathymetry 
extracted  from  the  National  Geophysical  Data  Center  (NGDC) 
and  multibeam  bathymetry  extracted  from  the  Lamont-Doherty 
Earth  Observatory  RIDGE  (Ridge  InterDiscipIinary  Global 
Studies)  multibeam  data  synthesis  [Keeley  et  at.,  1994].  Between 
ship  tracks  the  bathymetry  was  supplemented  with  DBDB5  data 
(Digital  Bathymetric  Data  Base  5).  Where  multibeam  bathy¬ 
metry  existed  the  NGDC  bathymetry  was  not  used.  All  of  the 
data  were  resampled  at  a  grid  spacing  of  1  km  and  then  the 
surface-fitting-under-tension  algorithm  of  the  Generic  Mapping 
Tools  (GMT)  program  [Smith  and  Wessel,  1990]  was  applied  to 
obtain  the  final  bathymetric  map  (Figure  1 ).  The  accuracy  of  this 
map  varies  vridcly,  being  accurate  to  ~50  m  over  the  ridge  crest, 
which  is  entirely  mapped  with  multibeam  data,  and  degrading 
into  regions  that  were  filled  using  the  DBDB5  data. 

Shipboard  free  air  gravity  data  from  the  NGDC  database  were 
compiled  in  two  subregions  which  are  centered  on  the  Juan  de 
Fuca  and  Gorda  Ridges,  12S'’-132®W,  44®-49®N  and  125°- 
129®W,  41.5®-44®N  respectively  (Plate  I).  To  improve  the 
internal  consistency  of  the  gravity  data,  we  use  the  method  of 
Prince  and  Forsyth  [1984]  to  minimize  discrepancies  in  gravity 
measurements  at  ship  track  crossings.  We  estimate  the  data 
uncertainty  to  be  5  mGal  based  on  the  rms  crossover  error.  After 
applying  appropriate  shifts  to  continuous  segments  of  ship  data, 
we  gridded  the  gravity  data  using  the  surface-fitting-under- 
tension  algorithm  of  the  GMT  program  [Smith  and  Wessel^  1990] 
with  a  grid  spacing  of  1  arc  min  (=~0.0167®).  The  resulting 
gravity  maps  are  shown  in  Plates  la  and  lb.  These  maps  are 
accurate  near  the  ship  tracks  which  are  shown  in  black  dots  since 
no  gravity  data  exist  between  ship  tracks. 

Gravity  Reduction 

The  calculation  of  mantle  Bouguer  anomalies  (MBA)  and 
residual  mantle  Bouguer  anomalies  (RMBA)  follow  the 
procedure  developed  by  Prince  and  Forsyth  [1988]  and  Kuo  and 
Forsyth  [1988].  T^Ie  1  summarizes  the  values  of  the  parameters 
assumed  in  these  calculations.  The  theoretical  gravity  field  was 
calculated  over  the  entire  region  of  the  bathymetry  map.  Because 
mirroring  of  the  grids  introduces  bias  from  the  nonperiodic 


structure  due  to  the  trench  in  calculating  the  Fourier  transforms 
for  the  theoretical  gravity  field,  the  grid  was  padded  using  a 
linear  gradient.  Along  the  ship  tracks  the  theoretical  gravity 
signature  was  subtracted  from  the  free  air  gravity  to  obtain  the 
MBA  and  RMBA.  These  anomalies  were  then  gridded  within  the 
two  subregions  in  the  same  way  as  the  gravity  data  (1  arc  min 
grid  spacing)  to  obtain  maps  of  the  MBA  and  RMBA.  This 
procedure  ensures  that  the  long  wavelength  gravity  signature  of 
the  topography  is  retained  in  the  theoretical  correction  and  that 
the  MBA  and  RMBA  over  the  two  ridge  segments  are 
comparable  having  been  calculated  relative  to  the  same  reference 
levels. 

Mantle  Bouguer  anomaly  maps  for  the  Juan  de  Fuca  and 
Gorda  Ridges  are  shown  in  Plates  Ic  and  Id.  Again,  these  maps 
are  accurate  near  the  ship  tracks  which  are  shown  with  black 
dots.  Calculation  of  the  MBA  removes  most  of  the  short 
wavelength  variations  caused  by  local  topography  leaving  a 
broad,  oval-shaped,  negative  anomaly  aligned  along  the 
spreading  axis  of  the  Juan  de  Fuca  Ridge  (Plate  Ic).  The 
anomaly  is  most  negative  ( — 60  mGal)  over  Axial  Seamount  at 
46®N  and  gradually  becomes  more  positive  south  toward  the 
Blanco  Transform.  The  MBA  on  the  Gorda  Ridge  is  also  most 
negative  along  the  ridge  axis,  becoming  more  positive  toward  the 
Blanco  Fracture  Zone. 

Off  the  axes  of  the  Juan  de  Fuca  and  Gorda  Ridges  the  MBA 
becomes  more  positive,  especially  west  of  the  ridge  crest.  The 
increase  in  MBA  off  axis  is  primarily  due  to  increases  in  density 
of  the  lithosphere  which  thickens  and  cools  as  it  ages.  However, 
part  of  the  gradient  in  MBA  seen  to  the  southwest  of  the  Juan  de 
Fuca  Ridge  reflects  crustal  age  differences  due  to  the  southward 
propagation  of  the  ridge  from  18-1.5  Ma  [Wilson,  1988].  Neither 
the  effects  of  lithospheric  cooling  nor  ridge  propagation  are 
clearly  observed  on  the  eastern  side  of  the  ridge  where  the 
basaltic  crust  becomes  rapidly  buried  with  sediments  that  flood 
the  plate  from  the  continental  margin  [Ludwig  arid  Houtz,  1976]. 
The  low  density  of  the  sediments  (relative  to  the  igneous  crust) 
lead  to  an  artificially  low  MBA  value.  Owing  to  the  effects  of 
both  extensive  sediment  coverage  off  axis  and  ridge  propagation 
on  the  MBA  signal,  we  only  analyze  along  axis  gravity  profiles  in 
this  paper. 

The  RMBA  is  calculated  by  removing  the  effects  of 
lithospheric  cooling  with  age  and  heat  conduction  across  ridge 
offsets  from  the  MBA.  The  along  axis  bathymetry,  free  air 
gravity,  MBA,  and  RMBA  are  shown  in  Figure  3  plotted  versus 
distance  north  of  the  Blanco  Fracture  Zone.  In  an  along  axis 
profile  the  RMBA  corrects  the  MBA  for  cooling  at  ridge  offsets. 
When  the  effect  of  cooling  against  the  Bianco  Fracture  Zone  is 
removed  from  the  MBA,  the  RMBA  signature  along  the  Cleft 
and  Vance  Segments  is  almost  flat.  The  RMBA  beneath  Axial 
Seamount  is  -50  mGal  and  is  localized  within  100  km  of  the 
seamount.  It  increases  to  the  north  and  south  where  the  anomaly 
flattens.  The  RMBA  on  the  Gorda  Ridge  is  elevated  by  10-20 
mGal  relative  to  the  Cleft  and  Vance  Segments  indicating  thinner 
crust  and/or  colder  mantle. 

Modeling  of  Residual  Mantle  Bouguer  Anomalies 

We  interpret  the  gravity  fields  over  the  Juan  de  Fuca  and 
Gorda  Ridges  according  to  two  models  for  the  density  stmeture. 
In  the  first  mode!  the  density  variations  creating  the  observed 
RMBA  are  attributed  only  to  changes  in  the  thickness  of  a 
constant  density  crust  overlying  a  constant  density  mantle.  In  the 
second  model,  lateral  variations  in  both  crustal  thickness  and 


34 


HOOFT  AND  DETRICK:  JUAN  DE  FUCA  AXIAL  MORPHOLOGY  VARIATIONS 


48”N 


46”N 


44“N 


128“W 


132”W 


130”W 


128”W 


44“N 


>90  -80  -70  -60  -50  -40  -30  -20  -10  0  10  20  30  40  50  60  70  80  90 

mGal 


Plate  1.  (a)  and  (b)  Free  air  gravity  and(c)  and  (d)  mantle  Bouguer  anomalies  (MBA)  over  the  Juan  de  Fuca  and 
Gorda  Ridges.  Grid  spacing  is  1  arc  min;  contour  interval  is  10  mGal.  Shipboard  free  air  gravity  data  from  the 
NGDC  data  base  were  compiled  in  two  subregions  of  the  bathymetry  map  centered  over  the  Juan  de  Fuca  and  the 
Gorda  Ridges.  These  free  air  gravity  maps  are  only  accurate  near  the  ship  tracks  which  are  shown  by  black  dots. 
The  free  air  gravity  mimics  the  topography.  The  MBA  is  most  negative  over  Axial  Seamount,  gradually  becoming 
more  positive  toward  the  Blanco  Fracture  Zone.  The  MBA  over  the  Gorda  Ridge  is  also  most  negative  along  the 
ridge  axis,  becoming  more  positive  toward  the  Blanco  Fracture  Zone.  MBA  over  the  Gorda  Ridge  are  more 
positive  than  over  the  southern  Juan  de  Fuca  ridge  indicating  thinner  crust  and/or  colder  mantle. 
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Table  1.  Values  of  Model  Parameters 


Paraineier 

Value 

Gravitadoixal  coostant 

6.6673  X  10'^^  N  m^/kg^ 

Giavitadonal  accclcratioQ 

9.80  ra/s^ 

Wato/cru;^  density  cootrast 

1700  kg/m^ 

Crust/masQtk  dcarity  contra^ 

600kg/m^ 

Mean  crustal  duckness 

6  km 

Half  spreading  rale 

28  mm/yr 

Mande  viscosity 

10*®  Pas 

Mean  asth^ospheric  tci^>erature 

1350®C 

Coefficient  of  diermai  expansion 

3.4  X  lO'^rc 

mantle  density  are  omtroUed  by  temperatures  in  the  upper  mantle 
where  d<xompression  melting  occurs.  Model  1  is  an  example  of 
dynamic  upwelling  [e.g,.  Sparks  and  Parmentier,  1993]  where 
along  axis  variations  in  upwelling  rate  generate  crustal  thickness 


variations;  in  the  extreme  case  there  need  be  no  changes  in 
mantle  temperature.  Model  2  is  representative  of  a  passive  flow 
mode!  where  crustal  thickness  changes  arc  the  result  of  along  axis 
variations  in  mantle  temperature. 

Model  1:  Crustal  Thickiiess  Variations  Only 

RMBA  arc  traditionally  downward  continued  to  obtain  the 
crustal  thickness  variations  required  to  explain  the  observed 
gravity  field  [Kuo  and  Forsyth,  1988],  We  downward  continue 
the  RMBA  to  6  km  below  the  seafloor  afto-  applying  a  low-pass 
filter  with  a  cosine  taper  between  35  and  45  km.  The  crustal 
thickness  variations  arc  thus  referenced  to  a  "normal"  crustal 
thickness  of  6  km.  The  results  are  shown  in  Figure  4a.  The  crust 
beneath  Axial  Seamount  is  inferred  to  be  ~2  km  thicker  than  the 
assumol  mean  thickness  of  6  km.  The  crustal  root  under  Axial 
Seamount  extends  over  100  km  along  axis,  slightly  wider  than  its 
bathymetric  expression  (~80  km).  The  crust/mantle  interface  is 
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Figure  3,  Along-axis  profiles  of  bathymetry,  firee  air  gravity.  MBA,  and  RMBA  plotted  versus  distance  north  of  the 
Blanco  Fracture  Zone,  measured  along  the  ridge  as  on  Figure  I .  The  Blanco  Fracture  Zone  is  indicated  by  the  gray 
bar.  The  Vance  and  Cleft  Segments  have  a  very  flat  along  axis  topography.  In  contrast,  the  northern  Gorda 
Segment  has  a  bowed  along  axis  topography,  shallower  at  the  segment  center  than  at  its  ends.  The  corresponding 
MBA  is  most  negative  (—60  mGai)  over  Axial  Seamount  and  gradually  becomes  more  positive  towards  the  Blanco 
Transform.  The  MBA  on  the  Gorda  Ridge  also  increases  sharply  towards  the  fracture  zone.  When  the  effect  of 
cooling  against  the  Blanco  FrzKtoe  Zone  is  removed  from  the  MBA,  the  RMBA  signature  is  almost  flat  along  the 
Qeft  and  Vance  Segments.  The  RMBA  on  the  Gorda  Ridge  is  elevat«i  by  10-20  mGal  relative  to  the  Qcft  and 
Vance  Segments  indicating  thinner  crust  and/or  colder  mantle. 
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Figure  4.  (a)  Model  I:  Predicted  variations  in  crustal  thickness  assuming  that  the  observed  RMBA  is  due  only  to 
changes  in  the  thickness  of  a  constant  density  crust  overlying  a  constant  density  mantle.  Beneath  Axial  Seamount 
the  crust  is  ~2  km  thicker  than  the  normal  crustal  thickness  of  6  km.  The  Moho  is  very  flat  along  the  Cleft  and 
Vance  Segments  as  is  typical  of  fast  spreading  ridges.  Along  the  northern  Gorda  Ridge,  crust  is  thickest  at  the 
center  of  the  segment  and  thinner  at  the  ends  of  the  segment  (by  -500  m),  as  is  commonly  observed  along  slow 
spreading  ridge  segments.  Crust  on  the  Cleft  Segment  is  thicker  than  that  on  the  northern  Gorda  Segment  by  0.5-1 
km.  (b)  Model  2:  The  predicted  variation  in  crustal  thickness  is  shown  assuming  that  both  mantle  density  and 
crustal  thickness  variations  are  controlled  by  mantle  temperature  [Ito  and  Lin,  1995].  Crustal  thickness  variations 
account  for  70%-75%  of  the  RMBA  signal;  the  remainder  is  due  to  temperature-related  mantle  density  variations. 
(c)  Shows  the  along  axis  mande  temperature  variations  required  to  explain  the  observed  RMBA  (solid  line)  and 
calculated  by  assuming  that  the  axial  topography  is  isostatically  compensated  (dashed  line).  Note  the  small 
differences  in  crustal  thickness  and  mantle  temperature  between  the  Cleft  and  Vance  Segments,  which  have  an  axial 
high,  and  the  northern  Gorda  Ridge,  which  has  a  rift  valley. 


very  flat  along  the  Qeft  and  Vance  Segments  as  is  typical  of  fast 
spreading  ridges.  Along  the  northern  Gorda  Ridge,  crust  is 
thickest  at  the  center  of  the  segment  and  thinner  at  the  ends  of  the 
segment  (by  -500  m),  as  is  commonly  observed  along  slow 
spreading  ridge  segments  (e.g.,  Lin  et  aL,  1990].  Crust  on  the 
Cleft  and  Vance  Segments  is  0.5-1  km  thicker  than  on  the 
northern  Gorda  Segment.  This  difference  in  crustal  thickness  is  a 
maximum  estimate  and  assumes  that  none  of  the  observed 
density  variations  are  due  to  variations  in  crustal  or  mantle 
density.  Thus  at  most,  the  large  difference  in  axial  morphology 
between  the  southern  Juan  de  Fuca  and  Gorda  Ridges  is 
associated  with  a  difference  in  average  crustal  thickness  of  about 
a  kilometer. 

Model  2:  Crustal  Thickness  and  Mantle  Temperature 
Variations 

Crustal  thickness  variations  may-  result  from  mantle  thermal 
anomalies  which  affect  the  amount  of  decompression  melting  and 
also  the  density  structure  of  the  upper  mantle.  The  presence  of  a 


mantle  thermal  anomaly  is  suggested  by  analyses  of  lavas  along 
the  Juan  de  Fuca  Ridge.  The  composition  of  basalts  from  the 
Juan  de  Fuca  Ridge  and  Axial  Seamount  are  transitional  between 
normal  (N-type)  and  enriched  (E-type)  mid-ocean  ridge  basalt 
and  quite  variable  along  axis,  indicating  a  heterogeneous  mantle 
source  [Morgan,  1985].  There  is  no  evidence  in  the  seamount 
lavas  for  geochemical  or  isotopic  enrichment  typical  of  hotspot  or 
mantle  plume  activity.  However,  elevated  values  of  Na20,  CaO 
and  Sr,  at  a  given  MgO  content,  differentiate  the  seamount  lavas 
from  those  of  the  adjoining  ridge  segments  [Rhodes  et  al,  1990], 
suggesting  that  it  is  associated  with  a  zone  of  anomalously  warm 
mantle.  These  observations  suggest  that  both  crustal  thickness 
variations  and  temperature-related  mantle  density  variations  may 
be  responsible  for  the  observed  gravity  anomalies. 

Ito  and  Lin  (1995]  developed  a  model  in  which  crustal 
thickness  and  mantle  density  variations  are  both  controlled  by 
mantle  temperature:  crustal  thickness  by  temperature  enhanced 
melting  and  mantle  density  by  thermal  expansion.  The  linear 
melting  function  in  their  model,  which  incorporates  passive 
mantle  flow,  heat  transport,  and  decompression  melting,  gives 
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rise  to  a  linear  relationship  between  mantle  temperature,  AT  at  a 
depth  of  100  km,  variations  in  MBA,  AMBA,  and  bathymetry, 
AH: 


AT(MBA)  =  .0.576  AMBA 

a  misfit  <  3®C 

(1) 

AT(H)=48.3  AH 

<j  misfit  <  3®C 

(2) 

Using  these  predicted  relationships  we  estimated  the  along 
axis  variations  in  crustal  thickness  and  mantle  temperature 
required  to  explain  the  RMBA  observol  along  the  Juan  de  Fuca 
and  Gorda  Ridges  (Figures  3c  and  3d).  We  use  the  RMBA 
instead  of  the  MBA  to  account  for  the  cooling  effect  at  segment 
ends.  Calculated  crustal  thickness  variations  (Rgure  4b)  now 
give  rise  to  70%-75%  of  the  RMBA  signal.  The  remaining  25%- 
30%  of  the  topography  and  gravity  signal  is  supplied  by  the 
anomalously  hot  and  less  dense  mantle  beneath  the  rise  axis.  The 
solid  line  in  Figure  4c  shows  the  along  axis  mantle  temperature 
vari^ons  (AT(RMBA))  required  to  explain  the  observed  RMBA 
calculated  from  (1).  An  independent  estimate  of  the  along  axis 
variations  in  mantle  temperature  (AT(H)),  shown  by  the  dashed 
line  in  Figure  4c,  was  made  from  (2)  which  includes  the 
assumption  that  the  axial  topography  is  isostatically 
compensated.  For  this  calculation  axial  depth  was  averaged  over 
a  30-km-wide  across-axis  zone  to  minimize  errors  that  might  be 
introduced  if  dynamic  stresses  are  contributing  to  axial  relief. 
The  two  estimates  show  the  same  long- wavelength  characteristics 
with  the  highest  mantle  temperatures  beneath  Axial  Seamount, 
progressively  cooler  temperatures  north  and  south  along  the  Juan 
de  Fuca  Ridge  away  from  Axial  Seamount,  and  the  lowest  mantle 
temperatures  beneath  the  Gorda  Ridge.  The  bathymetrically 
derived  AT  estimate  displays  more  along  axis  variation  than  the 
RMBA-derived  estimate.  This  is  because  the  bathymetry 
contains  more  short  wavelength  energy  than  the  gravity  field. 
The  difference  between  the  two  mantle  temperature  estimates 
(e.g.,  along  the  Cleft  segment)  could  indicate  some  dynamic 
support  for  the  along  axis  depth  variations. 

The  results  presented  in  Figure  4  provide  the  best  available 
constraints  on  the  differences  in  crustal  thickness  and  mantle 
temperature  along  the  Gorda/Juan  de  Fuca  ridge  system.  Mantle 
temperature  beneath  Axial  Seamount  may  be  elevated  by  as 
much  as  3(r-40®C,  and  the  crust  may  be  up  to  30^  thicker  than 
that  formed  elsewhere  along  the  Juan  de  Fuca  Ridge.  The  mantle 
thermal  and  crustal  thickness  anomaly  associated  with  Axial 
Seamount  decays  rapidly  along  the  Juan  de  Fuca  Ridge,  and  the 
Qeft  and  Vance  Segments  are  associated  with  crust  that  is  only 
300-700  m  thicker  and  with  mantle  temperatures  that  are  only 
10°-15°C  higher  than  that  estimated  beneath  the  northern  Gorda 
Ridge.  Model  1  represents  the  maximum  possible  crustal 
thickness  variations,  yet  the  difference  in  crustal  thickness 
between  the  northern  Gorda  and  southern  Juan  de  Fuca  Ridges  is 
only  0.5-1  km.  Thus  our  results  suggest  that  comparatively  small 
differences  in  crustal  thickness  and/or  mantle  temperature  can 
produce  dramatic  changes  in  axial  morphology  at  these 
intermediate  spreading  rates.  In  the  following  section  we 
compare  these  observations  with  model  predictions  of  how 
changes  in  crustal  thickness  affect  the  rift  vallcy/no  rift  valley 
transition.  We  also  consider  the  origin  of  the  vari^on  in  crustal 
thickness  and  mantle  temperature  inferred  along  the  Gorda/Juan 
de  Fuca  ridge  system. 

Discussion 

A  number  of  different  mechanisms  have  bo^n  proposed  to 
explain  why  some  mid-ocean  ridges  have  rift  valleys  and  others 


do  not  (e.g.,  Sleep,  1969;  Lauchenbruch,  1973;  Tapponnier  and 
Francheieau,  1978;  Sleep  and  Rosendahl,  1979;  Phipps  Morgan 
et  al.,  1987;  Chen  and  Morgan,  1990].  Over  the  past  few  years  a 
growing  consensus  has  developed  around  a  hypothesis  originally 
proposed  by  Tapponnier  and  Francheteau  (1978]  that  median 
valleys  form  as  a  result  of  stretching  and  thinning  of  the 
lithosphere  at  the  ridge  axis.  Chen  and  Morgan  (1990]  showed 
that  a  combination  of  brittle-plastic  and  viscous  rheologies 
undergoing  extension  at  a  slow  spreading  ridge  can  produce  an 
axial  valley  of  the  right  dimensions  provided  that  the  crust  is  cool 
enough  to  behave  in  a  brittle  fashion.  At  higher  spreading  rates, 
like  those  found  along  the  East  Pacific  Rise,  they  propose  that  a 
thin,  brittle  surface  layer  is  decoupled  from  the  stronger  mantle 
below  by  the  presence  of  a  crustal  magma  body  and  a  hot,  ductile 
lower  crust.  The  reduced  strength  of  the  crust  allows  local 
isostatic  forces  to  create  an  axial  high.  This  model  has  recently 
been  refined  by  Phipps  Morgan  and  Chen  (1993a,  b]  and 
Neumann  and  Forsyth  (1993]  and  extended  to  explore  the  effect 
of  variable  crustal  thickness  on  the  thermal  and  rheological 
properties  of  zero-age  lithosphere  and  axial  morphology. 
According  to  these  models  the  transition  in  axial  morphology 
from  rift  valley  to  axial  high  depends  critically  on  the  thermal 
structure  of  the  rise  axis,  which  vrill  be  controlled  by  the  balance 
between  heat  supplied  by  magmatic  intrusion  and  solidification 
and  heat  removed  by  hydrothermal  circulation. 

Our  observations  on  the  Juan  de  Fuca/Gorda  system,  which  is 
spreading  at  an  intermediate  half  rate  of  28  mm/yr,  support  the 
sharp  transition  in  ridge  morphology  predicted  by  the  Phipps 
Morgan  and  Chen  (1993b]  model.  The  gravity  data  show  no 
large  difference  in  density  structure  between  the  two  ridges. 
Variations  in  crustal  thickness  are  at  most  ~1  km.  If  a  mantle 
thermal  anomaly  is  postulated  as  the  source  of  these  crustal 
thickness  variations  and  the  density  anomalies  are  partitioned  in  a 
reasonable  way  between  mantle  thermal  density  variations  and 
crustal  thickness  variations  [I to  and  Lin,  1995],  then  the 
differences  in  average  crustal  thickness  between  the  southern 
Juan  de  Fuca  and  the  northern  Gorda  Ridge  may  be  as  small  as 
300-700  m.  Thus  our  results  support  predictions  that  the 
transition  from  rift  valley  to  axial  high  will  occur  abruptly  with  a 
small  change  in  crustal  thickness  and/or  mantle  temperature  at 
intermediate  spreading  rates. 

The  crustal  thickness  and  mantle  temperature  variations  we 
have  inferred  along  the  Gorda/Juan  de  Fuca  Ridge  system  may  be 
related  to  the  presence  of  the  Cobb  thermal  anomaly  beneath  the 
central  Juan  de  Fuca  Ridge  (Figure  5),  The  Cobb  hotspot  has 
been  proposed  as  the  source  of  the  Cobb-Eickelberg  seamount 
chain  which  intersects  the  Juan  de  Fuca  Ridge  at  46®N  [Duncan 
and  Clague,  1985;  Smoot,  1985].  Axial  Seamount  lies  at  the 
intersection  of  the  Cobb  Eickelberg  seamount  chain  and  the  Juan 
de  Fuca  Ridge,  and  it  is  believed  to  mark  the  present  location  of 
the  Cobb  hotspot  The  thermal  anomaly  associated  with  the 
Cobb  hotspot  which  we  estimated  in  this  study  to  be  about  30®- 
40®C,  appears  to  influence  a  large  portion  of  the  Juan  de  Fuca 
Ridge,  perhaps  through  along  axis  channeling  of  the  hotter, 
upwclling  mantle  along  the  ridge  axis  as  originally  proposed  by 
Vogt  and  Johnson  (1975].  The  Gorda  Ridge  may  be  isolated 
from  the  influence  of  the  Cobb  thermal  anomaly  by  the  400-km- 
long  Blanco  Transform  which  blocks  the  southward  flow  of  this 
material. 

Tlte  hotter  upwclling  along  the  Juan  de  Fuca  Ridge  results  in 
the  formation  of  thicker  crust  than  along  the  Gorda  Ridge 
because  of  a  greater  extent  of  decompression  ntelting.  This  effect 
is  most  pronounced  beneath  Axial  Seamount  where  the  crust  is 
~309&  thicker  than  normal.  Elsewhere  along  the  Juan  de  Fuca 
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Figure  5.  CartCK)n  illustrating  the  effect  of  the  Cobb  thermal 
anomaly  located  under  Axial  Seamount  on  the  Juan  de  Fuca 
Ridge.  Hotter  upwelling  mantle  results  in  warmer,  less  dense 
mantle,  a  thinner  lithosphere,  more  decompression  melting,  and  a 
thicker  crust  beneath  the  Juan  de  Fuca  Ridge,  The  weak 
lithosphere  forms  an  axial  high.  The  Blanco  Fracture  Zone 
isolates  the  Gorda  Ridge  from  the  influence  of  the  Cobb  thermal 
anomaly.  On  the  northern  Gorda  Ridge,  crust  is  thinner,  the 
mantle  is  colder,  and  the  lithosphere  is  thicker.  The  colder 
thermal  regime  results  in  a  stronger  lithosphere  which  deforms  to 
produce  a  rift  valley. 


Ridge  (except  perhaps  at  West  Valley),  the  crust  may  be  up  to 
700  m  thicker  than  that  forming  along  the  northern  Gorda  Ridge 
where  the  mantle  is  colder.  The  thicker  crust  and  hotter  mantle 
result  in  the  formation  of  thinner,  weaker  lithosphere  along  the 
Juan  de  Fuca  Ridge  than  at  the  Gorda  Ridge.  The  stronger 
lithosphere  found  along  the  Gorda  Ridge  can  deform  by  brittle- 
plastic  necking  leading  to  the  formation  of  a  rift  valley,  while  the 
weaker,  thinner  lithosphere  along  the  Juan  de  Fuca  Ridge  is 
locally  compensated  forming  an  axial  high. 

Conclusions 

The  southern  Juan  de  Fuca  Ridge  has  an  axial  high 
morphology,  typical  of  faster  spreading  ridges,  while  the  northern 
Gorda  Ridge  has  a  rift  valley  morphology,  as  is  found  at  slower 
spreading  ridges,  even  though  they  are  opening  at  the  same 
spreading  rate  (28  mm/yr  half  rate).  The  Gorda  Ridge  is  also 
characterized  by  much  higher  levels  of  seismicity  suggesting  the 
presence  of  a  thicker,  stronger  lithosphere.  Residual  mantle 
Bouguer  anomalies  along  the  northern  Gorda  Ridge  are  elevated 
10-20  mGal  relative  to  the  southern  Juan  de  Fuca  Ridge 
indicating  the  presence  of  a  thinner  crust  and/or  colder  mantle. 
We  attribute  these  systematic  differences  in  axial  morphology, 
gravity,  and  seismicity  to  along  axis  variations  in  crustal 
thickness,  mantle  temperature,  and  lithospheric  strength  due  to 
the  presence  of  the  Cobb  thermal  anomaly  at  Axial  Seamount 
along  the  central  Juan  de  Fuca  Ridge,  We  estimate  that  mantle 
temperatures  beneath  Axial  Seamount  are  elevated  and 

the  crust  is  -30%  thicker  than  that  forming  elsewhere  along  the 
Juan  de  Fuca  Ridge.  The  mantle  thermal  and  crustal  thickness 
anomalies  associated  with  the  Cobb  thermal  anomaly  decay 
rapidly  along  axis  and  the  Vance  and  Cleft  Segments  at  the 
southern  end  of  the  Juan  de  Fuca  Ridge  are  associated  with  crust 
that  is  only  300-700  m  thicker  and  mantle  temperatures  that  are 
only  10®-15°C  higher  than  that  estimated  beneath  the  Gorda 
Ridge.  Even  when  all  of  the  along  axis  gravity  anomaly  is 


attributed  to  crustal  thickness  variations,  the  maximum  difference 
in  crustal  thickness  between  the  Gorda  and  southern  Juan  de  Fuca 
Ridge  is  only  -1  km.  The  thicker  crust  and  hotter  mantle  along 
the  Juan  de  Fuca  Ridge  result  in  the  formation  of  a  thinner, 
weaker  axial  lithosphere  than  at  the  Gorda  Ridge.  Although  the 
differences  in  crustal  thickness  (-0.5  km)  and/or  mantle 
temperature  (10®-15°C)  that  we  infer  between  the  Juan  de  Fuca 
and  Gorda  Ridges  are  not  large,  they  are  sufficient  to  produce  a 
major  change  in  lithospheric  strength  and  axial  morphology  at 
intermediate  spreading  rates. 
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Chapter  3 

Seismic  Structure  and  Indicators  of  Magma  Budget 
ALONG  the  Southern  East  Pacific  Rise 
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Abstract 


We  re-examine  the  relationship  between  seismically-constrained  variations  in 
crustal  structure  and  indicators  of  ridge  magma  supply  along  the  southern  East  Pacific 
Rise.  Over  the  past  decade,  work  on  the  fast  spreading  East  Pacific  Rise  has  shown  the 
existence  of  significant,  segment-scale  variations  in  depth,  morphology,  crustal  structure, 
basalt  geochemistiy,  and  hydrothermal  activity  that  have  often  been  attributed  to  along-axis 
differences  in  magma  supply  to  the  ridge.  Broad,  shallow  ridge  segments  have  been 
assumed  to  reflect  a  high  magma  supply  and  are  predicted  to  have  larger,  more  robust 
magma  chambers,  erupt  more  primitive  basalts,  form  thicker  crust,  and  have  more 
vigorous  hydrothermal  activity  than  sections  of  the  East  Pacific  Rise  with  a  deeper, 
narrower  rise  crest.  Here  we  use  multichannel  seismic  (MCS)  reflection  data  from  an  800- 
km-long  section  of  the  southern  East  Pacific  Rise  (SEPR)  to  quantitatively  examine  the 
correlation  between  along-axis  variations  in  crustal  structure  and  other  indicators  of  magma 
supply  on  the  regional  scale.  Good  correlations  exist  along  the  SEPR  between  axial  depth, 
ridge  cross-sectional  area,  mantle  Bouguer  anomaly,  and  the  MgO  wt%  of  basalts 
recovered  from  the  rise  axis.  These  correlations  suggest  along-axis  changes  in  cmstal 
thickness  and  temperature  consistent  with  variations  in  magma  supply  on  time  scales  of 
-100  Kyr.  However,  we  show  that  the  magma  sill  properties  (depth  and  width  of  the  melt 
lens),  seismic  layer  2A  thickness,  and  the  intensity  of  hydrothermal  venting  are  poorly 
correlated  with  regional  variations  in  ridge  depth  and  cross-sectional  area.  We  suggest  that 
the  emplacement  geometry  (width  of  the  intrusion  zone  and  flow  lengths),  not  magm?i 
supply,  controls  extrusive  layer  (seismic  2A)  thickness.  The  depth  of  the  magma  sill  is 
very  sensitive  to  the  local  thermal  structure.  Magma  lens  properties  and  hydrothermal 
activity  are  closely  linked  to  spreading  events  (dike  intrusion,  eruptions,  faulting)  which 
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occur  on  much  shorter  time  scales  (ca.  10-100  yr)  than  the  longer-term  variations  in 
magma  supply  reflected  in  along-axis  changes  in  ridge  morphology. 


Introduction 

"As  axial  width  seems  to  reflect  the  magmatic  budget,  one  might  expect  a 
correlation  between  axial  width  and  the  presence  of  hydrothermal  activity. 

It  is  also  intriguing  to  speculate  that  axial  shape  may  reflect  the  width  or 
depth  of  axial  magma  chambers" .  [Langmuir  et  al.,  1986] 

Over  the  past  decade,  work  on  the  fast  spreading  East  Pacific  Rise  (EPR)  has 
shown  the  existence  of  significant,  segment-scale  variations  in  axial  depth  [Macdonald  et 
al,  1984;  Macdonald  and  Fox,  1988],  axial  volume  [Scheirer  and  Macdonald,  1993],  axial 
gravity  anomaly  [Madsen  et  al,  1990;  Magde  et  al,  1995],  magma  chamber  depth  and 
width  [Detrick  et  al,  1987;  Detrick  et  al,  1993],  basalt  geochemistry  [Langmuir  et  al, 
1986;  Sinton  et  al,  1991;  Mahoney  et  al,  1994],  and  hydrothermal  activity  [Haymon  et 
al,  1991;  Baker,  1996].  This  along-axis  variability  is  largely  independent  of  spreading 
rate  and  has  often  been  attributed  to  differences  in  magma  supply  to  the  ridge  within  or 
between  segments  [Macdonald  et  al,  1984;  Macdonald  and  Fox,  1988;  Macdonald  et  al, 
1992].  Broad,  shallow  sections  of  ridge  are  assumed  to  reflect  a  high  magma  supply  and 
are  predicted  to  have  larger,  more  robust  axial  magma  chambers,  erupt  more  primitive 
basalts,  form  thicker  crust,  have  an  axial  summit  caldera  (ASC),  and  have  more  vigorous 
hydrothermal  circulation  than  sections  on  the  EPR  with  a  deeper,  narrower  rise  crest 
(Figure  1). 
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The  strongest  evidence  for  this  simple  magma  supply  model  has  come  from  the 
systematic,  segment-scale  variations  in  axial  depth  and  morphology  observed  along  the 
EPR  [Macdonald  et  ai,  1984;  Macdonald  and  Fox,  1988;  Lonsdale,  1989;  Scheirer  and 
Macdonald,  1993],  and  the  apparent  correlation  of  these  variations  with  the  presence  or 
absence  of  a  reflection  from  what  has  been  interpreted  as  a  magma  sill  at  the  roof  of  the 
magma  chamber  [Detnck  et  ai,  1987;  Harding  et  ai,  1989;  Kent  et  al,  1990;  Detrick  et 
ai,  1993],  The  axial  high  at  the  EPR  is  typically  shallowest  near  the  middle  of  segments 
and  deepens  significantly  toward  segment  offsets  [Macdonald  et  al,  1984;  Lonsdale, 
1989].  These  variations  in  axial  depth  are  associated  with  systematic  changes  in  the  shape 
of  the  axial  high  [Macdonald  and  Fox,  1988;  Scheirer  and  Macdonald,  1993].  Where  the 
ridge  is  shallowest  the  axial  high  generally  has  a  broad,  flat-topped  cross-section;  where 
the  ridge  axis  is  deeper  the  axial  high  usually  has  a  narrow,  triangular  shape.  Based  on  this 
morphological  evidence,  Macdonald  etal  [1984]  and  Macdonald  and  Fox  [1988]  suggest 
that  the  magma  supply  to  the  crust  is  greater  near  segment  centers  resulting  in  a  larger, 
more  robust  magma  chamber  and  a  shallower,  more  inflated  ridge  than  near  segment  ends. 

This  interpretation  is  supported  by  the  correlation,  first  noted  by  Macdonald  and 
Fox  [1988]  for  the  9-13°N  EPR  area,  between  axial  morphology  and  the  presence  or 
absence  of  the  AMC  reflector  mapped  by  Detnck  et  al  [1987].  Where  the  magma  sill  is 
present,  such  as  along  the  ridge  segment  immediately  south  of  the  Clipperton  transform, 
the  axial  high  is  broad  and  flat-topped,  and  a  narrow  (<200m),  shallow  (<50  m)  axial 
summit  caldera  (ASC)  is  recognized  along-axis  (except  in  the  case  where  the  ASC  has  been 
filled  by  recent  lava  flows).  In  contrast,  where  the  rise  axis  is  deep,  narrow  and  more 
triangular  in  cross-section,  such  as  along  the  ridge  segment  immediately  north  of  the 
Clipperton  transform,  a  magma  lens  reflector  is  not  detected  and  an  ASC  is  absent. 
Seismic  reflection  data  from  the  southern  EPR  are  generally  consistent  with  this  association 
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between  shallow,  inflated  sections  of  the  ridge  and  the  presence  of  a  mid-crustal  magma  siU 
[Scheirer  and  Macdonald,  1993].  For  example.  Mutter  et  al.  [1995]  report  that  between 
17°S  and  17°30'S,  the  melt  lens  becomes  shallower  and  wider  as  the  ridge  axis  becomes 
shallower  and  broader,  and  use  deviations  from  this  expected  relationship  to  suggest  very 
recent  or  ongoing  magmatic  activity  near  17°26'S.  Similarly,  it  has  been  suggested  that  the 
thickness  of  the  extrusive  layer  reflects  the  magma  budget  {Harding  et  al.,  1993].  The 
evidence  for  recent  volcanism  [Haymon  et  al,  1993]  and  abundant  hydrothermal  activity 
{Haymon  et  al.,  1991]  along  the  northern  EPR  near  9  SON,  and  the  general  good 
correlation  between  morphological  indicators  of  magma  supply  and  the  presence  of 
hydrothermal  plumes  [Baker,  1996],  are  also  consistent  with  the  magma  supply  model 
summarized  in  Figure  1. 

However,  new  results  suggest  that  the  relationship  between  morphological 
indicators  of  magma  supply  and  axial  crustal  stracture  is  more  complicated  than  this  simple 
model  would  predict.  Hussenoeder  et  al.  [1996]  have  shown  that  the  thickness  and 
inferred  crystal  content  of  the  magma  lens  along  the  northern  and  southern  EPR  are 
uncorrelated  with  axial  depth.  Variations  in  magma  sill  depth  have  been  reported  from  both 
the  northern  [Harding  et  al,  1993]  and  southern  [Detrick  et  al,  1993]  EPR  that  are  not 
related  to  changes  in  the  depth  or  cross-sectional  area  of  the  ridge.  Along  the  northern 
EPR,  the  widest  magma  sill  is  not  found  at  9°50'  N  where  the  ridge  is  shallowest  and 
broadest,  but  where  the  ridge  is  deepest  just  north  of  the  9°03'N  overlapping  spreading 
center  (OSC)  [Kent  et  al,  1993a].  Barth  and  Mutter  [1996]  report  that  the  crust  is  actually 
thinner  beneath  the  inflated  portion  of  the  ridge  at  9°50'N  than  elsewhere  along  this  ridge 
segment.  Similarly,  Begnaud  and  McClain  [1996]  and  van  Avendonk  et  al  [1997]  find 
little  difference  in  the  thickness  of  crust  produced  on  either  side  of  the  Clipperton 
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transform,  despite  large  differences  in  the  depth  and  morphology  of  the  EPR  on  either  side 
of  the  fracture  zone. 

We  examine  the  relationship  between  various  indicators  of  magma  supply  and 
seismically-constrained  variations  in  crustal  structure  along  more  than  800  km  of  the 
southern  EPR.  The  variables  we  compare  are:  axial  depth,  cross-sectional  area  of  the 
ridge,  layer  2A  (extrusive  layer)  thickness  on-  and  off-axis,  the  rate  at  which  the  extmsive 
layer  thickens  off-axis,  magma  sill  depth  and  width,  mantle  Bouguer  anomaly,  temperature 
of  the  erupted  lavas  (MgO  content),  hydrothermal  activity,  and  the  presence  of  an  axial 
summit  caldera.  We  will  show  that  the  thickness  of  the  extrusive  section  (seismic  layer 
2A),  the  properties  of  the  magma  sill  (depth,  width,  thickness),  and  the  intensity  of 
hydrothermal  venting  are  poorly  correlated  with  morphological  indicators  of  magma  supply 
such  as  axial  depth  or  ridge  cross-sectional  area.  We  explain  these  observations  by 
proposing  that  the  emplacement  geometry  (width  of  the  dike  intrusion  zone  and  flow 
lengths),  not  magma  supply,  controls  the  extmsive  (layer  2A)  thickness.  We  also  suggest 
that  the  melt  lens  properties  and  hydrothermal  activity  are  closely  linked  to  spreading  events 
(dike  intrasion,  emptions  and  faulting)  which  occur  on  much  shorter  time  scales  (ca.  10- 
100  yr)  than  the  variations  in  magma  supply  and  thermal  state  which  are  reflected  in  the 
ridge  morphology. 

Data  Processing  and  Analysis 

In  1991  a  multichannel  seismic  reflection  survey  imaged  an  -800  km  long  section 
of  the  southern  East  Pacific  Rise  (Figure  2)  from  the  northern  limb  of  the  20.7°S 
overlapping  spreading  center  (OSC)  to  the  Garrett  fracture  zone  [Detrick  etal,  1993].  This 
extensive  two  ship  reflection/refraction  survey  concentrated  on  three  areas  centered  at 
14°15'S,  17°20'S,  and  over  the  15°55'S  OSC  (Figure  2).  Outside  of  these  regions. 
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reflection  data  were  shot  along  the  ridge  axis  and  five  additional  across-axis  seismic  lines 
were  obtained.  The  reflection  data  were  acquired  aboard  the  research  vessel  Maurice 
Ewing  of  the  Lamont-Doherty  Earth  Observatory  using  a  4-km-long,  160-channel  digital 
streamer  with  25  m  group  separation.  The  source  was  a  20-gun,  8385-cubic-inch  air  gun 
array.  For  the  lines  north  of  17°20'S,  shots  were  fired  every  20  s,  yielding  a  shot  spacing 
of  50  m  and  ^fold  coverage  for  a  nominal  ship  speed  of  5  knots.  South  of  this  area,  port 
and  starboard  halves  of  the  array  were  alternately  fired  every  10  s,  yielding  a  shot  spacing 
of  25  m  and  80-fold  coverage. 

The  multichannel  seismic  data  were  processed  using  SIOSEIS  (written  by  Paul 
Henkart,  Scripps  Institution  of  Oceanography,  La  Jolla,  Ca.).  Shot  gathers  were  sorted 
into  12.5-m  bins  producing  common  mid-point  (CMP)  gathers  whose  ranges  extended 
between  263  and  4263  m.  The  along-axis  data  were  stacked  using  a  velocity  model  that 
varied  with  distance  along  the  ridge  crest.  The  normal  moveout  (NMO)  velocities  of  the 
seafloor,  layer  2A,  and  magma  sill  reflections  were  determined  from  constant  velocity 
stacks  of  the  data  constructed  at  locations  where  preliminary  stacks  showed  a  change  in 
crustal  structure.  Seafloor  and  postcritical  reflections  near  the  base  of  layer  2A  were  muted 
at  ranges  greater  than  3200  m  to  prevent  excessive  stretching  that  occurs  within  the  steep 
velocity  gradient  at  the  base  of  layer  2A.  Muting  of  the  farthest  offsets  (>3200  m)  of  the 
layer  2A  reflection  essentially  includes  only  those  rays  which  turn  near  the  base  of  layer  2A 
thereby  minimizing  any  errors  associated  with  stacking  these  wide-angle  reflections 
[Harding  et  al,  1993].  The  across-axis  reflection  profiles  were  stacked  using  a  velocity 
model  determined  from  constant  velocity  stacks  made  every  2.5  km.  Exact-log  dip 
moveout  (DMO)  corrections  [Liner,  1990]  were  applied  prior  to  stacking.  Finally,  the 
profiles  were  finite  difference  time  migrated  using  a  45°  algorithm  [Brysk,  1983]  to  account 
for  the  2-D  velocity  variations  produced  by  seafloor  and  layer  2A  topography.  On  the 
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resulting  profiles,  the  magma  sill  widths  may  be  slightly  over-estimated  due  to  energy 
diffracting  off  the  edges  of  the  sill  [Kent  et  ai,  1990],  however,  these  errors  are  minimized 
by  the  finite  difference  migration. 

The  two-way  travel  time  to  the  layer  2A  and  magma  sill  reflectors  were  picked  from 
the  along-  and  across-axis  seismic  profiles  at  100  m  spacing.  The  on-axis  data  picks  are 
not  significantly  affected  by  off-axis  wander  of  the  ship  or  streamer  during  the  experiment. 
There  is  a  good  agreement  between  picks  of  on-axis  layer  2A  thickness  and  magma  sill 
depth  made  from  the  across-axis  profiles  with  those  made  at  the  intersection  points  on  the 
along-axis  profile.  The  excellent  navigation  of  this  experiment  is  confirmed  by  plots  of  the 
ship  track  superimposed  on  the  multibeam  bathymetry  of  the  rise  crest.  Even  in  regions 
where  the  ship  deviates  from  the  center  of  the  axial  high  there  is  no  consistent  pattern  of 
disappearance  of  the  magma  sill,  or  thickening  of  layer  2A.  Two-way  travel  time  picks 
were  converted  to  depth  using  the  composite  velocity  model  of  Harding  etal.  [1993]  which 
includes  a  low  velocity,  low  gradient  layer  (2.2  to  2.3  km/s)  beneath  the  seafloor  underlain 
by  a  layer  with  a  high  velocity  gradient  (2.3  to  5.6  km/s).  The  ratio  of  the  thicknesses  of 
these  two  layers  is  held  fixed  (95:160  m)  while  their  total  thickness  is  varied  to  match  the 
observed  travel  time.  The  travel  time  between  the  base  of  layer  2A  and  the  magma  sill 
reflector  was  converted  to  thickness  by  assuming  a  linear  gradient  in  velocity  of  0.74  km/s 
per  km  beneath  the  base  of  layer  2  A. 

On  the  across-axis  multichannel  seismic  lines  we  measured  the  depth  of  the  magma 
sill  below  the  seafloor  and  its  width,  the  on-  and  off-axis  thickness  of  layer  2A  and  the  rate 
at  which  the  thickening  occurs.  The  water  depth,  cross-sectional  area  [Scheirer  and 
Macdonald,  1993],  latitude,  and  measures  of  seismic  structure  for  each  of  the  across-axis 
lines  are  given  in  Table  1.  In  this  table  we  have  also  compiled  the  information  from  the 
across-axis  seismic  profiles  shot  in  the  intensive  study  area  at  14°15'S  [Kent  et  ai,  1994], 
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and  at  9°N  on  the  northern  EPR  [Harding  et  al,  1993;  Kent  et  ai,  1993a&b].  All  the 
seismic  lines  from  the  northern  and  southern  EPR  were  processed  and  interpreted  using 
exactly  the  same  methods  so  inter-comparison  of  the  observations  is  particularly 
meaningful.  Correlation  coefficients  between  these  variables  for  the  along-  and  across- 
axis  lines  are  given  in  Table  2. 

Results 

Variation  in  Seismic  Structure  along  the  SEPR 

The  SEPR  has  a  relatively  uniform  depth  from  the  Garrett  transform  to  about  18°S, 
and  then  deepens  systematically  towards  the  large  overlapping  spreading  center  (OSC)  at 
20.7°S  (Figure  3a).  The  regional  variation  in  crustal  structure  along  this  portion  of  the 
SEPR  can  be  divided  into  northern,  central,  and  southern  ridge  sections  bounded  by  the 
Garrett  transform,  the  15°55'S  and  17°55'S  OSCs,  and  the  large  20.7°S  OSC,  respectively. 
These  sections  of  ridge  are  defined  based  on  differences  in  the  range  and  variability  of  the 
observed  crustal  structure  as  determined  from  the  regional  variability  in  seafloor  depth, 
extrusive  layer  thickness  and  magma  sill  depth.  The  magmatic  segments  defined  by  Sinton 
et  al.  [1991]  (Figure  3c)  are  contained  within  these  three  ridge  sections;  southern  section  = 
magmatic  segments  G,H,  I,  and  J;  central  section  =  magmatic  segments  L  and  K;  northern 
section  =  magmatic  segments  M,  N,  and  O.  Along  the  northern  section  of  ridge  the  axis  is 
broad,  shallow  and  remarkably  uniform  in  depth.  The  axis  remains  shallow  and  is  very 
broad  along  the  central  section  but  shows  a  greater  along-axis  variability  in  seafloor  depth. 
Between  17°S  and  18°S  the  ridge  has  the  most  inflated  morphology  of  the  entire  SEPR 
region.  Along  the  southern  section  the  ridge  deepens  and  narrows  rapidly  from  a  shallow 
point  in  the  "Hump"  region  at  18°27'S  to  the  large  20.7°S  OSC.  The  along-axis  variability 
in  axial  depth  increases  from  north  to  south.  The  intensive  seismic  studies  were  located  at 
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the  15°55'S  OSC,  and  at  latitudes  14°15'S  and  17°20'S  where  the  ridge  is  broad  and 
shallow  and  has  the  morphological  characteristics  of  a  high  magma  budget . 

Figure  4  shows  four  examples  of  across-  and  along-axis  multichannel  profiles  from 
the  SEPR.  These  seismic  profiles  illustrate  the  crustal  structure  of  parts  of  the  ridge 
ranging  from  inflated  to  magma-starved  in  morphology.  Figure  4a  is  from  the  17°20'S 
intensive  study  area,  a  portion  of  the  ridge  along  the  central  segment  with  an  extremely 
inflated  morphology  (cross-sectional  area  =  5.6  km^).  The  top  panel  shows  an  across-axis 
profile  at  17°26'S  (SEPR27)  and  illustrates  the  general  features  of  the  crustal  structure  of  a 
fast-spreading  ridge.  The  seafloor  is  shallowest  at  the  axis  and  deepens  on  either  side  of 
the  ridge.  In  cross-section  the  ridge  has  a  broad,  domed  shape.  Wide-angle  reflections 
from  the  base  of  layer  2A  are  imaged  ~0. 15  s  below  the  seafloor  at  the  ridge  axis.  Seismic 
layer  2A,  which  we  correlate  with  the  volcanic  extmsive  layer  [Herron,  1982;  Harding  et 
ai,  1989;  Christeson  et  al,  1992],  is  230  m  thick  at  the  ridge  axis  and  doubles  in  thickness 
(570  m)  within  3  km  on  either  side  of  the  ridge.  The  magma  sill  lies  at  3.85  s,  a  depth  of 
760  m  below  the  seafloor,  and  has  a  width  of  850  m. 

As  the  ridge  deepens  and  narrows  the  ridge  morphology  and  crustal  structure 
change.  Figure  4b  and  c  are  from  the  northern  and  southern  segments  (latitudes  14°  15'  S  - 
SEPR70  and  18°28'S  -  SEPR06)  respectively,  regions  with  a  morphology  indicative  of  a 
medium  magma  supply  (cross-sectional  areas  of  4.8  and  4.9  km^  respectively).  The 
quality  of  the  data  on  lines  SEPROl  through  SEPRIO  is  not  as  high  as  for  the  other  seismic 
lines  because  a  10  s  shooting  rate,  though  increasing  the  fold  of  the  data,  introduces  noise 
due  to  reverberation  of  the  previous  shot  in  the  water  colunm.  The  ridge  in  these  regions 
has  a  domed  but  narrower  axial  high  and  the  axial  depths  are  slightly  greater  than  in  Figure 
4a.  Note  the  small  off-axis  graben  on  SEPR70  and  the  axial  graben  on  SEPR06,  evidence 
for  an  increased  role  of  tectonism  relative  to  SEPR27.  The  on-axis  thickness  of  layer  2A 
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on  SEPR70  (220  m)  is  about  the  same  as  in  Figure  4a,  while  on  SEPR06  it  is  somewhat 
greater  (270  m).  The  off-axis  thicknesses  of  layer  2A  (510  and  520  m  respectively)  are 
again  about  twice  the  on-axis  values.  The  rate  of  2A  thickening  on  these  two  profiles  is 
somewhat  less  than  in  the  previous  case,  with  half-widths  of  thickening  of  3.4  and  3.2  km 
respectively.  The  magma  sill  on  SEPR70  is  relatively  shallow  (1080  m  below  the  seafloor) 
and  wide  (750  m),  while  on  SEPR06  it  is  deeper  (1340  m),  narrower  (-300  m),  and  is 
associated  with  a  weaker  reflection.  The  weak  reflection  on  the  across-axis  profile  may  be 
due  to  the  scattering  of  energy  from  seafloor  topography.  On  the  along-axis  line,  where 
topographic  effects  are  less  of  a  problem,  the  magma  sill  reflector  appears  quite  bright. 

Figure  4d  is  from  where  the  southern  segment  deepens  towards  the  20.7°S  OSC,  a 
region  with  an  extremely  magma-starved  morphology  (cross-sectional  area  =  0.6  km^, 
latitude  20°10'S).  The  ridge  is  significantly  deeper  (2846  m)  and  very  narrow  and 
triangular  in  shape.  On  this  profile  (SEPR02)  the  seafloor  shows  abundant  evidence  for 
tectonic  extension,  and  the  resulting  scattering  as  well  as  the  10  s  shooting  rate  on  this 
profile  reduce  the  quality  of  the  image  of  the  crustal  structure.  The  along-axis  data  reveal  a 
magma  sill  at  this  location,  but  it  is  so  weak  that  on  the  across-axis  profile  it  can  only  be 
seen  in  constant  velocity  stacks.  This  magma  sill  has  an  unusually  shallow  depth  given  the 
magma-starved  morphology.  On-axis  layer  2 A  is  thicker  on  this  profile  (330  m)  and 
thickens  more  rapidly  (half-width  of  thickening  =1.5  km)  than  on  other  profiles  to  attain  an 
off-axis  thickness  of  620  m. 

The  three  sections  of  ridge  defined  above  show  significant  differences  in  the  range 
and  variability  of  the  seismic  crustal  structure.  While  the  extrusive  layer  (seismic  layer  2A) 
has  a  comparatively  uniform  thickness  over  the  entire  region  (median  on-axis  thickness  is 
235  m;  range  200-300  m;  see  Figure  5a),  there  are  differences  in  the  thickness  and 
variability  of  layer  2A  along  the  ridge.  The  northern  section  has  a  very  uniform  2A 
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thickness  (222±21  m),  the  central  section  shows  more  variability  in  thickness  (233±48  m), 
while  the  southern  section  of  ridge  has  the  thickest  and  most  variable  extrusive  layer 
(284±56  m).  It  appears  that  some  thickening  of  layer  2A  occurs  towards  ridge  axis 
discontinuities  (shown  as  vertical  lines  on  Figure  3).  This  is  particularly  clear  near  the 
15°55'S  OSC. 

In  contrast  to  the  small  changes  in  extrusive  layer  thickness  along-axis,  the  depth  of 
the  magma  sill  is  quite  variable  along  the  study  area  (Figure  3a).  It  ranges  in  depth  from 
800  to  2000  m  below  the  seafloor  with  a  median  depth  of  1300  m  (Figure  5b).  The 
northern  and  central  ridge  sections  exhibit  a  general  along-axis  doming  of  the  magma  sill 
(shown  by  a  dashed  line  in  Figure  3a).  On  average,  the  magma  sill  is  deeper  near  the  ends 
of  the  ridge  sections  and  shallower  near  the  mid-point.  However,  there  are  also  local 
changes  in  magma  sill  depth  (-500  m)  superimposed  on  these  general  trends.  Along  the 
northern  ridge  section  the  depth  to  the  magma  sill  is  relatively  uniform  (1220±1 10  m).  The 
local  variability  in  magma  sill  depth  is  small  and  mirrors  the  relatively  uniform  water 
depths.  Along  the  central  section  of  ridge,  the  magma  sill  is  on  average  somewhat  deeper 
and  shows  a  much  greater  depth  range  (820  to  2000  m  with  a  median  depth  of  1350  m). 
The  shallowest  lens  is  observed  at  17°25'S  where  the  magma  sill  shoals  from  1350  m 
below  the  seafloor  to  820  m  over  a  distance  of  only  4  km  (Figure  4a).  Doming  of  the 
magma  sill  on  a  wavelength  of  tens  of  kms  (e.g.  between  17°  and  17°20'S)  may  be 
associated  with  individual  volcanic  constructs  [pers.  comm.,  Sinton,  1995].  Along  the 
southern  section  of  the  study  area  the  magma  sill  is  observed  more  sporadically.  The  depth 
range  is  similar  to  that  of  the  central  section  (1350±210  m)  but  changes  in  magma  sill  depth 
are  very  abrupt  and  show  no  consistent  pattern.  An  unusually  shallow  magma  sill  (896  m) 
is  imaged  at  20°15'S  (Figure  4d)  despite  the  magma-starved  ridge  morphology.  The 
magma  sill  appears  to  shoal,  rather  than  deepen,  towards  the  large  OSC  at  20.7°S. 
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Correlations  with  Indicators  of  Magma  Supply 


We  compare  the  variations  in  cmstal  stmcture  along  the  SEPR  described  above  with 
various  indicators  of  magma  supply  (Figure  3).  These  indicators  are;  the  depth  of  the 
seafloor  at  a  1  km  spacing  obtained  from  the  URl  Sea  Beam  data  synthesis  \Tighe  et  ai, 
1988];  cross-sectional  area  measured  at  1  km  intervals  along  the  ridge  from  Scheirer  and 
Macdonald  [1993];  the  along-axis  mantle  Bouguer  anomaly  calculated  at  2  km  intervals 
from  Magde  et  al.  [1995];  MgO  wt%  of  lavas  dredged  from  the  ridge  axis  from  Sinton  et 
al.  [1991]  and  samples  collected  from  the  youngest  on-axis  flows  by  the  submersible 
Nautile  [Auzende  et  al,  1996];  light  attenuation  in  the  water  column  integrated  from  100  to 
300  m  above  the  seafloor  in  bins  of  1.8  minutes  (~  18  km)  from  Baker  and  Urabe  [1996]; 
and  the  locations  along  the  ridge  where  Lonsdale  [1989]  detected  an  ASC  or  an  axial 
summit  graben  using  SeaBeam  data. 

Morphology  and  Gravity 

Axial  depth,  ridge  cross-sectional  area,  and  along-axis  mantle  Bouguer  anomaly 
are  all  well  correlated  (Figure  6).  The  ridge  axis  becomes  deeper  as  the  cross-sectional  area 
of  the  ridge  decreases  and  at  the  same  time  the  mantle  Bouguer  anomaly  increases 
indicating  higher  subsurface  densities.  The  ridge  topography  is  supported  by  a  low  density 
volume  underlying  the  ridge  crest,  which  is  partitioned  between  the  cmst  and  the  mantle 
[Wang  and  Cochran,  1993;  Magde  etal,  1995].  The  along-axis  variation  in  gravity  (range 
=  20  mGal)  and  part  of  the  axial  depth  changes  (range  =  400  m)  can  be  explained  by  a 
decrease  in  crustal  thickness  of  -500  m  between  18°S  and  20°S  [Cormier  and  Macdonald, 
1995].  The  variation  in  cross-sectional  area  (range  =  6  km^)  suggests  a  change  in  volume 
of  the  hot,  low  density  region  in  the  lower  cmst  accompanied  by  a  variation  in  the  density 
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anomaly  in  a  narrow,  deep  region  in  the  mantle  [Magde  et  al,  1995],  Thus,  as  is  predicted 
by  the  magma  supply  model,  along  axis  variations  in  axial  depth  and  ridge  cross-sectional 
area  are  indicative  of  variations  in  magma  supply  to  the  ridge. 

Basalt  Geochemistry 

The  variation  in  basalt  geochemistry  along  the  SEPR  is  shown  in  Figure  7  [Sinton 
etal,  1991].  The  MgO  content  of  basalts  reflects  the  amount  of  fractional  crystallization 
that  the  lavas  have  undergone  in  the  crust.  If  the  parental  compositions  are  similar,  MgO 
content  relates  to  eraption  temperature.  For  the  lavas  from  the  SEPR,  MgO  contents  range 
from  5  to  8.5  wt%  with  a  corresponding  change  in  Mg  number  from  0.35  to  0.64  (Figures 
7a&b).  Applying  the  olivine  thermometer  of  Sisson  and  Grove  [1993]  to  representative 
along-axis  major  element  compositions  [Sinton  et  al,  1991]  we  obtain  an  approximately 
linear  relationship  between  MgO  content.  Mg  number  and  lava  temperature  for  this  area 
(Figures  7a-c).  Emption  temperatures  range  from  1 100°C  to  i210°C.  The  MgO  content 
(and  thus  the  temperature)  of  the  erupted  lavas  correlates  well  with  cross-sectional  area 
(Figures  6c,  3b&c).  The  correlation  between  MgO  and  cross-sectional  area,  or 
alternatively  the  width  of  the  axial  high,  which  has  been  noted  previously  [Langmuir  et  al, 
1986;  Scheirer  and  Macdonald,  1993],  suggests  that  along  morphologically  inflated 
sections  of  the  ridge  the  empted  lavas  are  more  primitive  and  reside  in  the  crust  for  a 
shorter  time. 

The  variation  in  major  element  basalt  geochemistry  (Na^.o  and  Feg.o  -  Figure  7d&e) 
does  not  indicate  any  systematic  patterns  in  the  extent  or  depth  of  melting  along  the  SEPR. 
This  is  consistent  with  the  small  along-axis  variations  in  cmstal  thickness  inferred  from 
gravity  data  in  this  area  [Cormier  and  Macdonald,  1995;  Magde  et  al,  1995].  Diapiric 
feeding  of  fast  spreading  ridge  axes  has  recently  been  suggested  by  Nicolas  et  al  [1996] 
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and  Wang  et  al.  [1996].  This  melt  supply  model  would  predict  higher  Feg.o  and  lower 
Nag.o  near  segment  mid-points  due  to  a  greater  depth  and  extent  of  melting.  Feg.o  does  not 
vary  significantly  along  the  SEPR  indicating  that  the  depth  of  melting  is  relatively  uniform. 
There  is  more  variation  in  Nag.o  along  axis  suggesting  more  variability  in  the  extent  of 
melting.  However,  the  variations  in  Nag  o  do  not  correlate  with  variations  in  ridge  cross- 
sectional  area  or  mantle  Bouguer  anomaly.  For  example,  the  higher  Nag.o  values  observed 
along  the  inflated  ridge  section  between  14°S  and  15°S  are  opposite  of  what  would  be 
expected  if  a  mantle  diapir  with  greater  extent  of  melting  were  located  beneath  this  section 
of  ridge.  The  manner  in  which  melt  is  supplied  to  the  ridge  from  the  mantle  appears  to  be 
complicated.  However,  the  correlation  of  mantle  Bouguer  anomaly  variations  with  cross- 
sectional  area  and  axial  depth  suggests  that  on  a  broad  scale  variations  in  the  melt  supply  to 
the  SEPR  result  in  differences  in  the  mean  thermal  structure  and  morphology  of  the  ridge. 

Magma  Sill  Properties 

Purdy  et  al.  [1992]  has  proposed  that  on  average,  magma  sill  reflector  depth 
decreases  with  increasing  spreading  rate  beneath  fast  and  intermediate  spreading  centers. 
The  average  depth  of  the  magma  sill  is  indeed  greater  on  the  northern  EPR  than  on  the 
southern  EPR  (Table  1).  Spreading  rate  along  the  SEPR  varies  from  a  full  rate  of  146 
mm/yr  at  13.8°S  to  153  mm/yr  at  20°S  for  0  to  0.73  My  {Naar  and  Hey,  1989].  Along  the 
northern  EPR  spreading  rates  are  slower,  -112  mm/yr  [Carbotte  and  Macdonald,  1992]. 
The  difference  in  mean  magma  sill  depth  is  400  m  (1 140±150  m  to  1540±85  m  below  the 
seafloor  respectively),  which  is  greater  than  the  standard  deviation  of  the  averages.  On  the 
other  hand,  there  is  no  significant  difference  in  the  width  of  the  magma  sill  between  the 
northern  and  southern  EPR  (6801280  m  and  788±295  m  respectively.  Table  1.) 
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The  simple  magma  supply  model  predicts  that  for  a  given  spreading  rate,  variations 
in  magma  sill  depth  and  width  depend  on  changes  in  magma  supply  to  the  ridge. 
However,  along  the  SEPR,  where  spreading  rates  vary  only  slightly,  the  magma  sill 
properties  (depth  and  width)  are  not  well-correlated  with  axial  depth,  cross-sectional  area  or 
MgO  wt%  of  erupted  lavas  (Figure  8).  The  magma  sill  shows  a  large  range  of  depths 
along  the  SEPR  from  800-2000  m  (Figure  5b).  The  current  magma  supply  model  would 
predict  a  deeper  and  narrower  magma  sill  where  the  ridge  becomes  deeper  and  narrower 
[Macdonald  and  Fox,  1988;  Mutter  et  al.,  1995].  However,  we  find  that,  on  a  regional 
scale,  the  correlation  between  magma  sill  depth  and  water  depth  is  low  (Figure  8a).  This  is 
due  to  the  large  variability  in  magma  sill  depth  compared  to  the  relatively  uniform  water 
depth  in  the  region.  A  plot  of  magma  sill  depth  versus  cross-sectional  area  (Figure  8b)  also 
shows  a  large  amount  of  scatter.  Similarly,  lava  eraption  temperature  (MgO  content)  does 
not  correlate  with  magma  sill  depth  (Figure  8c).  It  has  been  suggested  that  magma  sill 
width  and  depth  are  correlated  and  vary  with  ridge  morphology,  and  that  deviations  from 
these  trends  indicate  recent  or  on-going  eruptions  [Mutter  et  al.,  1995].  However,  we  find 
that  magma  sill  widths  are  quite  variable  and  do  not  correlate  with  either  water  depth  or 
cross-sectional  area  (Figure  8d&e).  The  across-axis  lines  also  show  that,  on  the  regional 
scale,  there  is  a  great  range  in  magma  sill  width  for  any  given  magma  sill  depth  (Figure  8f). 
Thus  the  width  and  depth  of  the  magma  sill  are  extremely  variable  along  the  SEPR,  and 
show  little  relationship  with  morphological  indicators  of  magma  supply. 

The  observation  that  EPR  lavas  are  not  homogenized  has  been  interpreted  to  mean 
that  a  continuous  magma  chamber  can  exist  only  if  there  are  many  separate  regions  of 
supply  from  below  the  crust,  and  the  magma  chamber  is  not  well  mixed  [Langmuir  et  al, 
1986].  This  supports  the  rapidly  varying  nature  of  the  magma  sill  that  we  have  noted.  A 
magma  sill  which  is  not  fully  molten  [Hussenoeder  et  ai,  1996]  would  probably  be  less 
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well  mixed.  The  magmatic  system  at  fast  spreading  ridges  is  complicated,  consisting  of  a 
mostly  molten  magma  lens  or  sill  [Detrick  et  ai,  1987;  Harding  et  al.,  1989;  Kent  et  ai, 
1990]  overlying  a  broader  region  of  the  lower  cmst,  with  low  seismic  velocities  and  few 
percent  of  partial  melt  [Toomey  etal.,  1990;  Vera  et  ai,  1990;  Caress  et  al.,  1992].  There 
is  also  some  evidence  for  a  second  magma  sill  at  the  Moho  [Garmany,  1989;  Crawford  et 
al,  1996;  Kelemen  et  al,  1996].  In  addition,  it  been  suggested  that  the  lower  cmst  at  the 
EPR  may  be  composed  of  a  series  of  gabbroic  and  partially  molten  sills  [Kelemen  et  al, 
1996].  The  correlation  between  MgO  content,  mantle  Bouguer  anomaly,  and  ridge  cross- 
sectional  area  (Figure  6)  suggests  that  empted  melts  have  temperamres  representative  of  the 
lower  cmst.  This  implies  that  the  melts  that  empt  onto  the  seafloor  have  not  resided  very 
long  within  the  melt  lens.  Magmatic  events  (dike  intmsion  and  emptions)  may  thus  be 
closely  tied  to  replenishment  of  the  mid-cmstal  sill  from  a  lower-cmstal  or  sub-cmstal 
reservoir  [Kelemen  et  al,  1996)  Natland  and  Dick,  1996]. 

Extmsive  layer  (2A)  thickness 

While  the  magma  sill  varies  considerably  in  depth,  the  on-axis  layer  2 A  thickness  is 
relatively  constant  along  the  ridge  (Figure  9).  There  is  no  dependence  on  spreading  rate 
and  we  find  no  significant  difference  in  the  thickness  or  rate  of  thickening  of  the  extmsive 
layer  between  the  northern  and  southern  EPR  (Table  1).  The  high  correlation  coefficient 
between  on-axis  layer  2 A  thickness  and  water  depth  (Table  2)  is  due  to  the  uniform  value 
of  each  of  these  variables  and  is  not  very  informative.  The  on-axis  thickness  of  layer  2A  is 
constant  for  almost  the  entire  range  of  cross-sectional  areas.  However,  for  regions  with  an 
inferred  low  magma  budget,  i.e.  deep  ridges  with  small  cross-sectional  areas,  the  on-axis 
layer  2 A  thickness  appears  to  increase  (Figure  9a).  This  is  contrary  to  the  magma  supply 
model  which  predicts  a  thinner  extmsive  section  as  magma  supply  decreases.  Off-axis 
layer  2A  thickness  does  not  relate  to  cross-sectional  area  (Figure  9b)  and  seems  to  be 


58 


independent  of  the  magma  supply.  Both  on-  and  off-axis  layer  2A  thicknesses  are 
uncorrelated  with  magma  sill  depth  (Figures  9c&d).  The  width  of  the  zone  of  across-axis 
thickening  of  the  extrusive  layer  shows  a  large  range  (0.5  to  4.5  km)  (Figure  10),  but  there 
is  no  correlation  of  the  width  of  extrusive  thickening  with  magma  sill  width  (Figure  10a), 
and  only  a  weak  positive  correlation  with  ridge  cross-sectional  area  (Figure  10b).  In 
contrast  to  the  predictions  of  the  magma  supply  model,  the  extrusive  layer  thickness  shows 
little  correlation  with  either  morphological  indicators  of  magma  supply  or  magma  sill  depth. 

Hydrothermal  Activity 

Hydrothermal  activity  along  the  SEPR  is  not  strongly  correlated  with  indicators  of 
magma  supply  or  magma  sill  depth  (Figure  1 1).  Hydrothermal  activity  is  very  variable 
along  the  SEPR  and  is  particularly  vigorous  from  17°30'S  to  18°40'S  where  the  water 
column  survey  ended  [Baker  and  Urabe,  1996].  The  segment  centered  on  18°30'S  shows 
high  venting  activity  and  is  unusual  in  that  it  does  not  have  an  ASC  at  its  summit  but  rather 
a  summit  graben,  indicating  more  extensive  tectonic  activity  (across-axis  seismic  profile 
SEPR06  is  located  here.  Figure  4c).  In  a  study  of  the  distribution  of  hydrothermal  vents 
along  the  EPR,  Baker  [1996]  concludes  that  a  high  magma  supply  rate  and  an  observable 
magma  sill  are  necessary  but  not  deterministic  for  extensive  hydrothermal  venting. 
However,  our  observations  suggest  that  there  is  a  better  correlation  between  hydrothennal 
venting  and  regions  where  an  ASC  or  an  axial  graben  is  present  [Lonsdale,  1989]  than  with 
the  presence  or  depth  of  the  magma  lens  or  with  indicators  of  variations  in  magma  supply 
(Figure  3). 

It  is  generally  thought  that  hydrothermal  fields  begin,  or  are  renewed,  by  dike 
injection  events,  which  are  associated  with  both  heat  input  and  cracking  of  the  shallow 
crust  [Baker  et  al,  1987;  Haymon  et  al,  1993;  Embley  and  Chadwick,  1994;  Embley  et 
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a/.,  1995;  Wright  et  ai,  1995]  or  by  the  cooling  of  extremely  young  lavas  [Fomari  and 
Embley,  1995].  Charlou  et  al.  [1996]  and  Auzende  et  al.  [1996]  have  observed 
hydrothermal  vent  chimneys  predominantly  located  on  faults  at  the  edge  of  axial  summit 
graben  in  the  region  from  18°22'S  to  18°37'S.  In  contrast,  where  the  magma  sill  is 
observed  to  be  unusually  shallow,  at  17°25'S,  they  observed  mainly  diffuse  venting  of 
seawater  heated  by  cooling  lavas.  Our  observations  over  the  entire  region  suggest  that  at 
the  fast-spreading  SEPR  there  is  always  a  sufficient  heat  source  within  the  lower  crust  to 
drive  hydrothermal  circulation.  However,  a  permeability  structure  capable  of  focusing 
water  migration  through  the  crust  is  required  for  longer-lived  hydrothermal  venting.  This 
permeability  structure  can  be  generated  by  cracking  due  to  recent  diking,  or  by  tectonic 
fissuring.  The  relative  importance  of  these  two  processes  is  still  unclear. 

Discussion 

From  the  comparison  of  axial  seismic  structure  with  various  indicators  of  magma 
supply  along  the  SEPR  we  find  that;  (1)  as  expected  from  the  magma  supply  model,  axial 
depth,  cross-sectional  area  and  mantle  Bouguer  anomaly  are  well  correlated  with  one 
another  and  with  the  along-axis  variation  in  MgO  wt%  (Figure  6);  (2)  contrary  to 
predictions  of  the  simple  magma  supply  model,  magma  sill  properties  (depth  and  width)  are 
not  well-correlated  on  a  regional  scale  with  axial  depth,  cross-sectional  area  or  MgO  wt% 
of  erupted  lavas;  (3)  on-axis  layer  2A  thickness  is  relatively  constant  and  uncorrelated  with 
along-axis  variations  in  cross-sectional  area  or  magma  sill  depth  (Figure  9);  and  (4) 
hydrothermal  intensity  is  not  strongly  correlated  with  magma  sill  depth  or  ridge  cross- 
sectional  area  (Figure  11).  These  observations  lead  to  three  specific  questions  regarding 
the  temporal  and  spatial  processes  controlling  the  magmatic  construction  of  the  oceanic 
cmst:  Why  does  extrusive  thickness  not  seem  to  correlate  widi  morphological  indicators  of 
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magma  supply?  What  controls  the  depth  of  the  magma  sill?  Why  doesn't  magma  sill  depth 
correlate  with  morphological  indicators  of  magma  supply? 

Thickness  and  Construction  of  the  Extrusive  Layer 

The  uniform  along-axis  thickness  of  layer  2A  (Figure  5a),  relative  to  morphologic 
indicators  of  magma  supply  and  to  the  depth  of  the  magma  lens,  indicates  that  construction 
of  the  extrusive  section  is  controlled  by  eruption  dynamics  rather  than  magma  supply  (c.f. 
Kappus  et  al,  1995).  Models  of  the  emplacement  geometry  of  the  extrusive  section  show 
that  the  width  of  the  dike  intrusion  zone  and  the  length  of  lava  flows  control  the  thickness 
and  thickening  rate  of  layer  2A  [e.g.  Kidd,  1977;  Atwater,  1979;  Schouten  and  Denham, 
1979;  Pdlmason,  1980;  Hooft  et  al.,  1996].  The  relative  uniformity  of  on-axis  layer  2A 
thickness  along  the  ridge  axis  thus  indicates  that  the  dike  injection  zone  has  a  fairly  constant 
width  throughout  this  region. 

Near  ridge  offsets,  and  along  parts  of  the  ridge  with  small  cross-sectional  areas, 
seismic  layer  2A  thickens  (Figure  9a).  Near  OSCs  the  ridge  axis  is  not  well  localized  and 
the  dike  intrusion  zone  may  be  wider.  This  would  result  in  a  more  distributed  dike 
intrusion  and  thicker  average  extrusive  section  [c.f.  Hooft  et  al.,  1996].  An  alternative 
possibility  is  that  the  magma  supply  decreases  near  these  ridge  offsets,  and  increased 
tectonic  faulting  and  fissuring  reduces  the  surficial  seismic  velocities  of  the  crust. 
Interpretation  of  the  two-way  travel  time  thicknesses  without  accounting  for  these  reduced 
velocities  would  result  in  an  apparent  thickening  of  seismic  layer  2A,  which  does  not 
correspond  to  an  actual  thickening  of  the  extrusive  section . 

The  off-axis  thickness  and  thickening  of  seismic  layer  2A  is  more  variable  and  may 
be  controlled  by  the  shape  of  the  axial  high.  Emplacement  models  predict  that,  if  the  dikes 
are  rooted  in  the  magma  sill,  the  width  of  the  thickening  region  of  the  extrusive  section  is 
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related  to  magma  sill  width.  However,  our  observations  show  that  the  width  of  2A 
thickening  is  much  greater  and  more  variable  than  the  width  of  the  magma  sill  (Figure  10a). 
On  the  other  hand,  ridge  crest  topography  does  appear  to  control  the  construction  of  layer 
2A  to  some  extent.  The  width  of  the  thickening  zone  shows  a  weak  positive  correlation 
with  ridge  cross-sectional  area  (Figure  10b)  -  a  narrower  ridge  has  a  narrower  zone  over 
which  thickening  occurs.  These  observations  are  consistent  with  the  existence  of  a  narrow 
dike  intrusion  zone  and  lavas  that  flow  a  considerable  distance  off-axis  as  in  the  layer  2A 
construction  model  of  Hooft  et  al.  [1996]. 

Controls  on  the  Depth  of  the  Magma  Sill 

It  has  been  suggested  that  magma  ponds  at  a  level  of  neutral  buoyancy  (LNB)  in  the 
crust  where  the  melt  density  is  equal  to  the  bulk  density  of  the  surrounding  crust  [c.f. 
Walker,  1989;  Ryan,  1993].  At  mid-ocean  ridges  the  LNB  lies  at  an  approximately 
constant,  and  very  shallow  depth  below  the  seafloor  (-100-400  m)  [Hooft  and  Detrick, 
1993].  The  depth  of  the  seismically  observed  magma  lens  reflector  along  the  SEPR  is 
considerably  deeper  than  this  LNB.  The  median  depth  of  the  magma  sills  along  the  SEPR 
is  1300  m  below  the  seafloor,  and  even  the  shallowest  magma  sill  (800  m)  lies  well  below 
the  LNB  (Figure  5b). 

The  Hooft  and  Detrick  [1993]  calculation  of  the  depth  of  the  LNB  assumes  that  the 
magma  sill  is  completely  molten.  However,  if  the  magma  sill  contains  a  significant 
proportion  of  crystals,  it  will  be  neutrally  buoyant  at  a  greater  depth.  Detailed  waveform 
modeling  and  amplitude  versus  offset  analysis  of  multichannel  seismic  records 
[Hnssenoeder  et  al.,  1996]  show  that  in  some  areas  along  the  SEPR  the  magma  sill  is  not 
fully  molten  and  is  able  to  propagate  shear  waves,  though  the  shear  wave  velocity  is 
reduced.  Laboratory  measurements  suggest  that  these  lowered  shear  wave  velocities 
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correspond  to  a  melt  fraction  anywhere  between  30%  and  80%  [Murase  and  McBimey, 
1973].  Consistent  with  these  measurements,  numerical  models  of  the  observed 
compressional  wave  velocity  and  anisotropy  of  the  magma  sill  require  a  basaltic  melt 
fraction  of  50-70%  within  a  gabbro  matrix  [Mainprice,  1996].  To  calculate  the  density  of 
such  a  magma,  we  assume  (i)  that  magma  density  is  linearly  related  to  the  proportions  of 
crystals  and  melt,  (ii)  that  gabbro  crystal  densities  range  from  ~2.90-2.95  g/cm^ 
[Christensen  and  Salisbury,  1982],  and  (iii)  that  basaltic  melt  density  is  ~2.7  g/cm^  [Hooft 
and  Detrick,  1993].  Crystal  contents  of  30-50%  give  magma  densities  of  2.76-2.83  g/cm^ 
(Figure  12a).  For  a  seismically-derived  crustal  density  profile  [Christensen  and  Shaw, 
1970;  Harding  etal,  1993]  (Figure  12b),  such  a  magma  sill  is  neutrally  buoyant  between 
250  and  830  m  below  the  seafloor.  Thus,  even  if  the  magma  sill  has  a  significant  crystal 
content  it  still  lies  well  below  the  level  of  neutral  buoyancy,  although  some  more  evolved 
(heavier)  melts  may  lie  near  their  LNB. 

Another  line  of  evidence  indicating  that  the  magma  sills  do  not  pond  at  a  LNB  is  the 
large  range  of  magma  sill  depths  observed  along  the  SEPR.  Bulk  crustal  density  is  closely 
related  to  seismic  velocity.  Since  the  upper  cmstal  velocity  structure  is  relatively  invariant 
along  the  ridge  crest  [Tolstoy  et  ai,  1996],  the  upper  crustal  density  profile  remains 
approximately  uniform.  However,  the  depth  of  the  seismically  observed  magma  lens 
reflector  varies  considerably  ranging  from  800  to  2000  m  below  the  seafloor  (Figure  3). 
The  large  variability  of  magma  sill  depth  observed  along  the  SEPR  and  its  location  well 
below  the  LNB  are  inconsistent  with  the  hypothesis  that  the  magma  lens  forms  by  the 
ponding  of  melt  at  a  level  of  neutral  buoyancy  within  the  shallow  crust. 

We  believe  that  a  thermally  controlled  permeability  boundary,  such  as  the  solidus 
horizon,  controls  the  depth  of  the  magma  sill  at  the  EPR  [Hooft  and  Detrick,  1993;  Phipps 
Morgan  and  Chen,  1993a;  Natland  and  Dick,  1996].  The  thermal  structure  of  the  ridge  is  a 
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balance  between  melt  supplied  to  the  ridge  and  heat  removed  by  hydrothermal  cooling.  If 
magma  ponds  at  a  thermal  boundary,  the  observed  variations  in  magma  sill  depth  are  due  to 
changes  in  axial  thermal  stmcture.  This  concept  is  consistent  with  the  observed  decrease  in 
average  depth  of  the  magma  sill  with  increasing  spreading  rate  (Table  1,  Purdy  et  ai, 
1992). 

Lack  of  Correlation  of  the  Magma  Sill  Depth  with  Indicators  of  Magma  Supply 

Models  of  the  thermal  balance  at  mid-ocean  ridges  show  that  the  advection  of  heat 
from  variations  in  melt  supply  to  the  ridge  affects  the  crustal  thermal  structure  [Phipps 
Morgan  and  Chen,  1993b].  Excess  heat  from  an  increased  magma  budget  is  predicted  to 
result  in  a  shallower  magma  lens.  However,  the  lack  of  correlation  between  magma  sill 
depth  and  cross-sectional  area  along  the  SEPR  does  not  support  such  a  relationship.  The 
large  variability  in  magma  sill  depth  and  the  lack  of  correlation  with  morphologic  indicators 
of  magma  supply  (Figures  5  and  8)  require  temporal  and/or  spatial  variations  in  the  thermal 
strucmie  of  the  ridge  that  are  not  solely  dependent  on  spreading  rate  or  magma  supply. 

We  assume  that  the  magma  sill  lies  at  the  temperature  of  the  solidus  (~1200°C). 
The  lower  crust  is  assumed  to  be  relatively  isothermal  since,  to  satisfy  the  geophysical 
observations  [Toomey  etai,  1990;  Vera  etal.,  1990;  Caress  etai,  1992],  it  cannot  contain 
more  than  a  few  percent  partial  melt.  The  uppermost  crust  will  be  rapidly  cooled  by  the 
circulation  of  hydrothermal  fluids.  Temperamres  increase  rapidly  just  below  the  seafloor 
(0-350°C  in  the  first  few  hundred  meters),  but  below  this  the  upper  crust  remains  nearly 
isothermal  (~400°C)  to  the  maximum  depth  of  penetration  of  hydrothermal  fluids  [Nehlig 
and  Juteau,  1988;  Sleep,  1991].  There  is  a  conductive  lid  just  above  the  magma  sill  with 
very  steep  thermal  gradients,  from  400-1200°C  within  a  few  hundred  meters  [Lister,  1983; 
Nicolas,  1989].  Since  the  magma  sill  depth  ranges  from  800  to  2000  m  below  the 
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seafloor,  the  steep  thermal  gradient  and  the  depth  of  the  solidus  move  up  and  down 
considerably  within  the  mid  to  upper  crust  (Figure  13).  Rapid  variations  in  local  thermal 
structure  of  the  upper  crust  are  probably  due  to  an  imbalance  between  the  amount  of  crustal 
cooling  relative  to  the  heat  supplied. 

The  heat  supplied  to  the  magma  sill  from  below  is  probably  relatively  invariant, 
otherwise  the  thin  sill  would  quickly  be  frozen  by  hydrothermal  cooling.  Based  on  the  time 
required  to  construct  the  axial  high  area,  along-axis  variations  in  magma  supply  are 
probably  maintained  on  order  100  Kyr  [Scheirer  and  Macdonald,  1993].  Crustal  thickness 
[Barth  and  Mutter,  1996]  and  abyssal  hill  characteristics  [Gojf,  1991]  also  indicate  that  the 
average  rate  of  magma  supply  changes  little  over  several  100  Kyr.  A  record  of  the  thermal 
history  of  axial  magma  chambers  can  be  deduced  from  the  variations  in  chemistry  of 
volcanic  rocks  collected  along  flowlines  on  the  flanks  of  the  EPR  (9°30'N,  10°30'N  and 
11°20'N)  [Batiza  et  al.,  1996].  These  observations  show  that  where  the  ridge  is 
magmatically  robust  the  magma  lens  is  steady  state  and  there  are  no  significant  changes  in 
eruption  temperature  over  -800  Kyr.  On  the  other  hand,  in  magmatically  starved  areas, 
there  are  large  variations  in  eruption  temperature  with  time  (~160°C  over  600  Kyr) 
suggesting  that  the  magma  lens  departs  from  steady  state.  These  observations  are 
supported  by  a  study  of  lower  crustal  rocks  exposed  in  fracture  zones  and  rift  deeps  along 
the  fast-spreading  ridge  system.  Gabbro  samples  from  the  ends  of  segments  have  more 
evolved  compositions  than  those  from  more  stable  portions  of  the  ridge  [Constantin  et  al, 
1996]  indicating  less  frequent  replenishment  of  the  magma  system.  Relative  to  the  more 
constant  supply  of  melt  from  below  (-100  Kyr),  the  variations  in  the  amount  of  crustal 
cooling  may  depend  on  the  time  that  has  elapsed  since  a  diking  event  drained  the  magma  sill 
(1-10  yrs),  as  well  as  on  the  strength  of  hydrothermal  cooling  in  the  overlying  crust  (10s 
-100s  of  yrs). 


65 


The  cartoon  in  Figure  13  summarizes  two  possible  different  states  of  the  magma  sill 
for  two  parts  of  the  ridge  that  have  the  same  magma  budget,  i.e.  the  amount  of  melt 
supplied  from  below  is  constant.  Also  shown  are  non-steady-state  thermal  profiles 
resulting  from  temporal  variations  in  the  amount  of  cooling,  due  to  interaiittent  emptive  and 
tectonic  events  at  the  ridge  crest,  relative  to  the  more  constant  magma  supply  from  below. 
In  the  first  case,  the  crustal  permeability  is  high  allowing  deep  penetration  of  hydrothermal 
fluids.  The  magma  lens  has  also  recently  empted  onto  the  seafloor  thereby  draining  away 
the  melt.  The  amount  of  heat  removed  from  the  magma  sill  is  greater  than  that  being 
supplied  and  the  sill  lies  fairly  deep  within  the  crust.  This  scenario  represents  across-axis 
seismic  profile  SEPR06  (Figure  4c)  where  the  presence  of  an  axial  summit  graben 
enhances  the  penetration  of  hydrothermal  fluids  and  there  is  extensive  hydrothermal 
activity.  The  resulting  magma  sill  lies  deep  within  the  cmst.  On  the  other  hand,  if  the 
crustal  permeability  becomes  reduced  and  no  eruption  has  occurred  for  some  time,  then  the 
magma  sill  will  slowly  bum  its  way  to  shallower  levels  in  the  cmst,  perhaps  by  stoping  of 
the  overlying  sheeted  dike  section  [Nicolas  and  Boudier,  1991].  This  scenario  is 
representative  of  across-axis  seismic  profile  SEPR70  (Figure  4b),  where  the  magma  sill 
lies  much  shallower  than  on  SEPR06  even  though  both  ridge  segments  have  similar  cross- 
sectional  areas,  and  thus  their  longer  term  magma  budgets  are  the  same.  The  lower  panel 
of  Figure  13  shows  apparent  local  variations  in  the  depth  of  the  magma  sill  along  the  ridge 
axis;  a  shallow  and  almost  horizontal  magma  sill  fades  over  a  distance  of  ~7  km  giving  way 
to  a  bright  magma  sill  reflector  lying  ~1  km  deeper. 

The  relationship  between  these  local,  short-term  variations  in  cmstal  stmcture  and 
longer-term  changes  in  magma  supply  to  the  ridge  is  shown  in  an  along-axis  view  in  Figure 
14.  Over  time  scales  of  -100  Kyr  there  are  variations  in  magma  supply  to  the  different 
parts  of  this  section  of  ridge.  The  resulting  variations  in  cross-sectional  area  are  closely 
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mirrored  in  the  along-axis  patterns  of  the  mantle  Bouguer  anomaly  and  MgO  content  of  the 
lavas.  However,  the  magma  sill  depth,  extrusive  thickness  and  locations  of  hydrothermal 
venting  show  little  correlation  with  these  patterns.  Changes  in  average  magma  sill  depth 
are  on  a  scale  of  -200  km.  Short  wavelength  doming  of  the  magma  sill  on  the  scale  of  -50 
km  may  be  associated  with  individual  volcanic  constructs.  The  extrusive  section  thickens 
slightly  at  the  ends  of  these  constructs.  The  white  region  around  the  magma  sill  indicates  a 
mixed  zone  of  dikes  and  gabbros  generated  by  rapid  (10- 100s  yr),  local  changes  in  the 
depth  of  the  magma  sill  (Figure  13).  Hydrothermal  venting  is  diffuse  where  there  are 
recent  eruptions,  and  focused  where  there  is  an  ASC  or  axial  graben. 

Summary 

The  current  magma  supply  model  predicts  that  broad,  shallow  ridge  segments 
reflect  a  high  magma  supply  and  have  larger,  more  robust  magma  chambers,  erupt  more 
primitive  basalts,  form  thicker  crust,  and  have  more  vigorous  hydrothermal  activity  than 
sections  of  the  East  Pacific  Rise  with  deeper,  narrower  axial  highs.  To  test  this  paradigm 
we  examine  the  correlation  along  >800  km  of  the  SEPR  between  along-axis  variations  in 
seismic  crustal  stmcture  and  other  indicators  of  magma  supply.  Our  main  conclusions  are: 

•  Good  correlations  exist  between  axial  depth,  ridge  cross-sectional  area,  mantle 
Bouguer  anomaly,  and  the  MgO  wt%  of  basalts  recovered  from  the  rise  axis.  These 
correlations  suggest  along-axis  variations  in  lower  crustal  and  upper  mantle  density  which 
are  consistent  with  variations  in  magma  supply  on  time  scales  of  -100  Kyr. 

•  The  magma  sill  properties  (depth  and  width  of  the  melt  lens),  seismic  layer  2A 
thickness,  and  extent  of  hydrothermal  venting  are  poorly  correlated  with  regional  variations 
in  ridge  depth  and  cross-sectional  area. 
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•  The  emplacement  geometry  (width  of  the  intrusion  zone,  flow  lengths),  not 
magma  supply,  controls  extrusive  layer  (seismic  2A)  thickness. 

•  Hydrothermal  activity  is  closely  linked  to  spreading  events  (dike  intrusion, 
eruptions,  faulting)  which  occur  on  short  time  scales  (ca.  10-100  yr).  The  presence  of  a 
permeability  structure  that  allows  water  to  penetrate  to  hot  rock  is  key. 

•  Magmas  accumulate  at  a  depth  in  the  shallow  crust  where  their  buoyant  ascent  is 
arrested  by  a  thermally-controlled  permeability  boundary. 

•  On  average  magma  sill  depth  depends  on  spreading  rate  and  magma  supply. 

•  We  observe  a  large  local  variability  in  magma  sill  depth  (800-2000  m  below  the 
seafloor),  suggesting  that  on  the  short  term  the  thermal  structure  of  the  crust  is  strongly 
modified.  The  advection  of  heat  to  the  crust  by  the  ascent  of  melt  is  relatively  continuous. 
In  contrast,  the  removal  of  heat  from  the  crust  is  highly  variable  and  depends  on  the 
variability  of  hydrothermal  cooling,  on  intermittent  eruptive  and  intrusive  episodes  that 
drain  the  magma  sill,  and  on  tectonic  events  that  alter  the  cmstal  permeability  structure. 
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Table  1.  Ridge  morphology,  magma  sill  and  Layer  2A  characteristics 

*  From  the  along  axis  data 

2  Excluding  CDP  33.  Including  CDP33:  145011500 

*  Interpolated  from  Scheirer  and  Macdonald  [ 1 993] 

^  I  picked  and  calculated  from  Kent  et  al.  [1994] 

*  I  picked  and  calculated  from  Harding  et  al.  [1993] 

^Kentetal.  [1994] 

^  Harding  et  al.  [1993] 

^Kentetal.  [1993b] 

<^Kentetal.  [1993a] 
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Table  1:  Ridge  morphology,  magma  sill  and  Layer  2A  characteristics 


Latitude 

line  # 

water 

depth 

cross- 

sectional  _ 

magma  sill 

idtlCd 

Layer  2A  Thickness 

(m) 

area* 

(km^) 

width 

(m) 

depth 

(m) 

On-axis 

(m) 

Off-axis 

(m) 

Half-width 

(km) 

-20.16 

SEPR02 

2846 

Southern  East  Pacific  Rise 
0.6  10601 

330 

620 

1.5 

-19.44 

SEPR04 

2782 

2.3 

560 

1110 

240 

410 

0.7 

-18.46 

SEPR06 

2619 

4.9 

300 

1340 

270 

520 

3.2 

-18.21 

SEPR08 

2673 

5.0 

590 

1550 

200 

350 

1.2 

-17.54 

SEPR34 

2611 

5.2 

740 

1210 

260 

460 

4.1 

-17.46 

SEPR26 

2590 

5.1 

1540 

1000 

230 

390 

2.0 

-17.44 

SEPR33 

2588 

5.4 

950 

1190 

247 

530 

2.3 

-17.42 

SEPR27 

2591 

5.6 

850 

760 

230 

570 

3.0 

-17.38 

SEPR28 

2611 

5.5 

710 

1100 

240 

490 

2.1 

-17.33 

SEPR32 

2621 

5.0 

375 

1350 

210 

530 

1.3 

-17.29 

SEPR29 

2621 

5.4 

875 

1140 

260 

510 

1.5 

-17.24 

SEPR31 

2610 

6.2 

710 

1110 

250 

530 

1.3 

-17.20 

SEPR30 

2605 

5.8 

560 

1150 

270 

560 

1.6 

-17.16 

SEPRll 

2601 

5.3 

500 

1210 

260 

450 

1.9 

-17.07 

SEPRIO 

2628 

5.4 

1010 

1260 

240 

480 

2.2 

-14.49 

SEPR57 

2631 

4.3 

650a 

1100¥ 

240V 

570a 

1.8* 

-14.45 

SEPR68 

2628 

4.2 

400 

lOOOV 

270V 

540* 

1.5^ 

-14.41 

SEPR56 

2630 

4.4 

700 

990¥ 

190V 

540* 

2.3a 

-14.37 

SEPR67 

2634 

4.5 

1050 

940¥ 

270V 

540* 

2.7a 

-14.33 

SEPR69 

2636 

4.9 

950 

960V 

280V 

570* 

2.5a 

-14.29 

SEPR66 

2640 

5.0 

750 

1040V 

210V 

530* 

2.8a 

-14.25 

SEPR70 

2642 

4.8 

750 

1080^ 

220V 

510* 

3.4a 

-14.21 

SEPR65 

2638 

4.4 

700 

llOOV 

240V 

530* 

3.2a 

-14.17 

SEPR71 

2635 

4.6 

750 

1220^ 

270V 

540* 

3.4a 

-14.13 

SEPR64 

2636 

4.2 

375* 

1250V 

240V 

510* 

2.9a 

-14.09 

SEPR72 

2635 

4.7 

550 

1120V 

250V 

540* 

2.4a 

-14.05 

SEPR63 

2634 

4.8 

550 

1220V 

210V 

530* 

2.8a 

-14.01 

SEPR46 

2636 

4.9 

400 

1220V 

270V 

550* 

2.7a 

-13.81 

SEPR61 

2640 

4.5 

690 

1200 

200 

520 

1.9 

mean 

2640 

4.7 

680 

1140 

245 

515 

2.3 

Standard  deviation 

50 

1.1 

260 

150 

30 

57 

0.8 

+08.92 

CDP37 

2580 

Northern 

3.5 

East  Pacific  Rise 
750<<  1660<i 

3801 

560b 

2.1b 

+09.31 

CDP33 

2590 

3.2 

4100 

1500c 

2001 

510b 

1.6b 

+09.49 

CDP31 

2580 

3.6 

1200 

1560<= 

1801 

450b 

4.1b 

+09.66 

CDP29 

2560 

3.6 

700 

1550c 

2301 

440b 

2.3b 

+09.83 

CDP27 

2510 

4.0 

500 

1430c 

2801 

460b 

2.5b 

mean 

2570 

3.6 

7882 

1540 

254 

484 

2.5 

Standard  deviation 

40 

0.3 

2952 

85 

80 

50 

0.9 

79 
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Table  2:  Correlation  coefficients  for  across-axis  and  along-axis  data  from 

the  southern  EPR 


a)  Correlation  coefficients  for  the  Southern  East  Pacific  Rise:  N=29 


water 

cross- 

sectional 

magma  sill 

Layer  2AThickness 

depth 

area* 

width  depth 

On- 

axis 

Off-axis  Half-width 

water  depth 

LOO 

0.93 

0.21  0.11 

-0.36 

-0.18  0.23 

cross-sectional 

0.93 

1.00 

0.26  -0.03 

-0.33 

-0.06  0.35 

area* 

magma  width 

0.21 

0.26 

1.00  -0.45 

-0.03 

-0.20  0.08 

sUl  depth 

0.11 

-0.03 

-0.45  1.00 

-0.17 

-0.44  -0.12 

Layer  2A  On-axis 

-0.36 

-0.33 

-0.03  -0.17 

1.00 

0.40  0.02 

Off-axis 

-0.18 

-0.06 

-0.20  -0.44 

0.40 

1.00  0.20 

Half- 

0.23 

0.35 

0.08  -0.12 

0.02 

0.20  1.00 

width 

b)  Correlation  coefficients  for  the  Along-Axis  data:  N 

>4071 

water 

cross- 

sectional 

magma  sill 

Layer  2A  Thickness 

depth 

area* 

width  depth 

On-axis 

Off-axis  Half-width 

water  depth 

1.00 

0.79 

-0.05 

-0.96 

cross-sectional 

0.79 

1.00 

0.11 

-0.70 

area* 

magma  width 

sill  depth 

-0.05 

0.11 

1.00 

0.34 

Layer  2 A  On-axis 

-0.96 

-0.70 

0.34 

1.00 

Off-axis 

Half¬ 

width 

*  Interpolated  from  Scheirer  and  Macdonald  [1993] 
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Figure  1.  The  current  paradigm  for  magmatic  activity  on  fast  spreading  ridges  predicts  a 
contrast  in  ridge  morphology  and  crustal  structure  between  ridge  sections  with  a  high 
magma  budget  and  those  that  have  a  low  magma  budget 
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SIMPLE  MAGMA  SUPPLY  MODEL 
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-  ASC  present  .  hydrothermal  venting 

-  hydrothermal  venting 


Figure  2.  Map  of  the  southern  East  Pacific  Rise  (SEPR).  Major  segment  boundaries  are 
indicated.  A  multichannel  seismic  line  was  shot  along  the  ridge  crest  form  the  northern 
limb  of  the  20.7°S  propagator  to  the  Garrett  fracture  zone.  Boxes  indicate  the  intensive 
multichannel  seismic  study  areas  at  14°S,  15°55'S  and  17°S.  The  remaining  five  across- 
axis  seismic  lines  are  shown  by  dark  lines,  as  well  as  the  location  of  SEPR27  and  SEPR70 
in  the  intensive  study  areas  at  17°S  and  14°S  respectively. 
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Figure  3.  Along-axis  summary  figure  of  the  seismic  stmcture  and  other  indicators  of 
magma  budget  along  the  SEPR  from  the  20.7°S  propagator  to  the  Garrett  fracture  zone 
(13.4°S),  >800  km.  Vertical  lines  show  overlapping  spreading  centers  (OSCs).  The 
northern,  central  and  southern  ridge  sections  are  defined  based  on  differences  in  the  range 
and  variability  of  the  seafloor  depth,  extrusive  thickness,  and  magma  sill  depth,  and  are 
bounded  by  the  Garrett  transform,  the  15°55'S  and  17°55'S  OSCs,  and  the  large  20.7°S 
OSC.  (a)  seafloor  topography  (light  grey)  [Tighe  et  al,  1988],  depth  to  the  base  of  the 
extrusive  layer  (layer  2A)  (dark  grey),  and  the  magma  sill  reflector  (black).  Note  the 
relatively  uniform  thickness  of  layer  2A  and  the  large  viuiability  in  magma  sill  depth.  The 
general  variation  of  magma  sill  depth  along  three  sections  of  the  ridge  is  outlined  (dashed 
line).  Lxx^ations  of  the  across-axis  seismic  profiles  shown  in  Figure  4  are  indicated  in  light 
grey  at  the  bottom  of  the  figure,  (b)  cross-sectional  area  of  the  ridge  (dark)  [Scheirer  and 
Macdonald,  1993]  and  the  mantle  Bouguer  anomaly  (light  grey)  [Magde  et  al,  1995]. 
Note  the  excellent  correlation  between  these  two  measures  of  the  volume  of  the  crustal  low 
velocity  zone,  (c)  MgO  wt%  for  basalt  samples  dredged  from  the  ridge  axis  (dark  stars) 
[Sinton  et  al,  1991].  T-MORBs  are  indicated  by  dark  circles.  Samples  collected  from  the 
youngest  on-axis  flows  on  Nautile  submersible  dives  (light  circles)  [pers.  comm.  Sinton, 
1995]  do  not  deviate  much  from  the  dredge  samples  collected  in  the  same  area.  The  curve 
represents  a  mnning  mean  made  through  the  data  to  show  the  local  trend  (variability  is  ~±1 
wt%).  The  scale  on  the  left-hand  side  shows  the  corresponding  lava  temperatures.  The 
letters  indicate  the  nomenclature  for  different  magmatic  segments  [Sinton  etal,  1991].  (d) 
Light  attenuation  in  the  water  column  integrated  from  100  to  300  m  above  the  seafloor  in 
bins  of  0.03  degrees  [Baker  and  Urabe,  1996].  Along  the  x-axis  the  dark  bars  show  the 
portions  of  the  ridge  along  which  Lonsdale  [1989]  located  an  ASC  or  axial  graben  using 
SeaBeam  data.  Note  the  large  variability  in  the  location  of  hydrothermal  fields  and  the 
correlation  with  the  presence  of  an  ASC. 
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Along  Axis  Seismic  Profile 
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Latitude,  South 


Figure  4.  Along-  and  across-axis  seismic  lines  from  four  areas  of  the  ridge  that  range  from 
inflated  to  magma-starved.  The  across-axis  data  have  had  prestack  DMO  applied,  then 
were  stacked  and  finite  difference  migrated  to  correct  for  seafloor  and  layer  2A  topography. 
Note  the  range  in  widths  of  the  ridge  crest  and  amount  of  topography.  Also  the  differences 
in  the  thickening  profiles  of  layer  2A.  The  along-axis  lines  show  the  stacked  multichannel 
seismic  section.  Note  the  uniformity  in  the  depth  and  nature  of  the  layer  2A  reflection 
everywhere  along  the  axis.  In  contrast  there  are  significant  along-axis  differences  in  the 
depth  of  the  magma  sill  reflector,  (a)  Across-axis  profile  SEPR27  at  17°25'S  and  along- 
axis  profile  SEPR35  from  17°32'S  to  17°15'S.  The  ridge  is  shallow  and  has  a  wide  domed 
morphology.  At  17°24'N,  just  north  of  the  shallowest  and  broadest  portion  of  the  ridge, 
the  melt  lens  shoals  from  an  average  depth  of  1 150  m  to  800  m  below  the  seafloor  over  a 
distance  of  4  km.  (b)  Across-axis  profile  SEPR70  at  14°15'S  and  along-axis  profile 
SEPR45  from  14°23'S  to  14°07'S.  The  magma  sill  is  shallow  and  relatively  uniform  in 
depth  (1080  m  below  the  seafloor).  Note  the  graben  on  west  side  of  the  ridge  on  the 
across-axis  profile,  (c)  Across-axis  profile  SEPR06  at  18°28'S  and  along-axis  profiles 
SEPR05  and  SEPR07  from  18°37'S  to  18°2rS.  Note  the  axial  summit  graben  on  the 
across-axis  profile.  The  magma  sill  Ues  at  a  deeper  depth  (1340  m  below  the  seafloor)  and 
deepens  towards  the  OSCs  at  the  18°39'S  and  18°21'S.  Note  that  SEPR70  and  SEPR06 
have  similar  cross-sectional  areas  and  yet  there  is  a  large  difference  in  magma  sill  depth 
between  the  two.  (d)  Across-axis  profile  SEPR02  at  20°10'S  and  along-axis  profiles 
SEPROl  and  SEPR03  from  20°17'S  to  20°01'S.  The  ridge  is  deep  and  very  narrow  in  this 
region  and  there  is  evidence  for  significant  tectonism  on  the  flanks.  From  20°10-15'S  on 
the  along-axis  line  a  small,  and  surprisingly  shallow  (1060  m  below  the  seafloor)  magma 
lens  is  imaged  for  the  magma-starved  morphology. 
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SEPR  27: DM0  and  FDMIGR 


90 


SEPR  70:FDMICR 
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SEPR  2:DM0  and  FDMIGR 

SEPR  1 


Figure  5.  Histograms  of  2 A  and  magma  sill  depth  from  the  SEPR.  Data  is  binned  at  100 
m.  This  emphasizes  the  uniformity  of  layer  2A  thickness  on  the  ridge  crest  and  the  much 
larger  variabihty  of  magma  sill  depth. 
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Figure  6.  Plots  showing  water  depth,  mantle  Bouguer  anomaly,  and  MgO  content  of 
dredged  lavas  for  the  SEPR  plotted  as  a  function  of  cross-sectional  area.  Correlation 
coefficients  for  these  variables  are  shown  on  the  figure.  Variations  in  magma  supply  are 
accompanied  by  crustal  thickness  changes  of  -500  m  and  are  maintained  over  -100  Kyr . 
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Figure  7.  Major  element  geochemist^  of  lavas  dredged  from  the  ridge  axis  (representative 
glass  analyses  from  Table  1  -  Sinton  et  ai,  1991).  OSCs  are  shown  by  black  lines  as  in 
Figure  3.  a)  MgO  content  in  wt%.  b)  Mg  number  =  (Mg/Mg+Fe  molar),  c)  Lava 
temperature  calculated  from  the  major  element  composition  using  the  oUvine  thermometer 
of  Sisson  and  Grove  [1993]  assuming  atmospheric  pressure  and  no  water  content,  d&e) 
Nag.o  and  Feg.o  calculated  using  the  regression  of  Sinton  et  al.  [1991]  (Table  2). 
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a) 


b) 
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Figure  8.  Relationships  between  indicators  of  magma  budget  and  magma  sill  depth  and 
width.  All  the  along-axis  data  points  are  shown  as  dots,  the  measurements  based  on  the 
across-axis  lines  are  shown  as  grey  diamonds.  The  correlation  coefficients  for  these 
relationships  are  given  in  Table  2.  Note  the  large  variability  in  magma  sill  depth  which 
does  not  correlate  with  indicators  of  magma  budget.  The  lack  of  correlation  between  MgO 
of  the  lavas  and  magma  sill  depth  indicates  that  the  residence  time  of  the  magmas  in  the  sill 
must  be  short. 
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Figure  9.  Relationship  between  indicators  of  magma  budget  and  layer  2A  thickness  on- 
and  off-axis.  The  symbols  and  correlation  coefficients  are  given  as  in  Figure  8.  The 
correlation  coefficients  for  these  relationships  are  given  in  Table  2.  The  on-axis  layer  2A 
thickness  is  uniform  except  for  an  increase  in  thickness  where  the  ridge  cross-sectional  area 
is  low.  The  off-axis  thickness  of  layer  2A  is  not  related  to  morphological  indicators  of 
magma  supply. 
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Figure  10.  Relationship  between  on  and  off-axis  thickness  of  layer  2A  and  the  width  of  the 
region  of  thickening,  based  on  the  across  axis  seismic  lines,  and  indicators  of  magmatic 
budget.  The  symbols  are  as  in  Figure  8.  The  correlation  coefficients  for  these 
relationships  are  given  in  Table  2.  The  width  of  layer  2A  thickening  may  increase  with 
cross-sectional  area. 
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Figure  11.  Relationship  between  light  attenuation  in  the  water  column,  an  indicator  of 
hydrothermal  venting,  ridge  cross-sectional  area,  and  magma  sill  depth  along  the  ridge 
crest.  There  is  little  correlation  between  these  parameters. 
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Figure  12.  a)  Contours  of  crystal  content  (%)  as  a  function  of  the  density  of  the  crystals 
and  the  density  of  the  magma  sill.  We  fix  the  melt  density  at  2.7  g/cm^,  though  this  value 
may  lie  between  2.68  and  2.73  g/cm^  {Hooft  and  Detrick,  1993].  Crystal  density  ranges 
from  2.8-3.0  g/cm^,  most  gabbro  crystals  lie  between  2.90  and  2.95  g/cm?  [Christensen 
and  Salisbury,  1982].  We  assume  a  linear  relationship  between  the  density  of  the  magma 
sill,  Pms,  and  the  crystal  and  melt  densities,  and  p„  respectively,  as  a  function  of 
percentage  of  crystals  present, 

P«.=P»«  +  -^XP. 

b)  Density  structure  for  the  crust,  Pcmsu  determined  from  seismic  measurements  of 
compressional  wave  velocity,  Vp,  [Harding  etai,  1993]  using  the  empirical  relationship  of 
Christensen  and  Shaw  [1970] 

Pc™.,=l-85  +  0.165y^ 

For  a  fully  molten  magma  sill  the  level  of  neutral  buoyancy  (LNB)  is  l(X)-400  m  below  the 
seafloor  [Hooft  and  Detrick,  1993].  For  a  crystal  content  of  30-50%,  the  LNB  is  250-830 
m  below  the  seafloor.  The  observed  magma  sill  reflectors  lie  deeper  (800-2000  m  below 
the  seafloor)  than  the  LNB. 

c)  Flistogram  of  the  distribution  of  magma  sill  densities  calculated  from  the  measured 
magma  sill  depths  (Figure  5b)  assuming  that  they  lie  at  neutral  buoyancy  and  the  crustal 
density  structure  in  (b).  This  requires  a  magma  sill  crystal  content  of  80%  to  completely 
solid  which  is  too  high  for  the  observed  seismic  properties  of  the  magma  sill. 
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Figure  13.  Upper  Panel:  Cartoon  showing  two  ridge  crest  thermal  stmctures,  and  the 
accompanying  crustal  stmctures;  the  extmsives  and  sheeted  dikes,  the  depth  of 
hydrothermal  venting,  and  the  magma  sills.  The  depth  of  the  magma  sill  is  determined  by  a 
thermally  controlled  permeability  boundary.  The  thermal  stmcture  of  the  cmst  depends  on 
the  balance  of  heat  advected,  by  the  supply  of  melt,  and  cooling,  by  the  removal  of  melt, 
conduction,  and  hydrothermal  convection.  The  two  scenarios  shown  are  for  a  constant 
magma  supply.  Hydrothermal  cooling  is  very  sensitive  to  the  permeability  stmcture  of  the 
cmst  and  fluctuates  due  to  intermittent  emptive  and  tectonic  events.  Hydrothermal  cooling 
maintains  the  upper  cmst  at  ~400°C  down  to  the  depth  of  penetration  of  fluids.  Below  this 
temperature  increases  very  rapidly  to  reach  the  solidus  at  the  magma  sill.  The  lower  crust  is 
also  relatively  isothermal.  On  the  left,  heat  removal  by  hydrothermal  cooling  and  the  recent 
emption  of  the  magma  sill  exceeds  the  heat  supplied  from  beneath.  The  isotherms  are 
moving  down  in  the  cmst  leading  to  a  deeper  magma  sill.  On  the  right,  there  have  been  no 
recent  emptions  and  hydrothermal  cooling  is  weak.  Heat  supplied  from  below  exceeds  that 
removed  from  above  and  the  magma  sill  is  burning  its  way  up  to  shallower  levels  in  the 
cmst.  These  scenarios  simulate  the  observations  on  SEPRC^  (Figure  4c)  and  SEPR70 
(Figure  4d),  both  of  which  have  similar  magma  budgets  as  indicated  by  their  cross- 
sectional  areas.  However,  on  SEPR06  the  magma  sill  is  deep  and  has  a  weak  reflection 
and  there  is  an  axial  summit  graben  with  related  hydrothermal  venting.  On  the  other  hand, 
on  SEPR70,  the  magma  sill  is  shallow  and  bright  and  there  is  no  evidence  for  recent 
tectonic  activity.  Lower  Panel:  Along-axis  seismic  section  SEPR09  from  18°13'S  to 
18°02'S  illustrating  ^parent  local  depth  variations  of  the  magma  sill.  Note  how  the  magma 
sill  is  almost  horizontal,  however,  over  ~7  km  the  shallower  magma  sill  fades  giving  way 
to  a  bright  magma  sill  reflector  lying  ~1  km  deeper. 
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Figure  14.  An  along-axis  view  showing  the  relationship  between  the  local,  short-term 
variations  in  crustal  structure  in  Figure  13  and  longer-term  changes  in  magma  supply  to  the 
ridge.  The  upper  panel  shows  the  gabbro  section  (hatched  pattern),  magma  sill  (thick, 
double  daik  line),  sheeted  dike  section  (vertical  lines),  and  extmsive  section  (horizontal  zig¬ 
zag  pattern).  The  white  region  around  the  magma  sill  indicates  a  mixed  zone  of  dikes  and 
gabbros  generated  by  short-term  depth  changes  of  the  magma  sill  (Figure  13). 
Hydrothermal  venting  (wavy  pattern)  is  diffuse  where  there  are  recent  eruptions  (1-10  yrs), 
and  focused  and  longer-lived  (10- 100s  yrs)  where  there  is  an  axial  summit  caldera  (ASC). 
Short  wavelength  doming  of  the  magma  sill  on  the  scale  of  ~50  km  may  be  associated  with 
individual  volcanic  constmcts.  The  extrusive  section  thickens  slightly  at  the  ends  of  these 
units.  The  longer-term  changes  in  average  magma  sill  depth  are  on  a  larger  scale  (-200 
km).  The  lower  portion  of  the  figure  indicates  the  variation  in  cross-sectional  area  along 
this  portion  of  the  ridge  axis.  General  shoaling  of  the  magma  sill  is  accompanied  by  a  trend 
of  increasing  cross-sectional  area,  but  in  detail  the  patterns  are  different. 
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Scales  of  Variability 
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Chapter  4 


The  Role  of  Density  in  the  Accumulation  of  Basaltic 
Melts  at  Mid-Ocean  Ridges 
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THE  ROLE  OF  DENSITY  IN  TOE  ACCUMULATION  OF  BASALTIC  MELTS  AT 
MID-OCEAN  RIDGES 

E.E,  Hooft  and  R.  S.  Detrick 

Department  of  Geology  and  Geophysics,  Woods  Hole  Oceanographic  Institution 


Abstract.  It  is  commonly  assumed  that  magma  ponds  at  a 
level  of  neutral  buoyancy  in  the  shallow  crust  where  melt 
densities  are  equal  to  the  bulk  density  of  the  surrounding 
crust  At  the  East  Pacific  Rise  this  neutral  buoyancy  level  lies 
only  100-400  m  below  the  sea  floor,  significantly  shallower 
than  the  depths  (>  1-2  km)  of  the  magma  bodies  imaged  in 
multichannel  reflection  data,  suggesting  that  other  factors 
must  control  the  collection  of  melt  in  these  reservoirs.  The 
apparent  inverse  reladonship  between  magma  chamber  depth 
and  spreading  rate  at  intermediate  and  fast  spreading  ridges 
suggests  that  the  thermal  structure  of  the  rise  axis,  not  the 
buoyancy  of  melt,  is  the  primary  factor  that  controls  the  depth 
at  which  melt  ponds  in  crustal  magma  chambers  beneath  mid¬ 
ocean  ridges. 

Introduction 

Multichannel  seismic  studies  along  the  fast  spreading  East 
Pacific  Rise  (EPR)  show  the  existence  of  a  narrow  (<l-4  km 
wide),  thin  (10s- 100s  of  m  thick)  sill-like  magma  body  1-2 
km  below  the  sea  floor  (Detrick  et  al.,  1987;  Detrick  ct  al., 
1993).  Despite  its  narrow  width,  reflections  from  the  roof  of 
this  magma  sill  are  remarkably  continuous  along  the  rise  axis 
and  relatively  constant  in  two-way  travel  time  below  the  sea 
floor  over  distances  of  many  tens  of  kilometers.  Along  the 
northern  EPR  between  9®N  and  13°N  this  reflector  typically 
occurs  >  1.2- 1.6  km  below  the  sea  floor  (Detrick  et  al., 
1987).  Along  the  faster  spreading  southern  EPR  between  the 
Garrett  FZ  and  20-7°S  the  magma  sill  is  somewhat  shallower 
with  a  median  depth  of  1 .3  km,  and  in  some  locations  (e.g. 
14°14’S,  17°25’S)  it  rises  to  within  1000  m  of  the  sea  floor 
(Detrick  et  al.,  1993).  Similar  axial  reflectors  have  also  been 
observed  at  intermediate-rate  spreading  ridges  including  the 
Juan  de  Fuca  (Morton  et  al.,  1987)  and  Valu  Fa  (Monon  and 
Sleep,  1985;  Collier  and  Sinha,  1990)  ridges  where  they  are 
significantly  deeper  (>3  km)  than  the  magma  bodies  found 
along  the  EPR  (Figure  1). 

What  controls  the  level  of  magma  ponding  beneath  these 
spreading  centers  and  the  striking  along-axis  uniformity  in  the 
depth  of  magma  sills  apparent  in  multichannel  seismic 
reflection  data?  It  is  commonly  assumed  that  the  ascent  of 
magmas  through  the  crust  beneath  magmatic  systems  is 
primarily  driven  by  the  buoyancy  of  the  melt  relative  to  the 
surrounding  country  rock  (c.f.  Ryan,  1987;  Walker,  1989). 
When  the  density  of  the  melt  and  that  of  the  local  crustal 
rocks  are  equal,  the  melt  is  neutrally  buoyant  and  will  tend  to 
spread  laterally  instead  of  rising  or  sinking.  The  collection  of 
melt  at  a  level  of  neutral  buoyancy  (LNB)  could  explain  the 
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formation  of  the  shallow,  sill-like  crustal  magma  bodies 
observed  along  mid-ocean  ridges.  However,  factors  other 
than  density  can  affect  the  upward  migration  of  melt  within 
the  crust.  For  example,  the  ascending  melt  may  stall  at  a 
rheological  boundary,  like  the  brittle-ductile  transition  (Smith 
and  Cann,  1992)  or  the  freezing  horizon  (Phipps  Morgan 
and  Chen,  1993).  The  mechanics  of  dike  formation  and 
propagation  and  the  state  of  stress  in  the  shallow  crust  may 
also  influence  the  accumulation  of  shallow  magma  reservoirs 
(Rubin,  1990;  Gudmundsson,  1990). 

In  this  paper  we  test  the  hypothesis  that  the  magma  bodies 
observed  beneath  mid-ocean  ridges  form  at  a  LNB  in  the 
shallow  crust  by  comparing  the  seismically  constrained 
depths  of  ridge  crest  magma  chambers  with  the  depth  of  the 
LNB  determined  from  estimates  of  melt  and  crustal  density. 

Estimation  of  melt  and  crustal  densities 

The  LNB  has  been  defined  by  Ryan  (1987)  as  the  depth  in 
the  crust  where  the  magma  density  is  equal  to  the  local  bulk 
density  of  the  surrounding  country  rock.  Below  this  level 
magma  will  be  less  dense  than  its  surroundings  and  tend  to 
rise,  while  above  the  LNB  magma  will  be  negatively  buoyant 
and  will  sink.  Therefore  the  LNB  corresponds  to  the  depth, 
znb.  of  the  cross-over  of  the  effective,  large-scale  in  situ 
density  of  the  crust,  p(z),  with  the  magma  density,  pm*  To 
estimate  the  LNB  at  the  EPR  we  must  know  both  the  magma 
density,  pm.  and  the  variation  of  density  with  depth  in 
oceanic  crust,  p(z). 

Melt  Density:  The  density  of  the  liquid  from  which  a 
basaltic  glass  formed  can  be  computed  using  the  empirical 
method  of  Bottinga  and  Weill  (1970).  They  determined  the 


Fig  1.  Depth  to  the  top  of  the  AMC  versus  spreading  rate 
for  the  northern  (NEPR)  and  southern  (SEPR)  East  Pacific 
Rise,  and  Juan  de  Fuca  (JUFU)  and  Valu  Fau  (LAU)  ridges 
(replotted  from  Phipps  Morgan  and  Chen,  1993).  The  neut¬ 
ral  buoyancy  level  and  critical  depth  calculated  for  the  EPR 
are  also  shown.  The  dashed  curve  represents  the  predicted 
depth  to  the  solidus  (assumed  to  be  1200®C)  as  a  function  of 
spreading  rate  from  Phipps  Morgan  and  Chen  (1993). 
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partial  molar  volume  of  the  major  oxide  components  in  a 
liquid  and  thence  the  density  a  silicate  melt  as  a  function  of 
major  oxide  composition  and  temperaturc. 

Using  these  equations  melt  densities  were  calculated  for 
the  entire  Smithsonian  Institution  Glass  file  consisting  of 
2004  samples  of  basaltic  glasses  from  around  the  wtM'ld’s 
oceans.  A  plot  of  melt  density  versus  Fe#  (Fe/(Fc+Mg), 
-shown  in  Figure  2,  follows  the  theoretical  curves  predicted 
for  the  evolution  of  residual  basaltic  liquids  (Stolper  and 
Walker,  1980).  Most  oceanic  basalts  lie  in  the  density 
minimum  of  this  evolution  curve  suggesting  that  melt  density 
is  an  important  factor  in  determining  the  eruptability  of 
magmas  (Stolper  and  Walker,  1980).  The  average  melt 
density  is  2700  kgm"^;  73%  of  the  samples  lie  within  a  range 
of  ±20  kgm*3  from  this  value. 

Wc  will  adopt  a  value  of  2700  ±  20  kgm'^  for  pm  in  this 
paper.  Density  increases  due  to  continued  crystal  fractiona¬ 
tion  (few  the  range  of  cruptaWe  basalts)  are  only  a  few  percent 
of  this  average  melt  density  and  will  not  significantly  effect 
our  results.  Volatiles  can  lower  melt  densities,  especially  if 
exsolution  of  these  gases  occurs.  But,  the  volatile  content  of 
MORE  is  in  general  very  low  and  the  effect  of  volatiles  on 
density  will  be  small  while  remaining  dissolved  (e.g.  0.25% 
H2O  will  reduce  the  melt  densities  by  only  20  kgm“^).  Thus 
the  effects  of  fractional  crystallization  and  volatiles  on  melt 
density  are  second  OTder  and  largely  offsetting. 

Crustal  Density  Structure:  Seismic  studies  at  the  EPR 
(e.g.  Vera  et  al.,  1990;  Christeson  ct  al.,  1992)  document  the 
presence  of  a  near-surface  layer,  ~150  m  thick  at  the  rise  axis, 
of  extremely  low  (<  3  km/s)  compressional  wave  velocities 
(Figure  3a).  Below  this  layer  p-wavc  velocities  increase 
rapidly  to  ~5  km/s  at  300  m  and  km/s  1-1.2  km  below  the 
sea  floor.  The  low  seismic  velocities  observed  in  the  upper 
kilometer  of  the  crust  are  due  to  high  bulk  porosities  caused 
by  the  presence  of  fractures  and  voids  too  large  to  be  sampled 
in  the  laboratory.  An  accurate  estimate  of  the  crustal  density 
structure  at  a  ridge  crest  thus  requires  a  knowledge  of  the 
variation  in  bulk  porosity  with  depth  in  oceanic  crust. 

The  only  in  situ  measurements  of  bulk  porosity  to 
substantial  depths  within  young  oceanic  emst  have  been  made 
in  Hole  5()4B  on  ~5.9  Ma  old  sea  floor  on  the  flanks  of  the 
Osta  Rica  Rift  in  the  eastern  equatorial  Pacific.  Large-scale 
electrical  resistivities  measured  in  this  borehole  (e.g.  Becker, 
1985)  have  been  used  to  calculate  apparent  bulk  porosities 
using  Archie’s  law  (Figure  3b).  Bulk  porosities  are  high  (10- 


Fig  2.  Plot  of  melt  density  for  MOR  basaltic  glasses 
versus  the  degree  of  fractionation  represented  by  the  Fe#= 
Fc/(Fc+Mg).  The  curve  shows  the  evolution  of  residual  melt 
during  fractional  crystallization  (Stolper  and  Walker,  1980). 
The  mean  melt  density  is  2700  ±  20  kgm“^  (shaded  bar). 
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Fig  3.  (Left)  Seismic  velocity  models  of  the  EPR  crest 
near  9°30‘N  after  CSiristcson  ct  al,  (1992)  (dashed  line)  and 
Vera  ct  al.  (1990)  (solid  line).  (Right)  Variation  of  crustal 
bulk  porosity  with  depth  in  young  oceanic  crust:  based  on  in 
situ  measurements  in  Hole  504B  (from  Becker  et  al.,  1985), 
theoretical  velocity-porosity  relationships  based  on  extended 
Walsh  and  Kuster-Toksoz  theories  (from  Berge  ct  al.,  1992), 
and  an  exponential  decrease  in  bulk  porosity  with  depth  (<^o 
=34.5%;  i  =  3.3*10-3  m-l). 


15%)  in  the  shallowmost  crust,  decrease  to  7-10%  within  the 
lower  part  of  the  extrusive  section  (200-600  m  sub-basement) 
then  drop  sharply  to  values  of  <1-2%  in  the  underlying 
sheeted  dike  section.  Bulk  porosities  remain  extremely  low 
(<1%)  throughout  the  sheeted  dike  section  to  the  base  of 
present  hole  (1726  m  sub-basement). 

In  situ  measurements  of  crustal  physical  properties  have 
not  been  made  at  the  EPR.  However,  an  estimate  of  crustal 
porosity  and  bulk  density  can  be  obtained  by  using  observed 
seismic  velocities  at  the  rise  axis  and  theoretical  predictions  of 
the  variation  in  seismic  velocity  as  a  function  of  pore  volume 
and  aspect  ratio.  Qassical  Hashin-Shtrikman  limits  on  the 
relationship  between  porosity  and  seismic  velocity  ignore 
void  shapes  and  are  too  broad  to  provide  useful  constraints. 
However,  Berge  et  al,  (1992)  have  devised  two  hybrid 
theories,  which  they  call  the  extended  Walsh  and  extended 
Kuster-Toksoz  theories  which  provide  upper  and  lower 
bounds  on  effective  medium  properties  for  porosities  of  up  to 
30%  and  a  distribution  of  pore  aspect  ratios  representative  of 
the  upper  oceanic  crust  (Wilkens  et  al.,  1991).  Figure  3b 
shows  the  variation  in  porosity  with  depth  predicted  using 
these  theories  for  the  (Thristeson  et  al.  (1992)  seismic  model 
of  the  EPR  shown  in  Figure  3a.  Bulk  crustal  porosities  arc 
very  high  (26-43%)  in  the  upper  100  m  then  decrease  to 
values  of  less  than  10%  200  m  into  the  crust.  Note  that 
below  100  m  into  basement  the  results  are  in  close  agreement 
with  in  situ  bulk  porosity  measurements  made  at  Hole  504B. 
Below  about  700  m  depth  large-scale  crustal  porosity 
decreases  exponentially  with  depth  (Figure  3b). 

Porosity  and  density  are  inversely  related  and  the  bulk 
density  at  any  depth  z  below  the  sea  floor  can  be  directly 
calculated  from  the  bulk  porosity,  <^(2),  using  the  expression: 

p(z)  =  pg  -^(z)  (pg  -  pw)  (1) 

where  pg  =  2950  kgm'^  is  the  matrix  grain  density  and  the 
density  of  water,  pw  =  1000  kgm'^.  Figure  4a  shows  p(z) 
for  the  two  EPR  crustal  porosity  models  of  Berge  et  al, 
(1992).  Also  shown  is  p(z)  calculated  assuming  an 
exponential  decrease  in  porosity  with  depth: 
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2000  2200  2400  2600  2800  3000 


Mean  overburden  density,  kg/m^ 

Fig  4.  (Top)  Crustal  density  versus  depth  in  young 
oceanic  crust  based  on  the  bulk  porosity  estimates  shown  in 
Fig.  3.  Also  shown  is  a  density  model  based  on  the  EPR 
seismic  velocity  model  of  Vera  et  al  (1990)  calculated  using 
the  empirical  velocity-density  equation  of  Christensen  and 
Shaw  (1970):  p  =  1.85  +  0.165  Vp.  The  shaded  bar 
represents  the  mean  magma  density  from  Fig.  2.  The  LNB  at 
the  EPR  is  lies  at  100-400  m  depth.  (Bottom)  Variation  in 
mean  overburden  density  versus  depth.  These  curves  were 
calculated  from  the  four  density  models  shown  above.  The 
critical  depth,  zc,  at  the  EPR  is  located  at  400-900  m  depth. 

(t>(z)=<|>oe-^  (2) 

with  an  assumed  surface  porosity,  (j)©,  of  34.5%  and 
A.=3.3*10"^  m‘l.  This  function  is  consistent  with  the  high 
near-surface  crustal  porosities  inferred  seismically  at  the  EPR 
and  provides  a  good  fit  to  the  decrease  in  in  situ  porosities 
measured  in  Hole  504B  below  600  m  depth  (Figure  3b). 

An  independent  estimate  of  the  density  structure  of  oceanic 
crust  at  the  EPR  can  be  obtained  by  converting  measured 
seismic  crustal  velocities  to  density  using  empirical 
relationships  determined  from  laboratory  measurements  of 
compressional  wave  velocity  and  density  made  on  oceanic 
dredge  samples  (Christensen  and  Shaw,  1970).  Because 
laboratory  velocity  measurements  don’t  take  into  account 
large-scale  crustal  porosity,  this  approach  will  be  most 
accurate  for  the  deeper  parts  of  the  crustal  section.  The 
seismic  structure  of  the  rise  axis  was  taken  from  an 
expanding  spread  profile  shot  along  the  crest  of  the  EPR  at 
9°N  (V era  et  al.,  1990).  The  resulting  p(z)  curve  is  plotted  in 
Figure  4a.  Note  that  the  decrease  in  compressional  wave 
velocity  between  1.2  and  1.6  km  depth  in  the  Vera  et  al. 
(1990)  model  marks  the  top  of  the  axial  magma  chamber;  the 
inferred  decrease  in  crustal  density  below  this  depth  is  due  to 
the  presence  of  magma  and  is  not  an  intrinsic  feature  of  the 
density  structure  of  oceanic  crust 

The  three  different  methods  of  estimating  the  density 
structure  of  oceanic  crust  at  the  EPR  described  above  all 
provide  remarkably  consistent  results  (Figure  4a).  Bulk 


crustal  densities  are  expected  to  be  very  low  (<2500  kgm"^) 
in  the  upper  100  m  of  the  crust.  They  increase  rapidly  to 
about  27(X)  kgm"^  -150-400  m  into  the  crust,  and  reach 
values  of  2800-2900  kgm"^  at  1000  m  depth. 

Discussion 

A  LNB  will  exist  in  the  shallow  oceanic  crust  where  the 
average  melt  density  (2700  +  20  kgm"^)  intersects  the  crustal 
density  curve.  For  the  data  shown  in  Figure  4a  the  LNB  is 
located  about  100-4(X)  m  below  the  sea  floor  at  the  EPR.  The 
LNB  cannot  lie  significantly  deeper  than  this  at  the  EPR 
without  being  in  serious  disagreement  with  the  known 
seismic  velocity  structure  of  the  EPR,  theoretical  and 
empirical  relationships  between  seismic  velocity,  bulk 
porosity  and  density,  and  in  situ  measurements  of  the 
physical  properties  of  young  oceanic  crust 

If  density  were  the  principal  factor  controlling  the 
accumulation  of  melt  beneath  mid-ocean  ridges,  as  some 
investigators  have  suggested  (e.g.  Ryan,  1987),  we  would 
expect  crustal  magma  bodies  to  form  at  the  depth  of  this 
LNB,  i.e.  about  100-400  m  below  the  sea  floor.  However, 
the  axial  magma  chambers  documented  seismically  along  the 
EPR  occur  1-2  km  below  the  sea  floor,  well  below  the 
calculated  depth  of  the  LNB  (Figure  1).  Even  along  the  fast 
spreading  southern  EPR,  where  the  shallowest  crustal  magma 
bodies  have  been  documented,  the  top  of  the  magma  sill  is 
typically  more  than  1000  m  below  the  sea  floor  (Detrick  et 
al.,  1993).  Along  the  Juan  de  Fuca  (Morton  et  al.,  1987)  and 
Valu  Fa  (Morton  and  Sleep,  1985;  Collier  and  Sinha,  1990) 
ridges  magma  chambers  lie  deeper  than  3  km.  We  thus 
conclude  that  magma  chambers  beneath  intermediate  and  fast 
spreading  ridges  do  not  form  at  a  LNB  and  other  factors  must 
control  the  accumulation  of  melt  beneath  mid-ocean  ridges. 

All  of  the  mid-ocean  ridge  magma  chambers  observed  to 
date  lie  at  or  below  what  we  call  the  critical  depth,  zc,  the 
depth  below  which  lithostatic  pressure  alone  provides  enough 
hydrostatic  head  to  force  magma  from  the  reservoir  to  the 
surface.  At  the  critical  depth  pav(zc)  =  Pm  where  pav(z)  is 
the  mean  overburden  density  at  depth  z. 

Figure  4b  shows  the  mean  overburden  density  versus 
depth,  pav(z)»  calculated  for  the  four  crustal  density  models 
shown  in  Figure  4a.  Our  best  estimate  of  the  critical  depth  at 
the  EPR  is  that  it  lies  400-900  m  below  the  sea  floor.  The 
axial  magma  chambers  imaged  seismically  at  the  EPR  occur 
close  to,  but  below,  the  critical  depth  while  those  found  along 
the  intermediate  spreading  Valu  Fa  and  Juan  de  Fuca  Ridges 
lie  well  below  the  critical  depth  (Figure  2).  Magma  that 
ponds  above  the  critical  depth  may  still  be  cniptable  since 
other  mechanisms  may  increase  magma  reservoir  pressures 
above  this  minimum  value,  but  for  magma  bodies  lying 
below  the  critical  depth  lithostatic  pressure  alone  is  sufficient 
to  insure  eruptability  provided  a  continuous  pressure 
connection  (e.g.  a  dike)  exists  with  the  surface. 

If  magma  ponds  at  mid-ocean  ridges  in  crustal  reservoirs 
well  below  the  LNB,  as  the  evidence  presented  in  this  paper 
strongly  suggests,  what  stops  the  buoyancy-driven  ascent  of 
melt  through  the  crust?  The  answer  to  this  question  is  not 
entirely  clear.  Gudmundsson  (1990)  has  proposed  that  rapid 
injections  of  dikes  at  a  fast  spreading  ridge  create  a  stress 
barrier  in  the  mid-crust  where  the  horizontal  compressive 
stress  is  higher  than  the  vertical  stress  causing  sill  formation 
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that  leads  to  the  development  of  a  mid-crustal  magma 
chamber.  Rubin  (1990)  showed  that  in  regions  undergoing 
active  extension,  the  brittle-ductile  boundary  may  become  a 
depth  of  ^effective**  neutral  buoyancy  due  to  the  accumulation 
of  non-hydrostatic,  tectonically-induced  stress  differences. 
Phipps  Morgan  and  Chen  (1993)  have  proposed  that  melt  will 
be  stopped  in  its  upward  ascent  when  it  reaches  the  freezing 
horizon  due  to  viscous  stresses  that  favor  magma  ponding. 

Whether  it  is  the  brittle-ductile  boundary  or  the  freezing 
horizon  that  effectively  stops  the  upward  ascent  of  melt,  the 
depth  of  both  will  be  primarily  controlled  by  the  thermal 
structure  of  the  crust  at  Ae  ridge  axis.  The  strongest  evidence 
for  a  thermal  control  on  the  depth  at  which  magma 
accumulates  beneath  mid-ocean  ridges  is  the  variation  in 
depth  of  crustal  magma  chambers  with  spreading  rate,  first 
noted  by  Purdy  ct  al.  (1992),  and  illustrated  in  Figure  1.  The 
hypothesis  that  the  depth  magma  bodies  beneath  mid-ocean 
ridges  are  primarily  controlled  by  thermal  processes  has 
recently  been  tested  quantitatively  by  Phipps  Morgan  and 
C!hen  (1993).  They  show  that  the  depth  of  the  1200®C 
isotherm  varies  with  spreading  rate  in  a  way  that  is  consistent 
with  the  variation  in  the  depth  to  axial  magma  bodies 
determined  seismically  along  intermediate  and  fast  spreading 
ridges  (Figure  1).  Their  results,  and  the  seismic  evidence  that 
axial  magma  chambers  lie  below  the  LNB,  suppon  the  hypo¬ 
thesis  that  the  buoyancy  of  melt  is  not  the  primary  factor 
controlling  the  depth  of  magma  bodies  beneath  ridge  crests. 

Conclusions 

1.  The  crustal  magma  bodies  observed  beneath  mid-ocean 
ridges  do  BSil  form  at  a  level  of  neutral  buoyancy  in  the 
shallow  crust.  The  neutral  buoyancy  level  lies  only  100- 
400m  below  the  sea  floor,  much  shallower  than  the  depth  (1- 
3  km)  of  the  axial  magma  bodies  imaged  in  multichannel 
reflection  data  at  intermediate  and  fast  spreading  ridges. 

2.  The  magma  bodies  found  along  mid-ocean  ridges  are 
located  at  or  below  the  critical  depth  where  lithostatic  pressure 
due  to  the  integrated  weight  of  the  overlying  crustal  column 
alone  is  sufficient  for  magma  to  be  erupted  onto  the  sea  floor. 

3.  The  apparent  inverse  relationship  between  magma 
chamber  depth  and  spreading  rate  at  intermediate  and  fast 
spreading  ridges  suggests  that  the  thermal  structure  of  the  rise 
axis  is  the  primary  factor  that  controls  the  depth  at  which  melt 
ponds  in  crustal  magma  chambers  beneath  mid-ocean  ridges. 
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Abstract 


A  stochastic  m^l  for  the  emplacement  of  dikes  and  lava  flows  at  a  fast  spreading  ridge  can  generate  an  upper  oceanic 
crustal  structure  sirndar  to  that  observed  in  seismic  data  from  the  East  Pacific  Rise  (EPR),  in  ocean  drill  holes,  and  in 
ophiolites.  In  this  model  the  location  of  successive  dike  intrusion  events  relative  to  the  ridge  axis  is  determin^  by  a 
Gausstan  probabihty  function, and  the  cumulative  flow  lengths  of  the  erupted  lavas  ate  chosen  to  build  a  Gaussian-shaped 
lava  pile.  We  interpret  wide-angle  seismic  reflections  from  the  steep  velocity  gradient  at  the  base  of  seismic  layer  2A  to  be 
the  extrusrve/sheeted  dike  contact.  Seismic  data  from  the  northern  and  southern  EPR  place  constraints  on  the  on-axis 
extrusive  layer  thickness  (230  ±50  m),  the  distance  over  which  the  thickening  of  the  extrusive  layer  occurs  (width  of  the 
accretion  zone=  1-3  km)  and  its  off-axis  thickness  (300-800  m).  Ophiolites  and  ocean  drill  holes  (DSDP  Hole  504B) 
provKfc  additional  estirnates  of  the  thickness  of  the  extrusive  Uyer  and  constrain  the  thickness  of  the  transition  region  from 
extrusives  to  sheeted  dikes  (~  100-200  m).  A  simple  stochastic  emplacement  model,  where  the  lavas  ate  described  by  one 
mean  flow  length,  fits  the  thickening  of  the  extrusive  layer  off-axis  inferred  from  the  deepening  of  seismic  layer  2A.  but  the 
predicted  transition  from  sheeted  dikes  to  extrusives  is  too  thick.  In  order  to  match  the  dimensions  and  flat-topped  shape  of 
the  seismic  layer  2A  boundary  as  well  as  the  thickness  of  the  extrusive-sheeted  dike  transition,  a  bimodal  distribution  of 
Uva  flows  IS  used.  Short  flows,  confined  within  the  axial  summit  caldera  (ASC),  build  up  approximately  half  the  extrusive 
volume.  Occasional  voluminous  flows  spill  out  of  the  ASC.  or  erupt  outside  of  the  ASC.  and  pond  at  a  considerable  distance 
off-ai^  to  build  up  the  remamder  of  the  extrusive  section.  The  upper  part  of  the  final  extrusive  section  will  be  dominated  by 
the  off-axis  flows,  while  the  lower  portions  wiU  be  primarily  composed  of  short  flows  erupted  within  the  ASC  Magnetic 
transition  widths  predicted  from  the  overlap  of  lavas  (~  2  km)  in  this  model  are  similar  to  those  measured  in  deep-tow 
studies.  Assuming  a  smoothing  function  which  acts  over  one  seismic  wavelength,  the  upper  crustal  velocity  structure 
^ic^  by  the  bimodal  lava  emplacement  model  is  consistent  with  the  shallow  seismic  velocity  structure  measured  on  the 
EPR.  The  ages  of  seafloor  lavas  in  this  model  are  younger  than  die  tectonic  spreading  model  ages  by  ~  30-70  kyr,  in 
agreement  with  anomalously  young  lava  ages  obtained  from  radioisotope  dating  of  seafloor  basalts  near  the  EPR. 
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1.  Introduction 

Tlie  upper  portion  of  the  oceanic  cnist  is  formed 
by  dike  intrusion  and  the  eruption  of  lava  flows. 
These  two  magmatic  processes  generate  an  upper 
crustal  structure  composed  of  lavas  piled  on  top  of 
one  another,  underlain  by  a  sheeted  dike  complex. 
Constraints  on  the  processes  of  upper  crustal  con¬ 
struction  have  come  from  ofAiolite  studies  [1,2], 
observations  in  active  rifts  (e.g.,  Iceland  [3D,  model¬ 
ing  of  marine  magnetic  anomaly  shapes  and  transi¬ 
tion  thicknesses  [4-8],  mapping  [9-14],  dating  [15], 
and  drilling  (e.g.  DSDP  Hole  504B  [16D  of  young 
ocean  crust  Oj^iolite  studies,  on  the  basis  of  chilled 
margin  statistics  [1,17,18]  and  measurements  of  the 
thickness  of  the  lava/dike  transition  [1]  predict  nar¬ 
row  intrusion  zone  widths  of  10-50  m.  Modeling  of 
marine  magnetic  anomaly  transition  widths  suggests 
a  lava  emplacement  zone  1-3  km  wide.  Geological 
mapping  and  dating  of  lava  flows  at  the  East  Pacific 
Rise  (EPR)  suggest  that  volcanic  activity  is  primarily 
confixed  within  the  narrow  ( 50-200  m  wide) 
axial  summit  caldera  (ASC),  but  there  is  evidence  for 
very  young  flows  on  the  flanks  of  the  axial  high 
outside  the  ASC  that  suggest  eruptive  activity  over  a 
much  wider  zone  or  transport  of  lavas  a  considerable 
distance  off-axis  (e.g.  [19D.  The  relative  importance 
of  off-axis  eruptive  activity  and  the  processes  that 
control  the  thickness  of  the  extrusive  section  are  still 
poorly  understood. 

Advances  in  seismic  imaging  of  the  shallow  ve¬ 
locity  structure  at  the  East  Pacific  Rise  allow  the 
thickness  of  what  is  commonly  interpreted  as  the 
extrusive  layer  to  be  directly  mapped.  In  this  paper, 
we  use  the  variation  in  this  seismically  determined 
extrusive  thickness  to  provide  new  constraints  on 
crustal  emplacement  processes.  We  show  that  the 
dike  intrusion  zone  must  be  extremely  localized  about 
the  ridge  axis.  The  extrusive  section  is  made  up  of 
two  units  of  about  ^ual  thickness:  a  lower  unit, 
composed  of  short  flows  confined  within  the  axial 
summit  caldera,  and  an  upper  unit,  of  flows  which 
spill  oiit  of  the  axial  summit  caldera  (ASC),  or  are 
channelled  across  the  shoul<fcrs  of  the  rise  axis  in 
lava  tubes  and  flow  long  distances  down  the  fljmks 
of  the  axial  high,  or  erupt  off-axis.  We  discuss  the 
implications  of  this  emplacement  model  on  the  seis¬ 
mic  velocity  structure  of  the  shallow  oceanic  crust. 


subsidence  rates,  dip  of  isochrons  in  the  extrusive 
section,  magnetic  anomaly  transition  widths,  and  the 
age  of  lavas  exposed  at  the  seafloor. 

2.  Observational  constraints  on  upper  oceanic 
crustal  structure  at  the  EPR 

The  typical  seismic  velocity  structure  of  the  shal¬ 
low  oceanic  crust  has  been  characterized  from  high 
quality  expanding  spread  profiles  (ESP)  [20,21],  to¬ 
mographic  images  [22],  refraction  data  [23],  wide 
aperture  profiles  (WAP)  [24,25]  and  on-bottom  re¬ 
fraction  work  (NOBEL)  [26].  These  studies  have 
identified  a  small  velocity  layer  which  is  100-200  m 
thick,  has  velocities  ^  2.6  km/s  and  a  low  vertical 
velocity  gradient  This  surficial  low  velocity  layer  is 
underlain  by  a  steep  gradient  zone  where  velocities 

>  ^  5  km/s  are  reached  in  a  few  hundred  meters. 
The  crustal  section  with  the  surficial  low  velocity 
layer  and  the  underlying  steep  velocity  gradient  is 
called  layer  2A.  The  underlying  crust,  with  velocities 

>  ^  5  km/s  and  a  very  gradual  velocity  gradient, 
corresponds  to  layer  2B  [27].  The  base  of  the  high 
gradient  region  is  imaged  in  multichannel  seismic 
reflection  (MCS)  experiments  as  a  shallow  wide 
angle  reflection  formed  by  refracted  energy  turning 
in  the  steep  gradient  [24,28,29].  No  normal  incidence 
reflection  is  seen,  since  there  is  no  impedance  con¬ 
trast,  rather  the  steep  gradient  generates  diving  rays 
which  are  observed  at  large  offsets  on  the  streamer. 
In  Fig.  lb,d,  the  wide  angle  reflection  from  the  base 
of  layer  2A  is  outlined  with  a  gray  line  for  two 
seismic  profiles  shot  across  the  rise  crest  of  the 
southern  EPR  (SEPR31  at  iri4'S,  and  SEPR06  at 
18®28'S).  To  obtain  the  thickness  of  layer  2 A  the 
two-way  travel  times  to  the  base  of  layer  2A  are 
converted  to  depth  using  the  composite  velocity 
model  of  Harding  et  al.  [28]. 

The  shallow  velocity  structure  of  the  ocean  crust 
is  controlled  by  the  composition  and  porosity  of  the 
upper  crust.  We  interpret  layer  2 A  in  young  (<  1 
My)  oceanic  emst  to  be  die  extrusive  section  and  use 
measurements  of  its  thickness  to  map  the  variability 
in  the  accumulaticwi  of  extrusives.  The  surficial  low- 
velocity  zone  is  attributed  to  a  layer  of  high-porosity 
pillows  and  flows  (e.g.  [30,3  ID,  with  porosities  pos¬ 
sibly  greater  than  20%  [32-34].  The  higher-velocity 
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a) 


3,50 


4,00 


b) 


3,50 


4,00 


SEPR  6,  18*‘28’$,  finite  difference  migration 


3,50 


4,00 


d) 


3.50 


4.00 


Distance  from  ridge  axis,  km 

Rg.  1.  (a)  and  (c)  Across-axis  multichannel  seismic  sections  from  the  southern  East  Pacific  Rise,  SEPR31  (IT’14"S)  and  SEPR06  (I8®28"S) 
[40].  The  data  were  stacked  using  a  velocity^Iepth  function  that  varied  laterally  and  then  finite  difference  migrated  to  account  for  the 
seafloor  and  layer  2A  topograjAy.  The  axial  magma  diambcr  (AMC)  and  the  base  of  layer  2A  arc  indicated  on  the  figures,  (b)  and  (d)  Line 
drawings  showing  picks  of  the  seafloor  and  layer  2A  reflections  in  black  and  dark  gray  lines,  respectively.  The  light  gray  constant  thickness 
band  corresponds  to  the  on-axis  thickness  of  layer  2A.  This  emphasizes  that  the  on-axis  thickness  is  built  up  very  rapidly,  while  the  off-axis 
accretion  of  layer  2A  occurs  over  a  few  kilometers.  Both  profiles  show  similar  on-axis  and  off-axis  layer  2A  thicknesses,  -  250  m  and 
500-550  m,  respectively.  The  two  profiles  represent  end-members  of  the  width  of  the  region  over  which  the  accretion  of  layer  2A  takes 
place:  on  SEPR31  it  is  narrow  while  on  SEPR06  it  is  broad  and  the  widths  of  tfie  2A  accretion  zone  arc  1.2  km  and  3  km.  respectively.  Note 
tiiat  tiic  layer  2A/2B  boundary  has  a  flat-topped,  concave  shape  beneath  the  rise  axis. 
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Table  1 


Seismic  constraints  on  the  variability  of  layer  2A  along  tfie  EPR 


Experiment,  latitude 

On-axis  2A  thickness 

Off-axis  2A  thickness 

Thickness  of  2A 

accictioo  zone 

MCS  NEPR,  8*50'-.9®50'N  * 

~  230  ±  50  m 
—  400  m  near  9®N  OSC 

440-560  m 

0.8-2  km 

Mcs  SEPR.  ir-a)^s  * 

235  ±  40  m 

300-400  m  near  OSCs 

400-800  m 

0.5-3  km 

WAP.WN' 

~  150  ±  30  m 

250-280  m 

3km 

ESP.yayN  ** 

150  m 

310-350 m 

NOBEL,  9^'N  ' 

130  m 

~300m 

*  Harding  ct  al.  [28).  ^  Hooft  ct  al.  [41].  ^  Vera  EHcbold  [24].  ^  Vera  ct  al.  (20).  ^  Christeson  ct  al.  [36]. 


layer  2B  material  underlying  the  surficial  low- veloc¬ 
ity  layer  has  been  interpreted  to  be  either  the  top  of 
the  sheeted  dike  sequence  [28,35,36]  or  low-porosity 
pillow  lavas  and  flows  [20,21,23,37].  At  DSDP  Hole 
504B  —  6  My  old,  intermediate  spreading  crust  — 
high  seismic  velocities  arc  reached  within  the  extru¬ 
sive  section  [38].  At  Hess  Deep  —  young,  fast 
spreading  crust  —  the  transition  from  extrusives  to 
sheeted  dikes  occurs  at  a  depth  corresponding  to  the 
location  of  the  base  of  layer  2A  in  velocity-depth 
profiles  measured  on  fast-spreading  crust  [26].  We 
follow  the  reasoning  of  Christeson  et  al.,  1994  [27] 
in  considering  the  velocity  structure  measured  at 
DSDP  Hole  504B  to  be  the  result  of  the  alteration 
and  reduction  of  bulk  porosity  of  the  shallow  crust 


With  increasing  age  the  base  of  layer  2A  migrates 
from  die  sheeted  dike/extrusive  boundary  upwards 
to  a  porosity  contrast  within  the  extrusive  section.  By 
comparison  with  Hess  Deep,  we  interpret  layer  2A  to 
consist  of  both  the  extrusive  and  mixed  zones  (i.e.,  it 
is  a  transition  zone  where  lavas  and  sheeted  dikes 
intermingle).  The  high  velocity,  low  gradient  mate¬ 
rial  represents  the  top  of  the  sheeted  dike  complex 
(layer  2B).  The  wide  angle  multichannel  reflections 
are  due  to  rays  turning  in  the  steep-velocity  gradient 
in  the  transition  zone  at  the  base  of  layer  2A. 

Large  lateral  variations  in  layer  2A  thickness  are 
observed  in  multichannel  seismic  data  at  the  fast¬ 
spreading  northern  [24,28]  and  southern  [29,39-41] 
EPR,  in  wide  aperture  profiles  [24]  in  the  9°N  area. 


Tabic  2 


Geological  observations 


Geology 

Observation 

Type 

Location 

Spreading  rate 

Extrusive 

thickness 

Extrusive -dike 
transition 

duckness 

Ophiolites 

Betts  Cove,  Newfoundland 

slow 

1000  m 

-25m* 

Bay  of  Islands,  Newfoundland 

slow 

~700m*’ 

~  50  m  * 

Troodos,  Cyprus 

intermediate 

llOOm^ 

100-200 

Semail,  Oman 

fast 

400-1600 m  ^ 

few  lens  of  meters- 100  m 

Oceanic  sections 

504B  data 

intermediate 

575  ra® 

150-200 

Hess  I>ecp 

fast 

Alvm:  <  200-800  m  ^ 
Nautile:  100-200  m  * 

50-500  m  * 

Vema  Scarp 

slow 

280-420  m  “ 

60-300  m  ** 

Blanco  Scarp 

intermediate 

475  m  ‘ 

785  m  ‘ 

*  Kidd  [47].  ^  Kidd  [1].  ^  Uppaid  ct  ai.  [48].  ^  Nicolas  [42],  p.  61.  ®  Anderson  [16].  ^  Karson  ct  al.  [43]  and  pers.  commun.,  1994. 

*  Franchctcaui  ct  al.  [44].  **  Camiat  ct  al.  [45].  *  Naidoo  ct  al.  [34]. 
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in  expanding  spread  profiles  [20]  at  9®35'N,  and  in 
on-bottom  refraction  studies  [27]  at  9®30'N.  Table  1 
shows  die  on-axis  and  off-axis  thicknesses  and  the 
width  of  the  2A  accretion  zone  —  the  distance  from 
the  axis  over  which  layer  2A  thickens  to  die  mean 
off-axis  thickness  —  for  different  seismic  experi¬ 
ments  on  the  northern  and  southern  EPR,  Moving  off 
the  rise  crest,  layer  2A  is  observed  to  approximately 
double  in  thickness  over  a  distance  of  a  few  kilome¬ 
ters.  In  general,  on-axis  layer  2A  thicknesses  are 
remarkably  uniform,  230  ±50  m  except  near 
overlapping  spreading  centers  where  the  thicknesses 
reach  400  m  [283,41].  Off-axis  thicknesses  of 
400-800  m  are  reached  within  1-3  km  of  the  rise 
axis.  In  spite  of  the  differences  in  spreading  rate, 
layer  2A  thicknesses  are  similar  along  the  northern 
(55  mm/yr  half-rate)  and  southern  (75  mm/yr  half¬ 
rate)  EPR,  although  different  seismic  techniques  give 
somewhat  different  results,  due  to  the  differences  in 
inherent  resolution  [26],  If  the  velocities  in  the  extru- 
sives  and  dikes  are  assumed  to  be  constant  and  the 
steep  velocity  gradient  at  the  base  of  layer  2A  is  due 
to  the  interfrngering  of  these  two  units,  then  the 
thickness  of  the  high  velocity  gradient  zone  corre¬ 
sponds  to  the  thickness  of  the  transition  zone  from 
extrusives  to  sheeted  dikes.  In  seismic  data  from  die 
northern  EPR  this  high  velocity  gradient  zone  has  a 
thickness  of  ~  100-300  m  [20,26]. 

These  seismically  derived  constraints  on  extrusive 
layer  thickness  and  thickness  of  the  extrusive-dike 
transition  can  be  compared  with  geological  observa¬ 
tions  on  sections  of  oceanic  crust  exposed  in  Hess 
Deep,  in  ophiolites,  and  drilled  in  Hole  504B  (Table 
2).  Extrusive  layer  thicknesses  estimated  in  ophio¬ 
lites  range  from  400  to  1600  m  [1,42].  In  contrast, 
estimates  of  extrusive  thicknesses  from  oceanic 
scarps  are  only  300-500  m  [43-46]  while  at 
DSDP  Hole  504B  the  extrusive  section  is  less  than 
600  m  thick  [16],  comparable  to  layer  2A  thicknesses 
determined  from  off-axis  seismic  data  (Table  1).  The 
extrusive-dike  transition  zone  (defined  as  the  verti¬ 
cal  depth  from  20%  to  70%  sheeted  dikes)  is  quite 
thin  in  ophiolite  sections  (Table  2),  a  few  tens  of 
meters  to  100-200  m  [1,47,48].  Oceanic  scarps  give 
values  also  ranging  from  a  few  tens  of  meters  up  to  a 
few  hundred  meters  [44-46],  while  at  DSDP  Hole 
504B  die  extrusive-dike  transition  is  --  150  m  thick 
[16].  These  estimates  are  comparable  to  the  thickness 


of  the  layer  2A/2B  transition  estimated  seismically. 
Modeling  of  the  extrusive-dike  transition  thickness 
in  ophiolites  shows  that  the  dike  intrusion  zone  must 
be  very  narrow,  only  10-50  m  wide  [1].  Such  a 
narrow  dike  intrusion  width  is  supported  by  indepen¬ 
dent  statistical  observations  of  the  proportion  of  dikes 
in  ophiolites  that  exhibit  one-sided  cooling,  due  to 
having  been  split  by  a  subsequent  dike  [1,17,18]. 
Studies  of  marine  magnetic  anomaly  reversal  widths 
and  variability  [4-8]  suggest  a  relatively  wide  zone 
of  lava  emplacement,  1-3  km  wide.  However,  these 
data  can  not  resolve  die  relative  contributions  of  the 
width  of  the  intrusion  zone  and  lava  flow  distribution 
to  the  lava  emplacement  process.  Macdonald  [6] 
attributes  the  width  of  the  magnetic  anomaly  reversal 
to  a  combination  of  a  relatively  wide  intrusion  zone 
and  long  lava  flows. 

There  are  four  observables  which  help  constrain 
the  dimensions  of  the  upper  crustal  emplacement 
process:  the  on-axis  2A  thickness,  off-axis  2A  diick- 
ness,  width  of  the  2A  accretion  zone,  and  the  extru¬ 
sive-sheeted  dike  transition  thickness.  These  values 
are  obtained  from  seismic  and  geological  observa¬ 
tions.  In  this  paper  we  attempt  to  match  two  end- 
member  profiles  of  layer  2A  thickening  from  the 
southern  EPR.  At  the  same  time  we  maintain  an 
extrusive-dike  transition  thickness  that  is  consistent 
with  geological  observations  —  greater  than  50  m 
and  no  wider  than  150-200  m.  The  two  across-axis 
seismic  sections  from  the  Southern  East  Pacific  Rise, 
SEPR31,  irM'S,  and  SEPR06,  18^28'S,  (Fig.  la,c) 
were  chosen  because  they  characterize  the  variability 
in  the  width  and  shape  of  thickening  of  layer  2A. 
Although  both  of  these  profiles  (Fig.  lb,d)  have 
similar  on-axis  (250  m)  and  off-axis  layer  2A  thick¬ 
nesses  (550  m),  they  differ  in  the  width  of  the  zone 
over  which  the  diickening  takes  place;  SEPR31  is 
narrow  (width  of  tiie  2 A  accretion  zone  is  1.2  km) 
while  SEPR06  thickens  over  a  wide  zone  (width  of 
the  2A  accretion  zone  is  3  km).  Both  diickening 
profiles  are  flat-topped  (i.e.,  the  minimum  on-axis 
thickness  is  maintained  for  some  distance  before  the 
layer  starts  to  thicken). 
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3.  Emplacement  models 

3.L  Unimodal  lava  emplacement  model 

We  use  a  stochastic  model  to  construct  the  sheeted 
dike  and  extrusive  sections  of  the  oceanic  crust  (c.f. 
[4^]).  At  an  active  spreading  center,  steady  state  is 
only  achieved  statistically,  since  the  accretion  pro¬ 
cess  is  a  series  of  discrete,  successive  events.  Fig.  2 
shows  a  schematic  example  of  four  steps  in  the 
emplacement  process.  Dikes  are  injected  into  the 
crust  at  a  location  determined  randomly  jffom  a 
probability  function.  Each  dike  is  associated  with  a 
lava  flow.  Lava  flows  out  on  either  side  of  the  locus 
of  intrusion  to  a  distance  determined  randomly  from 
a  second  probability  function.  Earlier  flows  are  cov¬ 
ered  by  later  flows.  The  lava  pile  subsides  as  it  is 
covered  so  that  a  flat  seafloor  is  maintained.  The 
variables  used  in  the  emplacement  models  are  de¬ 
fined  in  Table  3.  The  dikes  are  assumed  to  have  a 
thickness,  w  —  1  m.  This  is  consistent  with  dike 
widths  of  0.5- 1.5  m  observed  at  Hess  Deep  [43]  and 
a  mean  of  1  m  for  the  dike  width  distribution 
measured  in  the  Troodos  ophiolite,  Cyprus  [1].  The 
source  region  of  the  dikes  lies  below  the  base  of  our 
model.  The  probability  function  of  the  location  of 
dike  intrusion  relative  to  the  ridge  axis  that  we  use  is 
a  Gaussian,  with  a  standard  deviation  of  intrusion, 
a^.  Of  the  dikes,  95%  are  intruded  within  a  distance 
of  ±2  a,  from  the  ridge  axis.  Lava  flow  lengdis  are 
such  that  cumulative  flows  erupted  at  the  same  loca¬ 
tion  build  a  Gaussian-shaped  extrusive  pile  with 


Time  1 


Time  2 


Time  4 


Fig.  2.  Cartoon  of  four  steps  in  the  stochastic  emplacement 
process  that  constructs  the  upper  crusL  The  location  of  dike 
intrusion  is  detOTnincd  from  die  probability  function  for  tiK 
intrusion  locaticm  relative  to  the  ridge  axis.  Pre-existing  crust 
spreads  to  the  side  a  distance  equal  to  half  a  dike  width.  The 
length  of  the  lava  flow  erupted  from  this  dike  is  determined  from 
the  probability  frmctioo  describing  the  flow  lengths.  The  lava  flow 
extends  equal  distances  on  either  side  of  the  point  of  intrusioo. 
Pre-cxisting  crust  underlying  the  new  flow  subsides  an  amount 
equal  to  one  flow  diickness.  The  cycle  then  repeats  itself  with  die 
selection  of  a  new  point  of  intrusion  and  the  determination  of  the 
accompanying  lava  flow  length. 


Table  3 

Model  variables 

Model  variable  Definition 


Unimodal  Lava  Emptacement  Model 

psO 

t 

w 


standard  deviatiem  of  dike  intrusion 
standard  deviation  of  lava  emplacement 
probability  diat  a  dike  erupts  to  form  a  lava  flow 
diickness  of  the  lava  flows 
width  of  die  dikes  =  1  m 


Bintodal  Lava  Em^acement  Model 

mini 

t 

psO 

w 


Standard  deviation  of  dike  intrusion  for  all  flows  —  10  m 

standard  deviation  of  lava  emplacement  of  the  short  flows  (95%  of  events) 

standard  deviati<^  of  lava  cmplacem^t  of  die  long  flows  (5%  of  events) 

distance  over  which  no  lava  accumulates  for  die  long  flows 

diickaess  of  the  lava  flows  (bodi  short  and  long) 

probability  that  a  dike  erupts  to  form  a  lava  flow  =  1 

width  of  the  dikes  *=  I  m 
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Standard  deviation  a^.  Thus,  the  half-flow  length 
probability  function  that  we  chose  is  the  derivative 
of  a  half-Gaussian,  which  is  the  Rayleigh  distribu¬ 
tion.  The  flows  have  a  uniform  thickness,  r,  and  die 


median  flow  length  is  (2Tr)‘/^cr^.  Observations  in 
Hawaii  and  Iceland  suggest  that  a  significant  propor¬ 
tion  (50-75%)  of  diking  events  do  not  form  a  sur¬ 
face  eruption  (Delaney,  pers.  commun.,  RIDGE  The- 


Unimodal  Lava  Emplacement  Model 

Isochrons:  ci  =  300  m 


“loool 


“^“^-2  -1  0  1  2345 


800  • - J - 1 - 1 - 1 - 1 _ 1 _ I _ I _ 1 _ I 

-5  -4  -3  -2  -1  0  1  2  3  4  5 

Distance,  km 

Fig.  3.  Two  nms  of  Ac  unimodal  lava  emplacement  model  of  upper  crustal  constnictioa  where  Ac  standard  deviation  of  intrusion  is  wide, 
a,  =  300  m  (a  and  b)  and  narrow,  o,  =  10  m  (c  and  d),  Ac  standard  deviation  of  cxnusion,  o,  =  500  m,  Ac  thickness  of  Ae  dikes  and  of  Ac 
lava  flows  is  1  m  and  40%  of  Ac  dikes  erupt  a  lava  flow,  (a)  and  (c)  The  isochrons  wiAin  Ac  extrusives  are  shown  by  changing  shades  of 
gray.  The  Akes  arc  shown  in  white.  No  vertical  exaggeration,  (b)  and  (d)  The  gray  line  is  Ac  boundary  between  Ac  sheeted  dikes  and 
extrusives.  The  hollow  Imc  is  Ac  result  of  spatially  Altering  Ae  extrusive-dike  boundary  wiA  a  triangular  Alter  of  wavcIengA  200  m  to 
simulate  a  seismic  observation  of  Ais  contact  The  solid  black  line  is  Ae  across-axis  layer  2A  thickness  Aom  line  SEPR3I.  For  Ac  model 
wiA  a,  =  300  m  Ae  character  as  well  as  Ac  on-axis  (200  m).  off-axis  thicknesses  (510  m),  and  widA  of  Ac  2A  accretion  zone  (2300  m) 
mat^  Ac  seismic  proflle  well.  However,  Ac  transition  zone  Aom  extrusives  to  sheeted  Akes  is  310  m,  thicker  than  fliat  observed  in 
ophiolite  sections  and  oceanic  scarps.  When  o,  =  10  m.  ^  transition  zone  Aom  extrusives  to  sheeted  dikes  is  50  m,  comparable  to  ophiolitc 
and  oceanic  scarp  observations.  However,  Ac  shape  of  the  layer  2A  boundary  is  pointed  rather  than  flat-topped,  wiA  a  very  small  on-axis 
thickness  (70  m)  and  does  not  resemble  Ae  seismic  profllc. 
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oretical  Institute,  June,  1995),  This  results  in  addi¬ 
tional  extension  during  lava  emplacement  or,  alterna¬ 
tively,  represents  a  change  in  the  effective  ratio  of 
the  thickness  of  the  dikes  to  that  of  the  lava  flows. 
We  incorporate  this  effect  by  introducing  the  vari¬ 
able  psO^  which  represents  the  probability  that  a  dike 
will  actually  erupt  to  form  a  lava  flow  on  the  seafloor. 
Any  extension  of  the  crust  due  to  tectonics  is  consid¬ 
ered  negligible  in  this  model.  At  fast  spreading 
ridges  faulting  accounts  for  only  3-10%  of  the  total 


spreading  rate  [49,50],  To  match  a  full  spreading  rate 
of  150  mm/yr  (southern  EPR)  aim  thick  dike  must 
be  intruded  on  average  every  7  years,  assuming  that 
all  the  extension  is  accommodated  by  construction  of 
the  sheeted  dike  section.  At  a  full  spreading  rate  of 
1 10  mm/yr  (northern  EPR)  a  diking  event  occurs  on 
average  every  9  years.  Note,  however,  that  diking 
events  are  likely  to  be  clustered  in  time  at  any 
particular  location  with  intervening  quiescent  peri¬ 
ods. 


Fig.  4.  Tte  change  in  sha{>e  of  the  2A/2B  boundary  when  varying  the  stMidard  deviation  of  dike  intrusion  a,  ,  and  the  standard  deviation  of 
lava  ^placement,  cr^,  fcw^  die  unimodai  lava  emplacement  model.  The  curves  shown  are  for  averages  of  50  model  runs  which  are 
subsequently  filtered  widi  a  diangular  filter  of  200  m  wavelength,  F<m*  these  nms,  psO  =  0.4,  r  =  1  m,  and  w  =  1  m.  (a)  The  intrusion  width 
is  varied,  tXj  =  1-400  m  m  100  m  increments,  while  dw  flow  Icugdi  is  held  fixed,  <j^  =  500  m.  (b)  The  flow  length  is  varied,  =  200-600 
in  100  m  increments,  while  the  mtrusion  width  remains  fixed,  o,  =  300  m.  The  standard  deviation  of  the  mean  of  the  50  model  runs,  a,,  is 
indicated  as  a  measure  of  the  variability  in  the  profiles,  which  reflects  the  diickness  of  die  transidon  zcwic  from  20%  to  70%  sheeted  dikes. 
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The  crustal  structure  generated  by  our  model  gen¬ 
erates  a  highly  variable  contact  between  the  extru- 
sives  and  sheeted  dikes  and  predicts  tfiat  isochrons 
within  the  extrusive  section  dip  towards  the  ridge 
axis  (Fig.  3a).  Fig.  3a  shows  the  result  of  over 
10,000  intrusion  and  emplacement  events  to  form  an 
upper  crustal  section  10  km  wide.  To  match  seismic 
profile  SEPR31,  the  model  was  run  using  =  300 
m,  and  =  500  m,  with  40%  of  the  dikes  erupting. 
On  average,  as  they  are  rafted  away  from  the  axis, 
the  tips  of  the  earlier  flows  are  not  covered  by 
subsequent  flows,  resulting  in  isochrons  within  the 
upper  crustal  section  which  dip  towards  the  ridge 
axis.  The  contact  between  sheeted  dikes  and  extru- 
sives  in  this  model  is  extremely  variable.  To  simulate 
a  seismic  observation  of  the  extrusive-sheeted  dike 
contact  that  is  generated  by  this  model,  we  have 
smoothed  the  contact  using  a  triangular  filter  with  a 
width  of  one  seismic  wavelength  (Fig.  3b).  At  a 
dominant  frequency  of  20  Hz  and  crustal  velocity  of 
4  km/s,  the  seismic  wavelength  is  200  m.  One 
seismic  wavelength  is  an  approximation  of  the  dis¬ 
tance  required  for  the  energy  to  turn  within  layer  2A 
and  generate  a  wide-angle  ‘reflection ’/refraction. 
The  variation  in  extrusive  layer  thickness  predicted 
by  the  model  compares  well  with  seismic  observa¬ 
tions  of  layer  2A  thickness  variations  at  the  EPR  in 
terms  of  the  on-axis  (200  m)  and  off-axis  (510  m) 
thicknesses  and  the  distance  over  which  this  thicken¬ 
ing  occurs  (2300  m  half-width).  However,  the  extru¬ 
sive-dike  transition  thickness  calculated  from  the 
model  is  310  m,  wider  than  most  geological  esti¬ 
mates  (Table  2).  If  the  model  is  run  with  the  same 
but  a,  =  10  m  (Fig.  3c,d),  the  transition  zone  from 
extrusives  to  sheeted  dikes  is  50  m,  which  is  compa¬ 
rable  to  ophiolite  and  oceanic  scarp  observations. 
However,  now  the  shape  of  the  layer  2A  boundary  is 
pointed  rather  dian  flat-topped  with  a  very  small 
on-axis  thickness  (70  m),  which  does  not  resemble 
the  seismic  profile.  When  the  extrusive— dike  transi¬ 
tion  thickness  is  included  as  a  constraint,  there  is  no 
longer  a  solution  of  the  unimodal  lava  emplacement 
model  which  satisfies  the  four  seismic  and  geologi¬ 
cal  observables. 

5.2.  Sensitivity  of  the  model 

We  examine  the  effect  of  changing  the  character¬ 


istic  dimensions  of  dike  intrusion  and  lava  emplace¬ 
ment  on  the  four  observables  that  we  can  measure: 
the  on-cods  extrusive  layer  thickness^  the  off-axis 
extrusive  layer  thickness,  die  width  of  2A  accretion 
and  the  transition  thickness  from  20%  to  70%  sheeted 
dikes.  The  unimodal  lava  emplacement  model  is 
characterized  by  five  parameters  (Table  3):  the  stan¬ 
dard  deviation  of  dike  intrusion,  a,-,  die  standard 
deviation  of  the  extrusive  pile,  a^,  the  probability 
that  a  dike  will  erupt  a  lava  flow,  psO,  the  thickness 
of  the  lava  flows,  t,  and  the  width  of  the  dikes,  w. 
Fig.  4  shows  the  effect  of  varying  the  two  parame¬ 
ters,  and  on  die  shape  of  the  2A/2B  boundary 
for  an  average  of  50  model  runs  which  are  subse¬ 
quently  filtered  with  a  200  m  wavelength  triangular 
filter.  For  these  mns,  psO  =  0.4,  r  =  1  m,  and  w  —  1 
m.  The  parameter  is  the  standard  deviation  of  die 
mean  of  the  50  model  runs  and  reflects  the  thickness 
of  the  transition  zone  from  20-70%  sheeted  dikes.  In 
Fig.  4a,  is  held  fixed.  If  a,,  is  extremely  narrow, 
then  the  2A/2B  boundary  will  have  a  pointed  shape 
with  a  very  thin  layer  2A  at  the  rise  axis  and  a  very 
narrow  transition  thickness,  or,  (Fig.  3c,d).  If  a,  is 
wider,  then  the  on-axis  thickness  of  layer  2A  and  the 
thickness  of  the  transition  zone  increase,  and  the 
2A/2B  boundary  has  a  broader,  flat-topped  shape. 
The  off-axis  thickness  remains  constant,  since  is 
constant,  and  thus,  when  the  statistical  variability  of 
the  model  is  smoothed  out,  the  total  volume  of 
extrusives  remains  the  same.  In  Fig.  4b,  cr,  is  held 
constant  and  is  increased.  This  results  in  an 
increase  in  the  off-axis  thickness  of  the  extrusive 
layer,  which  is  accompanied  by  an  increase  in  on-axis 
thickness.  Longer  flows  also  result  in  layer  2A  thick¬ 
ening  over  a  greater  distance  off-axis  and  lead  to  a 
small  increase  in  the  thickness  of  the  transition  zone, 
o-,. 

The  off-axis  thickness  of  the  lava  pile  can  be 
determined  from  mass  conservation  and  is  indepen¬ 
dent  of  o-f.  The  off-axis  extrusive  thickness  is: 


t'psO)  " 


(I) 


The  off-axis  extrusive  thickness  depends  on  the 
mean  lava  flow  length,  scales  linearly  with  the  flow 
thickness  and  with  the  percentage  of  dikes  that  erupt, 
and  scales  inversely  with  dike  width.  Fig.  5a  shows 
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a) 

Off-axis  Extrusive  Thickness,  m 


b) 

Ornaxis  Extrusive  Thickness,  m 


c) 

Transition  Zone  Thickness,  m 


Fig.  5.  Contour  plots  of  the  values  of  the  on-axis  extrusive 
thickness,  off-axis  extrusive  thickness,  and  transition  thickness  as 
a  function  of  a,  and  for  die  unimodal  lava  emplacement 
model.  Values  are  determined  from  the  mean  of  five  runs  of  the 
emplacement  model  witfi  o,  ranging  from  10  to  400  m  and 
ranging  from  100  to  400  m.  For  these  iterations  the  scaling 
variables  are  held  fixed;  /  =  1  m,  w  =  1  m,  and  =  1.  The 
tm-axis  thickness  of  die  extrusive  layer  and  die  extrusive-dike 
transidem  thickness  are  soisitive  to  bodi  a,  and  while  die 
off-axis  extrusive  layer  thickness  depends  only  on  cr^. 


the  off-axis  extrusive  thickness  as  a  function  of  the 
two  stochastic  variables,  a,-  and  determined  from 
the  mean  of  five  runs  of  the  emplacement  model, 
with  or,  ranging  from  10  to  400  m  and  ranging 
from  100  to  400  m.  The  cm-axis  extrusive  thickness 
depends  on  the  ratio  of  the  two  standard  deviations 
and  scales  in  the  same  way  as  the  off-axis  extrusive 
thickness: 


The  on-axis  extrusive  thickness  is  contoured  as  a 
function  of  a,,  and  in  Fig.  5b.  As  the  intrusive 
width  and  the  flow  length  tend  to  zero,  the  on-axis 
thickness  goes  to  zero.  For  an  infinitely  wick  intru¬ 
sion  zone,  the  extrusive  thickness  on  axis  is  ^ual  to 
that  off  axis.  The  thickness  of  the  transition  zone 
from  extrusives  to  sheeted  dikes  is  shown  in  Fig.  5c 
and  has  a  similar  dependence  on  f,  psO  and  w 
as  the  on-axis  extrusive  thickness.  Note  that,  for  the 
definition  that  we  have  chosen  (20-70%  sheeted 
dikes),  the  transition  width  is  always  greater  tiian  or 
equal  to  the  on-axis  extrusive  thickness.  The  width 
of  the  2A  accretion  zone  is  a  less  stable  measure¬ 
ment  and  much  more  sensitive  to  the  statistical  vari¬ 
ability  in  the  model.  As  can  be  seen  from  Fig.  4,  the 
width  of  the  2A  acen'etion  zone  is  positively  related 
to  but  also  varies  with  a,-.  The  width  of  the  2A 
accretion  zone  is  independent  of  r,  but  depends  on 
psO  and  w. 

This  stochastic  emplacement  model  with  a  uni¬ 
modal  distribution  of  flows  matches  the  variation  in 
layer  2A  thickness  observed  seismically  at  the  EPR 
(Fig.  3)  and  provides  some  intuitive  insight  into  the 
processes  that  control  extrusive  thickness  at  fast 
spreading  ridges.  For  example,  this  model  shows  that 
the  distance  off-axis  at  which  layer  2A  reaches  its 
final  thickness  is  primarily  affected  by  the  distance 
flows  extend  off-axis.  For  a  unimodal  distribution  of 
flows,  the  shape  of  the  layer  2A/2B  boundary  (flat- 
topped  rather  than  triangular)  requires  a  zone  of 
intrusion  on  the  order  of  a  few  hundred  meters. 
However,  this  unimodal  lava  emplacement  model 
cannot  explain  both  the  broad,  flat-topped  shape  of 
the  layer  2A/2B  boundary  and  at  the  same  time  the 
sharpness  of  the  extrusive-dike  transition  measured 
in  outcrop  and  drill  holes.  The  wider  intrusion  zone 
required  to  match  the  seismic  observations  will  result 
in  a  broad  transition  zone,  due  to  the  variability  in 
location  of  dike  intrusicMi.  Indeed,  for  this  unimodal 
lava  emplacement  model  the  transition  thickness  is 
greater  than  or  equal  to  die  on-axis  layer  2A  thick¬ 
ness,  which  at  the  EPR  is  about  250  m  (Fig.  5).  The 
transition  thickness  predicted  for  the  model  shown  in 
Fig.  3,  which  fits  the  observed  layer  2A  thickness,  is 
over  300  m.  This  is  larger  than  transition  thicknesses 
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of  a  few  tens  of  meters  to  100-200  m  measured  in 
ophiolites  [1,42]  and  150-200  m  drilled  at  Hole 
504B  [16]  (Table  2).  To  match  the  observed  sharp¬ 
ness  of  the  extrusive-dike  transition  thickness,  the 
standard  deviation  of  intrusion  must  be  narrower 
than  the  mean  flow  length  by  about  an  order  of 
magnitude  [1],  Solutions  to  the  unimodal  lava  em¬ 
placement  model  which  match  the  shape  of  the 
seismic  2A/2B  boundary  have  <y.  and  of  the 
same  order  of  magnitude.  In  order  to  match  the 
sharpness  of  the  extrusive-dike  transition  zone  ob¬ 
served  in  ophiolites,  a,  must  be  of  the  order  of 
10-50  m.  However,  this  will  result  in  a  very  small 
on-axis  2A  thickness  and  a  boundary  which  is  peaked 
rather  than  flat-topped  (Fig.  3c,d),  contrary  to  what 
is  observed  in  the  seismic  data. 

3.3,  Bimodal  lava  emplacement  model 

We  propose  a  stochastic  emplacement  model  with 
a  bimodal  distribution  of  flows  in  order  to  match 
both  the  dimensions  and  shape  of  the  seismic  bound¬ 
ary  and  the  geologically  observed  extrusive-dike 
transition  thickness.  The  idea  underlying  the  bimodal 
lava  emplacement  model  is  to  rapidly  build  up  the 
on-axis  thickness  of  extrusives  by  accumulating  a 
sequence  of  short  flows  erupted  at  the  ridge  crest 
(standard  deviation  of  emplacement,  0^5).  The  more 
progressive  thickening  of  the  extrusive  layer  over 
1-3  km  observed  in  the  seismic  sections  is  accom¬ 
plished  by  a  series  of  flows  that  are  about  as  long  as 


the  width  of  the  region  of  thickening  (standard  devia¬ 
tion  of  emplacement,  Geologically,  the  short 
flows  remain  confined  within  the  axial  summit 
caldera  (ASC).  The  long  flows  flow  down  the  flanks 
of  the  rise  crest,  are  channelled  across  the  shoulders 
of  the  rise  in  lava  tubes,  or  erupt  outside  the  ASC 
(Fig.  6).  The  upper  part  of  the  resulting  extrusive 
section  is  composed  of  flows  deposited  away  from 
the  axis  while  the  lower  portion  consists  of  flows 
deposited  on-axis.  Since  the  on-axis  thickness  of 
layer  2A  is  approximately  half  of  the  off-axis  thick¬ 
ness,  the  on-axis  and  off-axis  flows  must  contribute 
equal  volumes  to  the  extrusive  section.  To  accom¬ 
plish  this,  if  the  flows  all  have  the  same  thickness, 
the  short  flows  must  be  erupted  much  more  fre¬ 
quently  than  the  long  flows  (95%  of  the  time).  In  this 
case,  the  mean  volume  of  the  long  flows  is  -20 
times  that  of  the  short  flows.  On  the  other  hand,  the 
relative  number  of  events  decreases  if  the  short  flows 
are  much  thicker  than  the  long  flows,  possibly  due  to 
the  ponding  of  these  flows  against  the  walls  of  the 
ASC.  A  flow  thickness  /  =  2  m  implies  flows  with  a 
cross-sectional  area  of  2,000  m^  for  the  long  flows. 
To  compare  this  2-D  flow  volume  with  that  of  the 
mid-crustal  melt  sill,  we  assume  a  melt  sill  thickness 
of  50  m  and  a  sill  width  of  800  m  [5 1  ].  The  melt  sill 
has  a  volume  of  4,000  m^,  twice  that  needed  to  form 
the  rare,  voluminous,  long  flow  length  eruptions. 

A  model  to  build  the  extrusive  section  of  the 
ocean  crust  that  combines  an  initially  narrow  and 
subsequently  wide  zone  of  accretion  was  suggested 


Rg.  6.  Diagram  showing  the  bimodal  flow  length  emplacement  model.  Short  flows,  shown  in  dark  gray,  remain  confined  within  the  axial 
summit  caldera  (ASC)  and  build  up  half  the  volume  of  the  extrusive  section.  Occasional  long  flows,  shown  in  light  gray,  are  envisioned  to 
flow  down  the  flanks  of  the  rise  crest,  be  channeled  across  the  shoulders  of  the  rise  in  lava  tubes,  or  erupt  outside  the  ASC.  They  flow  far 
off-axis  and  build  up  the  remainder  of  the  extrusive  section.  The  accompanying  dikes  arc  shown  in  white. 
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Bimodal  Lava  Emplacement  Model 
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Fig,  7.  Two  examples  of  nms  of  tlw  bimodal  lava  emplacement  model  of  upper  crustal  construction  that  match  the  layer  2A  thicknesses  for 
the  across-axis  seismic  lines  SEPR06  (wide  region  of  2A  accretion)  and  SEPR31  (narrow  region  of  2A  accretion)  shown  in  Fig.  1.  (a)  Run 
of  the  bimodal  lava  emplacement  model  that  matches  SEPR06;  standard  deviation  of  mtrusion,  a,-  =  10  m,  lava  flow  distribution,  95%  of  the 
flows  have  0^5  =  75  m,  5%  of  the  flows  have  —  1500  m,  r  =  1.5  m.  For  the  long  flows  the  distance  over  which  no  lava  accumulates  is 
1000  m.  (c)  Run  of  the  bimodal  lava  emplacement  model  that  matches  SEPR31;  standard  deviation  of  intrusion,  a,  =  10  m,  lava  flow 
distribution,  95%  of  the  flows  have  =  50  m,  5%  of  the  flows  have  ==  1000  m,  f  =  2.2  m.  For  die  Icmg  flows  the  distance  over  which 
no  lava  accumulates  is  400  m.  In  both  figures  die  isochrons  in  the  extrusive  section  are  shown  with  a  gray  scale,  sheeted  dikes  are  in  whitc. 
No  verticad  exaggeration.  Note  d^  diallow  dip  of  die  upper  extrusive  section  (0-10®  and  15®)  and  the  steep  dips  in  the  lower  exuusives, 
50-70®.  The  boundary  between  die  two  units  is  sharp  and  dips  are  zero  at  the  contact  (b)  and  (d)  The  gray  line  is  the  boundary  between  the 
sheetoi  dikes  and  extrusives  for  the  two  model  runs.  The  hollow  line  shows  the  model  extrusive-sheeted  dike  contact  averaged  using  a  200 
m  wavelength  triangular  filter.  The  black  line  is  the  layer  2A  thicknesses  determined  from  (b)  line  SEPR06  and  (d)  line  SEPR3 1 .  The 
on-axis,  off-axis  thicknesses  and  die  widdis  of  the  2A  accretion  zone  match  well.  The  transition  zone  thickness  for  these  models  is  90  m  and 
70  m,  respectively,  compatible  with  measurements  made  on  ophiolites  and  oceanic  crustal  scarps.  Note  that  tiicsc  fits  arc  not  unique,  but 
provide  a  characterizatitMi  of  die  differences  and  types  of  length  scales  required  to  fit  the  observables  available  to  constrain  crustal  accretion. 
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by  Atwater  [52]  to  reconcile:  (a)  the  observation  of  a 
wide  emplacement  2X)ne  suggested  by  the  presence 
of  intact  volcanoes  on  the  surface  of  the  off-axis 
seafloor  in  the  FAMOUS  area;  and  (b)  ophiolite 
observations  of  a  very  localized  transition  from  ex- 
trusives  to  sheeted  dikes.  Our  attempt  to  constrain 
the  accretion  process  using  the  geometry  and  dimen¬ 
sions  of  the  seismic  layer  2A  boundary  and,  at  the 
same  time,  match  the  transition  thickness  observed  in 
ophiolites  has  also  led  us  to  require  a  bimodal  lava 
emplacement  model  for  Ae  fast  spreading  EPR.  Our 
ability  to  image  seismically  layer  2A  thickness  varia¬ 
tions  enables  us  to  constrain  the  relative  volumes 
contributed  by  these  two  modes  of  extrusion  and  the 
length  scales  involved. 

In  the  bimodal  lava  emplacement  model,  the  nar¬ 
row  intrusion  zone  and  short  flows  generate  the 
on-axis  extrusive  pile  and  result  in  a  narrow  transi¬ 
tion  zone  from  extrusives  to  sheeted  dikes.  The 
off-axis  flows  thicken  the  crust  to  the  final  thickness 
over  the  distance  of  a  few  kilometers.  However,  the 
shape  of  the  extrusive-sheeted  dike  boundary  ob¬ 
tained  remains  pointed  rather  than  flat-topped.  To 
obtain  a  flat-topped  2A/2B  boundary,  as  is  observed 
in  the  across-axis  seismic  lines  (Fig.  1),  zero  eruptive 
volume  from  the  long  flows  must  be  accumulated 
over  the  flat  region.  A  distribution  of  long  flows  that 
have  zero  thickness  until  some  distance,  m/n/,  from 
the  axis  (Fig.  2c)  gives  an  inflection  point  in  the 
cumulative  lava  volume.  This  simulates  eruptions 
that  flow  over  the  shoulders  of  the  ASC,  leaving  a 
negligible  thickness  of  lava  until  they  pond  in  deeper 
areas  off-axis.  Alternatively,  this  represents  flows 
that  are  channelled,  perhaps  in  subseafloor  lava  tubes, 
from  the  axis  across  the  rift  shoulders  to  an  eruption 
site  off-axis  [14],  or  voluminous  flows  that  are 
erupted  off-axis  from  a  small  number  of  dikes  (Fig. 
6).  Once  a  flow  leaves  the  summit  region  there  is 
little  to  impede  its  flow  until  it  reaches  the  base  of 
die  axial  high  3-4  km  off-axis.  The  variables  used  in 
the  bimodal  lava  emplacement  model  are  listed  in 
Table  3.  To  keep  the  number  of  variables  to  a 
minimum,  we  have  kept  the  same  for  both  the 
long  and  the  short  lava  flows  (a,  ~  10  m)  and  we 
fixed  psO  =  1  and  w  =  1  m.  As  in  the  unimodal  lava 
emplacement  model,  the  effect  of  decreasing  the 
thickness  of  the  lava  flows  or  decreasing  the  amount 
of  dikes  that  erupt  a  lava  flow  is  to  decrease  the 


thickness  of  the  extrusive  pile  and  of  the  transition 
zone.  Decreasing  the  amount  of  dikes  that  erupt  a 
lava  flow  also  broadens  the  width  of  the  2A  accre¬ 
tion  zone. 

Fig.  7  shows  two  examples  of  the  bimodal  lava 
emplacement  model  that  match  the  two  across-axis 
seismic  lines  which  represent  extremes  of  the  dis¬ 
tance  over  which  die  off-axis  thickness  is  attained, 
SEPR06  (broad  2A  accretion  zone)  and  SEPR31 
(narrow  2A  accretion  zone),  shown  in  Fig.  1.  Solu¬ 
tions  to  both  the  unimodal  and  the  bimodal  lava 
emplacement  models  are  non-unique,  but  the  crustal 
accretion  models  that  we  present  provide  a  character¬ 
ization  of  the  differences  and  types  of  length  scales 
required  to  match  the  available  observations.  To  fit 
layer  2A  thicknesses  for  line  SEPR06  (Fig.  7a,b),  the 
model  had  a  standard  deviation  of  intrusion,  =  10 
m;  lava  flow  length  distribution,  95%  of  the  flows 
have  0^5  ==  75  m,  5%  of  the  flows  have  =  1500 
m;  and  die  thickness  of  the  lava  flows,  r  =  1 .5  m.  As 
a  simplification,  we  assume  that  all  dikes  erupt  a 
lava  flow,  since  there  is  a  trade-off  between  lava 
flow  thickness  and  percentage  of  dikes  erupting 
flows.  For  the  long  flows,  the  distance  from  the  axis 
over  which  no  off-axis  lava  accumulates,  minl  = 
1000  m.  This  is  the  distance  from  the  axis  over 
which  layer  2A  is  not  accreting  new  material  beyond 
what  was  built  up  on  the  axis  (Fig.  lb).  Averaging 
the  geological  contact  obtained  from  this  model  with 
a  200  m  wavelength  triangular  filter  gave  model 
values  of  on-axis  thickness  —  140  m,  off-axis  thick¬ 
ness  —  460  m,  width  of  2A  accretion  zone  =  2170  m, 
and  transition  zone  thickness  =  90  m.  A  good  match 
was  also  obtained  for  line  SEPR31  (Fig.  7c,d).  The 
model  parameters  used  in  this  case  were:  standard 
deviation  of  intrusion,  =  10  m;  lava  flow  length 
distribution,  95%  of  the  flows  have  0-^5  =  50  m,  5% 
of  the  flows  have  =  1000  m;  and  the  thickness 
of  the  lava  flows,  t  —  2.2  m.  For  the  long  flows  the 
distance  over  which  no  off-axis  lava  accumulates, 
mini  —  400  m  —  layer  2A  does  not  thicken  beyond 
the  on-axis  thickness  over  this  distance  from  the 
ridge  axis  (Fig.  Id).  The  filtered  sheeted  dike-extru¬ 
sive  boundary  had  an  on-cods  thickness  =  160  m, 
off-axis  thickness  —  500  m,  width  of  2A  accretion 
zone  =  1700  m  and  a  transition  thickness  —  70  m. 
For  these  two  model  runs  the  filtered  extrusive— dike 
contact  matches  the  seismic  layer  2A  thickness  pro- 
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a)  Blmodal  Lava  Emplacenier^  Model  RItered  Velocity  Structure;  wavelength  =  200  m 
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files  well,  both  in  terms  of  its  shape,  variability,  and 
extrusive  thickness  values.  In  addition,  the  transition 
thicknesses  calculated  from  the  bimodal  lava  em¬ 
placement  model  are  comparable  to  those  measured 
in  ophiolites  and  oceanic  cmstal  sections.  The  differ¬ 
ences  in  the  width  of  the  zone  of  2A  accretion  for 
the  two  models  is  consistent  with  the  observed  varia¬ 
tion  in  across-axis  topography  for  the  two  seismic 
profiles  (Fig.  la,c).  SEPR06  has  a  broad  axial  high 
and  thus  off-axis  lava  flows  do  not  pond  till  further 
from  the  ridge  axis,  resulting  in  a  broader  zone  of 
thickening.  On  the  other  hand,  SEPR31  shows  a 
subdued  narrow  axial  high  and  thus  long  lava  flows 
rapidly  reach  relatively  flat  seafloor  and  pond  closer 
to  the  rise  axis,  causing  earlier  thickening  of  the  lava 
pile  and  a  narrow  2A  accretion  zone. 

4.  Implications  of  the  model 

In  this  section  we  explore  the  implications  of  the 
stochastic  emplacement  model  for  measured  seismic 
velocities  in  the  upper  crust,  magnetic  anomaly  tran¬ 
sition  widths,  and  seafloor  lava  ages.  From  the  model 
we  determined  the  velocity  structure  of  the  upper 
crust,  the  dip  of  isochrons  with  depth,  the  age  of 
seafloor  lavas  relative  to  the  underlying  crust,  and 
magnetic  transition  widths.  The  thickness  of  the 
shallow,  steep  velocity  gradient  reflects  the  thickness 
of  the  transition  region  from  extrusives  to  sheeted 
dikes  and  the  inherent  resolution  of  the  seismic 
technique  used.  The  length  of  the  lava  flows  and 
their  subsidence  rates  affect  the  dip  of  isochrons  in 


the  extrusive  section,  the  width  of  the  transition 
region  observed  in  magnetic  anomalies,  as  well  as 
the  age  of  lavas  on  the  seafloor. 

4.1.  Velocity-depth  profiles  of  the  upper  crust 

In  order  to  see  if  the  emplacement  model  that  we 
propose  is  consistent  with  the  shallow  velocity  struc¬ 
ture  observed  at  the  EPR,  we  have  calculated  a 
synthetic  velocity  structure  based  on  the  distribution 
of  extrusives  and  dikes  predicted  from  the  emplace¬ 
ment  model.  The  shallow  velocity  structure  of  the 
crust  is  primarily  determined  by  bulk  porosity.  Frac¬ 
tured  porous  lavas  have  a  higher  porosity  than  mas¬ 
sive  sheeted  dikes  so  the  relative  abundance  of  lavas 
and  dikes  in  a  vertical  section  will  have  a  major 
effect  on  the  velocity  structure.  We  assume  that  the 
velocities  in  the  extrusive  and  sheeted  dike  sections 
are  constant  and  equal  to  2,2  km/s  and  5.5  km/s, 
respectively.  The  constant  low  velocity  extrusive 
section  results  in  a  surficial  low  velocity  zone.  The 
underlying  high  gradient  region  is  due  to  the  mixed 
zone  where  dikes  and  extrusives  interfinger.  This  is 
an  approximation  which  does  not  take  into  account 
the  probable  variation  of  velocity  with  depth  within 
the  extrusive  pile  that  results  as  this  layer  gets  buried 
and  porosity,  decreases  due  to  the  closing  of  cracks 
and  voids.  Differences  due  to  any  morphological  or 
chemical  differences  between  the  on-axis  and  the 
off-axis  flows  (i.e.  the  upper  and  lower  extrusives) 
are  also  ignored. 

The  first  area  where  we  compare  the  measured 
with  the  predicted  seismic  velocity  structure  is  9°35'N 


Fig-  8-  A(a)  Velocity  structure  for  a  bimodal  lava  emplacement  model  run  that  fits  the  layer  2A  thicknesses  measured  on  the  wide  aperture 
profile  at  9^35"  N  EPR  [241.  <r^  =  10  m,  =  60  m.  =  1200  m,  /  =  I  m,  mini  =  1000  m.  Velocities  of  2.2  km/s  (dark  gray)  and  5,5 
km/s  (light  gray)  are  assigned  to  the  extrusive  and  sheeted-dike  sections  respectively.  A  velocity  filter  (a  cone  of  radius  one  wavelength. 
200  m)  was  applied.  The  base  of  the  steep  velocity  gradient  is  similar  to  that  obtained  by  the  linear  triangular  filter  of  the  geological 
boundary  (white  line),  (b)  and  (c)  Depth  sections  of  the  velocity  model  on-  and  off-axis  (70  kyr)  (solid  line),  their  locations  are  shown  in  (a) 
as  dark  vertical  lines.  These  arc  compared  to  the  velocity -depth  profiles  determined  at  this  latitude,  ESP5  (off-axis.  180  kyr)  and  ESPl 
(on-axis)  (dashed  line)  (20].  Note  the  similarity  in  the  thickness  of  the  surficial  low-velocity  zone  and  the  depth  over  which  the  steep 
velocity  gradient  extcnds.B(a)  Velocity  structure  for  a  bimodal  lava  emplacement  model  run  that  fits  the  layer  2A  thicknesses  measured  on 
the  CDP3 1  at  9®30'N  EPR  (281.  =  10  m,  =  75  m.  =  1500  m.  /  =  1  m.  mini  —  1000  m.  Velocities  of  2.2  km /s  (dark  gray)  and  5.5 

km/s  (light  gray)  arc  assigned  to  the  extrusive  and  sheeted  dike  sections,  respectively.  A  velocity  filter  (a  cone  of  radius  one  seismic 
wavelength,  67  m)  was  applied.  The  base  of  the  steep  velocity  gradient  is  similar  to  that  obtained  by  the  linear  triangular  filter  of  the 
geological  boundary  (white  line),  (b)-(d)  Depth  sections  of  the  velocity  model  on-axis  and  at  10  and  20  kyr  (solid  line),  their  locations  are 
shown  in  (a)  as  dark  vertical  lines.  These  arc  compared  to  the  velocity-depth  profiles  determined  at  this  latitude  from  the  on-bottom 
refraction  (NOBEL)  experiment  (36).  Note  the  similarity  in  the  thickness  of  the  surficial  low- velocity  zone  and  the  depth  over  which  the 
steep  velocity  gradient  extends. 
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on  the  EPR,  The  seismic  velocity  structure  at  this 
location  is  constrained  by  several  expanding  spread 
profiles  (ESP)  [20].  The  geometry  of  the  layer  2 A 
boundary  is  constrained  from  a  two-ship  wide  aper¬ 
ture  profile  (WAP)  shot  here  [24].  We  generate  a 
model  using  the  bimodal  distribution  of  flows  that 
fits  the  layer  2A  geometry  in  this  area  and  assign 
velocities  of  2,2  km/s  and  5.5  km/s  to  the  extni- 
sives  and  sheeted  dikes,  respectively.  The  source 
used  for  the  ESPs  had  a  dominant  frequency  of  20 
Hz  corresponding  to  a  wavelength  of  200  m  at  4 
km/s.  We  therefore  averaged  our  velocity  structure 
using  a  cone-shaped  filter  with  radius  200  m  (Fig. 
8Aa)  in  order  to  compare  it  with  the  seismically 
derived  velocity  structure.  The  base  of  the  resulting 
high  gradient  region  closely  resembles  the  result  of 
running  a  one-dimensional  triangular  filter  over  the 
extrusive-dike  contact.  Velocity -depth  profiles  were 
constructed  at  ages  of  0  kyr  and  70  kyr  (solid  lines  in 
Fig.  8A).  They  compare  very  well  with  the 
velocity-depth  profiles  obtained  from  the  ESPs 
(dashed  lines).  There  is  no  surficial  low- velocity 
zone  in  the  model  for  0  kyr,  while  one  is  observed  in 
ESP5.  ESPS  may  be  located  slightly  off-axis,  which 
would  account  for  this  discrepancy.  ESPl  actually 
lies  on  180  kyr  old  crust  but  velocity  profiles  at  180 
kyr  and  70  kyr  are  assumed  to  represent  off-axis 
sections  where  crustal  construction  has  been  com¬ 
pleted.  The  thickness,  shape,  and  slope  of  the  high 
gradient  region  is  very  similar  to  those  measured 
with  the  ESPs. 

The  same  procedure  was  applied  to  a  second  area 
at  9*^30' N  on  the  EPR,  where  the  upper  crustal 
velocity  structure  is  constrained  by  a  high-resolution 
on-bottom  refraction  experiment  [36].  The  geometry 
of  layer  2A  is  determined  from  MCS  line  CDP31 
[28].  The  on-bottom  refraction  experiment  used  a 
source  with  a  dominant  frequency  of  60  Hz.  A  model 
was  constructed  to  matched  the  line  CDP31  and  was 
filtered  at  a  wavelength  of  67  m  (Fig.  8Ba),  Again, 
depth  sections  through  the  resulting  velocity  model 
show  a  high  velocity  gradient  very  similar  in  thick¬ 
ness,  shape,  and  slope  to  those  measured  in  the 
on-bottom  refraction  experiment  on  crust  0,  10  and 
20  kyr  (Fig.  8Bb-d). 

The  fact  that  the  boundary  between  the  sheeted 
dikes  and  extrusives  is  not  seismically  sharp  is  con¬ 
sistent  with  the  multichannel  data.  The  base  of  layer 


2A  is  not  imaged  as  a  normal  incidence  reflection  in 
the  multichannel  seismic  data  but,  rather,  is  a  diving 
ray  which  turns  in  a  steep  velocity  gradient  and  is 
observed  at  larger  offsets  on  the  streamer  as  a  wide- 
angle  ‘reflection ’/refraction.  These  observations 
support  the  idea  that  the  boundary  between  the  extru¬ 
sives  and  sheeted  dikes  is  not  sharp  on  the  scale  of  a 
quarter-wavelength  ( --  50  m).  However,  on  the  scale 
of  a  full  wavelength  a  coherent  steep  velocity  gradi¬ 
ent  is  sensed,  causing  the  rays  to  turn  back  to  the 
surface  and  resulting  in  a  wide-angle  reflection. 
Therefore,  in  this  case,  one  seismic  wavelength  is  an 
estimate  of  the  distance  required  for  energy  to  refract 
in  a  steep  velocity  gradient  and  produce  a  turning 
ray. 

The  thickness  of  the  seismic  high  velocity  gradi¬ 
ent  predicted  from  the  bimodal  lava  emplacement 
model  depends  on  the  thickness  of  the  transition 
zone  from  extrusives  to  sheeted  dikes  and  on  the 
dominant  wavelength  of  the  experiment  that  is  mea¬ 
suring  the  velocity  structure.  For  a  model  with  a 
wide  transition  region,  as  would  be  obtained  for  the 
unimodal  lava  emplacement  model,  a  shorter  seismic 
wavelength  will  give  a  sharp  boundary  but,  laterally, 
it  will  be  highly  variable.  A  long  seismic  wavelength 
will  average  out  the  variability  and  give  a  thicker 
high  velocity  gradient  than  has  been  seismically 
observed.  For  a  narrow  transition  region,  the  long 
wavelength  energy  will  simply  see  a  somewhat 
broader  transition  region  than  the  shorter  wave¬ 
length.  Under  the  simple  assumptions  that  we  make 
here  the  velocity  structure  predicted  from  the  bi¬ 
modal  lava  emplacement  model  is  consistent  with 
the  shallow  crustal  velocity  structure  measured  at  the 
EPR.  At  the  same  time,  the  higher  resolution  on-bot¬ 
tom  refraction  experiment  resolved  a  thinner  on-axis 
thickness  of  layer  2A  than  the  MCS  experiment  (Fig. 
8Be  and  Ad,  respectively),  consistent  with  the  bi¬ 
modal  lava  emplacement  model. 

4.2.  Isochron  dips  and  magnetic  anomaly  transition 
widths 

Lavas  erupted  on  the  seafloor  arc  rapidly  buried 
by  subsequent  flows  and  subside.  In  the  unimodal 
lava  emplacement  model  subsidence  rates  are  20 
mm/yr,  assuming  full  spreading  rates  of  150  mm/yr 
for  the  southern  EPR.  For  the  bimodal  lava  emplace- 
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merit  model,  the  rates  of  subsidence  are  initially  very 
high,  ~  200  mm/yr  for  the  first  1.5  kyr,  and  then 
drop  to  40-10  mm/yr  during  the  following  8-27 
kyr.  The  dips  of  isochrons  are  sensitive  to  the  length 
of  the  lava  flows  and  the  amount  of  subsidence  that 
occurs.  In  the  unimodal  lava  emplacement  model, 
dips  are  relatively  steep  throughout  the  extrusive 
section  and  vary  from  ~  0°,  near  the  surface,  to 
~  30°  at  depth  (Fig.  3).  For  the  bimodal  lava  em¬ 
placement  models  that  we  propose  in  this  paper  (Fig. 
7a,c),  the  upper  lavas,  composed  of  the  long  flows, 
have  very  shallow  dips,  0°  to  ~  10-15°.  At  the 
boundary  between  the  two  lava  types  dips  are  0°.  In 
the  lower  extrusive  section  flow  lengths  are  short 
and  subsidence  rates  are  high  and  the  dips  are  ex¬ 
tremely  steep,  50-70°.  In  Iceland,  lava  flow  dips  are 
zero  near  the  surface  and  increase  with  depth  reach¬ 
ing  6°  1500  m  below  the  original  top  of  the  lava  pile 
[43].  Considering  the  differences  in  extrusive  thick¬ 
ness  between  Iceland  and  the  EPR,  these  observa¬ 
tions  cannot  differentiate  between  both  modes  of 
lava  emplacement.  The  two  models  are  distinguish¬ 
able  because,  in  the  bimodal  lava  emplacement 
model,  the  dips  are  reduced  to  zero  at  the  boundary 
between  the  upper  and  lower  lavas.  H.  Schouten  and 
E.E.E.  Hooft  (Two-layered  structure  of  the  pillow 
lavas  in  oceanic  crust  and  ophiolites,  manuscript  in 
prep.,  1996)  find  evidence  of  a  bimodal  dip  structure 
in  the  paleomagnetic  inclinations  of  several 
DSDP/ODP  drill  cores,  suggesting  that  a  bimodal 
emplacement  of  the  upper  oceanic  crust  also  occurs 
at  intermediate  and  slow  spreading  centers. 

Work  on  the  transition  width  of  marine  magnetic 
anomalies  from  deep-tow  measurements  on  the  Gorda 
Rise  [4]  suggest  that  the  portion  of  the  transition 
width  due  to  constructional  processes  is  1.7  km. 
Near-bottom  magnetic  profiles  from  ridges  in  the 
Pacific  covering  intermediate  to  fast  spreading  rates 
give  an  average  half-width  of  crustal  formation  of 
2-3  km  [8].  Measurements  made  by  Macdonald  et 
al.,  [7]  on  the  EPR  2rN  give  transition  widths  of 
1-1.4  km.  For  our  bimodal  lava  emplacement  model 
the  dips  of  the  isochrons  vary  with  depth  in  the 
extrusive  section.  We  have  calculated  the  magnetic 
anomaly  from  the  geometry  of  the  bimodal  lava 
emplacement  model  and  compared  it  with  the 
anomaly  that  would  be  obtained  assuming  the  simple 
geometries  of  earlier  models.  The  transition  widths 


obtained,  2-3  km,  are  similar  to  those  observed. 
However,  as  mentioned  before,  these  magnetic 
anomaly  reversal  data  cannot  resolve  the  relative 
contributions  of  the  width  of  the  intrusion  zone  and 
lava  flow  distribution  to  the  lava  emplacement  pro¬ 
cess, 

43.  Age  of  seafloor  lavas 

Goldstein  et  al.  [15]  dated  seafloor  lavas  on  a 
transect  across  the  EPR  at  9°3f  N  and  showed  that 
these  lavas  were  at  least  20-70  kyr  younger  than  the 
predicted  age  of  the  crust  based  on  spreading  rate 
(Fig.  9).  They  suggested  that  off-axis  eruptions  are 
probably  the  cause  of  these  anomalously  young  ages, 
implying  the  existence  of  a  comparatively  broad 
(several  kilometers  half-width)  zone  of  intrusion. 
However,  if  a  significant  proportion  of  the  dikes  are 
intruded  over  a  wide  zone,  the  extrusive-dike  transi¬ 
tion  becomes  much  broader  than  is  consistent  with 
geological  observations  (Table  2). 

In  the  emplacement  model  that  we  propose,  the 
intrusion  zone  is  very  narrow  (on  the  order  of  a  few 
tens  of  meters)  but  the  upper  pan  of  the  extrusive 
section  forms  by  flows  that  have  been  transported  a 
long  distance  (up  to  several  kilometers)  off-axis  from 
their  point  of  eruption.  These  flows  will  have  ages 
that  are  younger  than  the  predicted  age  of  the  crust 
based  on  the  spreading  rate.  From  runs  of  the  bi¬ 
modal  lava  emplacement  model  that  match  the  range 
of  the  seismically  observed  variation  in  layer  2A 
thickness  in  the  9°N  area  on  the  EPR  [24,28],  we 
predict  seafloor  ages  will  average  about  ~  30-70 
kyr  younger  than  the  crustal  spreading  age  (Fig.  9). 
This  is  in  agreement  with  the  anomalously  young 
lavas  dated  by  Goldstein  et  al.  [15]  (Fig.  9),  however, 
the  predicted  ages  are  in  general  younger  than  those 
observed.  The  discrepancy  between  the  model  and 
the  dating  can  be  explained  if  a  significant  amount  of 
the  off-axis  thickening  occurs  below  the  surface,  by 
the  intrusion  of  sills.  Another  possibility  is  that,  if 
the  time  from  melt  segregation  in  the  mantle  to 
eruption  on  the  surface  is  substantial,  then  the  time 
since  eruption  is  significantly  overestimated  by  the 
uranium-series  disequilibrium  dating  [15]. 

Due  to  the  narrowness  of  the  extrusive-dike  tran¬ 
sition  and  its  minimal  effect  on  the  seismic  layer  2A 
thickening  profile,  we  infer  that  the  number  of  dikes 
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Seafloor  Lava  Ages 


Fig.  9.  Age  offsets  predicted  by  this  model.  The  solid  lines  shows  the  ages  of  seafloor  lavas  determined  from  two  bimodal  lava  emplacement 
models  that  match  the  range  of  layer  2A  thickening  profiles  measured  in  the  9°N  area  [24,28]  (Fig.  8Aa  and  8Ba).  The  dashed  line  is  the 
predicted  age  of  the  seafloor  based  on  the  assumed  half-spreading  rate  of  55  mm/yr.  The  seafloor  ages  are  younger  than  predicted  because 
young  lavas  flow  out  over  older  crust.  Age  offsets  arc  ~  30-70  kyr.  The  age  offsets  of  20-70  kyr  were  observed  at  9®3rN  by  Goldstein  et 
al.  [15]  (circles  =^^°Th  ages,  squares  ages).  Note  that  the  bimodal  lava  emplacement  model  assumed  a  very  narrow  (~40  m) 

intrusion  zone.  The  anomalously  young  lava  ages  reported  by  Goldstein  et  al.  [15]  do  not  necessarily  require  a  significant  volume  of  eruptive 
activity  outside  the  ASC,  if  lavas  erupted  in  the  ASC  periodically  spill  out  of  the  ASC  and  flow  for  considerable  distances  down  the  flanks 
of  the  axial  high. 


feeding  any  off-axis  eruptions  is  minimal.  The  source 
region  and  crustal  residence  time  of  any  off-axis 
eruptions  may  well  be  different  from  the  on-axis 
eruptions,  and  so  even  small  volumes  of  lavas  erupted 
off-axis  can  explain  the  diverse  compositions  mea¬ 
sured  in  off-axis  sampling  of  seafloor  basalts  [10].  It 
is  important  to  bear  in  mind  that  rare,  and  volumetri- 
cally  insignificant,  lavas  erupting  outside  of  the  zone 
of  lava  emplacement  will  form  a  thin  cover  on  the 
seaflotx  and  have  a  high  probability  of  being  sam¬ 
pled,  Thus,  while  it  is  possible  that  some  lavas  erupt 
off-axis,  extensive  off-axis  eruptions  are  not  neces¬ 
sarily  required  to  explain  the  anomalously  young  age 
of  seafloor  lavas  observed  near  the  EPR. 

The  existence  of  a  narrow  dike  intrusion  zone  is 
consistent  with  the  observation  that  most  hydrother¬ 
mal  activity  along  the  northern  and  southern  EPR  is 
located  within  the  ASC  [9,14].  Hydrothermal  activity 


is  driven  by  a  heat  source,  on  short  time  scales  by  a 
dike,  over  long  times  by  an  underlying  magma  body. 
The  predominance  of  hydrothermal  activity  within 
the  ASC  supports  the  assertion  that  the  majority  of 
dikes  are  intruded  very  close  to  the  ridge  axis.  Lava 
flows  in  the  bimodal  lava  emplacement  model  ex¬ 
tend  up  to  a  few  kilometers  off-axis  and  are  also 
consistent  with  the  observation  that  sediment  cover 
in  the  9'"30'N  EPR  area  is  anomalously  low  up  to 
3-4  km  from  the  ridge  axis,  beyond  which  the 
sediment  thickness  increases  dramatically  [15],  The 
off-axis  seafloor  is  covered  by  the  lavas  that  have 
flooded  out  of  the  ASC.  The  ASC  is  filled  by  the 
more  common  events  —  short  flows.  This  could 
explain  the  observation  that  most  off-axis  seafloor  is 
covered  by  lobate  lava  flows,  while  the  morphology 
within  the  ASC  is  more  varied,  consisting  of  pillows 
and  sheet  flows  with  extensive  collapse  features 
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[9,53].  Extensive  collapse  within  the  ASC  also  sup¬ 
ports  the  high  subsidence  rates  predicted  by  the 
bimodal  lava  emplacement  model  near  the  rise  crest. 

5-  Conclusions 

A  stochastic  model  of  the  emplacement  of  dikes 
and  lava  flows  generates  an  upper  oceanic  crustal 
structure  similar  to  that  observed  in  ophiolites,  in 
ocean  drill  holes  and  in  seismic  data  from  the  East 
Pacific  Rise.  An  emplacement  model  with  a  narrow 
intrusion  zone  (95%  of  the  dikes  are  intruded  within 
a  zone  40  m  wide,  centered  on  the  ridge  axis)  and  a 
bimodal  emplacement  of  lavas  can  explain  the  varia¬ 
tion  in  thickness  of  layer  2A  observed  seismically  at 
the  EPR  and  the  sharpness  of  the  extrusive-dike 
transition  region  observed  in  outcrop  in  ophiolites 
and  on  the  seafloor.  A  thin,  on-axis  extrusive  section 
is  rapidly  built  up  by  accumulating  a  sequence  of 
short  flows  erupted,  and  probably  confined,  within 
the  ASC.  The  progressive  thickening  of  the  extrusive 
layer  inferred  from  the  deepening  of  the  seismic 
2A/2B  boundary  over  1-3  km  from  the  ridge  axis  is 
caused  by  longer  flows  that  spill  out  of  the  ASC  and 
down  the  flanks  of  the  axial  high  or  erupt  off-axis 
(Fig.  6).  In  the  resulting  lava  pile  the  upper  part  of 
the  extrusive  section  is  composed  of  flows  deposited 
off  the  rise  crest  while  the  lower  portion  is  flows 
deposited  on-axis.  This  stochastic  model  for  the  con¬ 
struction  of  the  upper  crustal  section  is  consistent 
with  the  velocity  structure  of  the  upper  crust  (pro¬ 
ducing  a  near-surface  low  velocity  layer  underlain  by 
a  steep  velocity  gradient),  paleomagnetic  inclinations 
in  drill  cores  (a  bimodal  dip  structure),  magnetic 
transition  widths  ( '^  2  km),  and  age  offsets  of 
seafloor  lavas  relative  to  the  underlying  crust  (  ^  30- 
70  kyr).  Extensive  off-axis  eruptions  are  not  neces¬ 
sarily  required  to  explain  these  age  discrepancies. 
The  similarity  in  seismic  structure  between  the 
northern  and  southern  EPR  (full  spreading  rates  of 
110  and  150  mm/yr,  respectively)  indicate  that  the 
bimodal  lava  emplacement  model  is  applicable  to 
fast  spreading  ridges  over  a  range  of  spreading  rates. 
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Review  of  Results 


As  melt  rises  from  the  mantle  at  mid-ocean  ridges,  it  freezes  and  erupts  to  build  the 
oceanic  crust.  The  thermal  regime  is  important  for  determining  the  depth  at  which  magma 
freezes  within  the  crust  but,  through  its  control  on  the  rheology  of  the  crust,  also 
determines  the  morphology  of  the  axis  of  the  mid-ocean  ridge.  This  Thesis  underlines  the 
importance  of  the  crustal  thermal  structure  in  controlling  local  variations  in  magma  sill 
depth  and  ridge  morphology,  and  of  eruption  dynamics  in  the  emplacement  of  the 
uppermost  volcanic  layer  of  the  crust. 

Analysis  of  gravity  over  the  Juan  de  Fuca  and  Gorda  ridges  in  Chapter  2  shows  that 
even  small  differences  in  crustal  thickness  (300-700  m)  and  mantle  temperature  (10-15°C) 
are  sufficient  to  produce  major  changes  in  lithospheric  strength  and  axial  morphology.  The 
presence  of  the  Cobb  thermal  anomaly  beneath  Axial  Seamount  may  cause  the  systematic 
differences  in  axial  morphology,  crustal  thickness,  mantle  temperature,  and  lithospheric 
strength  observed  between  the  Juan  de  Fuca  and  Gorda  ridges.  The  hotter  mantle  under 
the  Juan  de  Fuca  Ridge  results  in  a  greater  amount  of  decompression  melting  and  the 
formation  of  a  thicker  crust  and  a  thinner,  weaker  lithospheric  than  along  the  Gorda  Ridge. 
These  results  are  consistent  with  the  predicted  sensitivity  of  recent  thermo-mechanical 
models  for  rift  valley  formation  to  small  changes  in  crustal  thickness  and  mantle 
temperature  at  intermediate  spreading  rates  [Phipps  Morgan  and  Chen,  1993b]. 

However,  recent  results  from  other  intermediate  spreading  ridges  show  that  the 
transition  from  axial  valley  to  axial  high  is  more  complicated  than  these  observations 
suggest.  Canales  et  al.  [1996]  show  that  along  the  Galapagos  ridge  (spreading  rate 
somewhat  faster  than  along  the  Juan  de  Fuca  ridge:  31  vs.  28  mm/yr  half-rate)  the 
transition  occurs  over  a  distance  of  200-300  km  and  is  accompanied  by  a  change  in  mantle 
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temperature  of  ~20°C  and  >1-2  km  in  crustal  thickness  [Ito  and  Lin,  1995],  On  the  other 
hand,  Ma  and  Cochran  [1996]  show  that  along  the  intermediate  spreading  southeast  Indian 
Ridge  there  are  three  distinct  ridge  morphologies  not  two;  a  rift  valley,  transitional 
morphology,  and  axial  high.  The  transitions  between  these  different  morphological  types 
occur  over  short  lengths  of  the  ridge  axis.  Thus,  though  intermediate  spreading  ridges 
display  morphologies  ranging  form  axial  valley  to  axial  high,  the  processes  controlling 
axial  topography  may  be  more  complicated  than  accounted  for  by  recent  thermo-mechanical 
models  for  rift  valley  formation. 

Chapter  3  presents  a  quantitative  examination  of  the  correlation  between  along-axis 
variations  in  crustal  structure  and  other  indicators  of  magma  supply  on  a  regional  scale. 
There  are  several  scales  of  spatial  and  temporal  variability  that  modify  the  simple 
correlations  predicted  by  the  current  magma  supply  model  for  fast  spreading  ridges.  Good 
correlations  exist  along  the  SEPR  between  axial  depth,  ridge  cross-sectional  area,  mantle 
Bouguer  anomaly,  and  the  MgO  wt%  of  basalts  recovered  from  the  rise  axis.  These 
correlations  suggest  along-axis  changes  in  crustal  thickness  and  temperature  consistent  with 
variations  in  magma  supply  on  time  scales  of  -100  Kyr.  However,  we  show  that  the 
magma  sill  properties  (depth  and  width  of  the  melt  lens),  seismic  layer  2A  thickness,  and 
the  intensity  of  hydrothermal  venting  are  poorly  correlated  with  regional  variations  in  ridge 
depth  and  cross-sectional  area.  We  suggest  that  the  emplacement  geometry  (width  of  the 
intrusion  zone  and  flow  lengths),  not  magma  supply,  controls  extrusive  layer  (seismic  2A) 
thickness.  The  depth  of  the  magma  sill  is  very  sensitive  to  the  local  thermal  structure. 
Magma  lens  properties  and  hydrothermal  activity  are  closely  linked  to  spreading  events 
(dike  intrusion,  eruptions,  faulting)  which  occur  on  much  shorter  time  scales  (ca.  10-100 
yr)  than  the  longer-term  variations  in  magma  supply  reflected  in  along-axis  changes  in  ridge 
morphology. 
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The  observed  variation  in  magma  sill  depth  is  compared  to  different  models 
proposed  for  the  accumulation  of  magma  in  the  mid-crust.  In  Chapter  4,  it  is  shown  that 
the  seismically  constrained  depths  of  ridge  crest  magma  chambers  (>l-2  km)  are 
considerably  deeper  than  the  level  of  neutral  buoyancy  (100-400  m),  suggesting  that  other 
factors  must  control  the  collection  of  melt  in  these  reservoirs.  The  apparent  inverse 
relationship  between  magma  sill  depth  and  spreading  rate  [Purdy  etal,  1992]  suggests  that 
a  thermally  controlled  permeability  boundary,  such  as  the  solidus  horizon,  controls  the 
depth  at  which  magma  ponds  in  crustal  magma  chambers  beneath  mid-ocean  ridges  [c.f. 
Phipps  Morgan  and  Chen,  1993a;  Natland  and  Dick,  1996].  In  Chapter  3,  a  more  refined 
approach,  which  takes  into  account  new  evidence  that  the  magma  sill  is  not  fully  molten 
[Hussenoeder  et  ai,  1996],  reaches  the  same  conclusion.  Observation  of  large  local 
variations  in  magma  sill  depth  (800-2000  m  below  the  seafloor)  suggest  that  on  shorter 
time  scales  the  local  thermal  structure  of  the  crust  is  strongly  modified.  The  advection  of 
heat  by  the  ascent  of  melt  from  the  mantle  is  relatively  continuous.  However,  the  removal 
of  heat  from  the  crust  is  highly  variable  and  depends  on  the  variability  in  hydrothermal 
cooling,  on  intermittent  eruptive  and  intrusive  episodes  that  drain  the  magma  sill,  and  on 
tectonic  events  that  alter  the  crustal  permeability  structure. 

A  stochastic  model  for  the  emplacement  of  dikes  and  lava  flows  is  presented  in 
Chapter  5.  This  model  generates  an  upper  oceanic  crustal  structure  similar  to  that  observed 
in  seismic  data  from  the  East  Pacific  Rise  (EPR),  in  ocean  drill  holes,  and  in  ophiolites.  In 
order  to  satisfy  these  geological  and  geophysical  constraints,  a  bimodal  distribution  of  lava 
flows  is  required.  Short  flows,  confined  within  the  axial  summit  caldera,  build  up 
approximately  half  the  extrusive  volume.  Occasional  flows  spill  out  of  the  ASC,  or  erupt 
outside  the  ASC,  and  pond  at  a  considerable  distance  off-axis  (up  to  2  km)  to  build  up  the 
remainder  of  the  extrusive  section.  Most  dikes  are  intruded  within  a  narrow  zone  at  the 
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ridge  axis.  The  upper  part  of  the  final  extrusive  section  will  be  dominated  by  the  off-axis 
flows,  while  the  lower  portions  will  be  primarily  composed  of  short  flows  erupted  within 
the  ASC.  Predicted  magnetic  transition  widths  are  similar  to  those  measured  in  deep-tow 
studies  [Macdonald  et  ai,  1980].  The  ages  of  seafloor  lavas  in  this  model  are  younger 
than  the  tectonic  spreading  model  ages  by  30-70  Kyr,  in  agreement  with  anomalously 
young  lava  ages  obtained  from  radioisotope  disequilibrium  dating  of  seafloor  basalts  from 
the  EPR  [Goldstein  etai,  1994]. 

Further  Work 

Numerous  questions  about  the  volcanic  processes  of  crustal  construction  and  the 
rheology  controlling  axial  morphology  and  magma  migration  remain  to  be  addressed.  The 
crustal  magmatic  system  may  represent  the  transition  between  semi-continuous  melt 
percolation  through  the  mantle  and  lower  crust  to  intermittent  dike  emplacement  in  the  crust 
and  eruptions.  How  does  the  composite  magma  chamber  work  exactly?  What  is  the 
relationship  between  the  more  steady-state  low-velocity  zone  and  the  overlying,  rapidly 
varying  melt  sill?  Since  there  are  not  significant  amounts  of  melt  within  the  crust,  are  there 
magma  reservoirs  within  the  mantle?  How  is  the  melt  delivered  to  the  crust  from  the 
mantle?  Detailed  seismic  and  electromagnetic  studies  would  improve  our  understanding  of 
the  structure  of  the  lower  crust  and  upper  mantle  at  the  ridge  crest.  Also,  there  is  still  no 
definitive  explanation  for  how  the  composite  volcanic  system  constructs  the  oceanic  crust, 
in  particular,  the  layered  gabbro  sequence. 

The  spatial  and  temporal  relationships  of  the  processes  controlling  the  magmatic 
system  are  still  not  well  understood.  Repeating  a  seismic  survey  in  a  well  studied  area 
(e.g.  9°N  EPR)  would  give  great  insight  into  the  time  dependence  of  the  magmatic  system, 
particularly  the  temporal  and  spatial  variations  in  the  properties  (depth,  width,  crystal 


152 


content  and  thickness)  of  the  magma  sill.  Greater  insight  into  the  relationship  of 
hydrothermal  venting  to  magma  injection  events,  faulting,  and  recent  eruptions,  will  better 
constrain  the  importance  of  hydrothermal  cooling  versus  intermittent  eruptions  in 
controlling  the  local  thermal  structure  of  the  upper  crust.  A  detailed  study  of  the  temporal 
variability  of  hydrothermal  vents  and  their  longevity  will  improve  our  understanding  of  the 
processes  controlling  their  formation  and  the  time  scales  over  which  they  are  active. 

How  does  the  thermal  structure  of  the  axial  lithosphere  control  its  strength  and  how 
closely  is  the  morphology  of  the  ridge  related  to  the  presence  or  absence  of  a  melt  sill. 
Further  work  using  more  realistic  rheological  models  of  the  axial  lithosphere  are  necessary 
for  an  improved  understanding  of  the  controls  on  axial  morphology. 

For  an  improved  understanding  of  the  construction  of  the  extrusive  section  of  the 
oceanic  crust  more  constraints  on  the  size  and  composition  of  lava  flows  and  the  location  of 
their  eruption  and  temporal  frequency  is  required.  A  detailed  study  of  the  ridge  crest  to 
date  and  sample  lava  flows  in  order  to  determine  the  typical  frequency  and  location  of 
eruptions,  the  evolution  of  lava  composition  during  an  eruption,  and  the  lateral  extent  and 
volume  of  lava  flows.  Mapping  of  the  extrusive  section  in  ophiolites  is  another  way  of 
testing  the  model  for  the  emplacement  of  the  volcanic  section  proposed  in  Chapter  5. 
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